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Abstract 

This thesis describes a new a measurement of the tau lepton polarization 
and its forward-backward asymmetry at the Z° resonance using the OPAL 
detector. This measurement is based on analyses of the r— >ePei^r , T^iiv^Vr 
and T— >7r(K)i/^ decays from a sample of 30663 e'^e" — >t+t~ events collected 
in the polar angle range of | cos6'| < 0.68 during the 1990-1992 data taking 
period. Taking then the Standard Model with the V-A structure of the tau 
lepton decay, we measure the average r polarization to be 

(P^) = (-13.5 ±2.9 ±2.2)% 

and the r polarization forward-backward asymmetry to be 

K^B ^ (-11.0 ±3.5 ±0.5)% 

where the first error is statistical and the second systematic. Combining these 
figures with the OPAL (P^) and A^^f measured with the r^pur channel we 
get an average r polarization, 

(P,) = (-14.9 ± 1.9 ± 1.3)% 

and for the asymmetry in the polarization, 

Ajf = (-8.9 ± 2.2 ± 0.9)% 

These results are consistent with lepton universality. Combining the two 
results we obtain for the electroweak mixing angle the value 

sin^ e^= 0.2331 ± 0.0023 



within the context of the Standard Model, where the error includes both 
statistical and systematic uncertainties. 
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Chapter 1 
Introduction 



One of the central phenomena characterizing weak interactions is the non- 
conservation of parity. This effect, originally established for weak charged- 
current interactions, is also embedded in the Standard Model [|I], ||, |^ to exist 
in neutral-current interactions, resulting in different TP couplings to left- 
handed and right-handed fermions. Consequently, the Z° boson produced by 
e"'"e~ annihilation is expected to be polarized because of its different couphngs 
to the incoming left-handed and right-handed electrons. Similarly, fermions 
produced in T? decay are expected to have a degree of polarization depending 
on their coupling constants. Another consequent of this Z° polarization is 
a forward-backward asymmetry in the polarization of the outgoing fermions. 
Hence, some of the best Standard Model tests in the annihilation of e^e~ at 
the Z° pole provided by the following asymmetries measurements: 

• The Forward Backward asymmetry, Kp^-, which can be measured 
in all the e+e^ — > // channels. 

• The Left Right Asymmetry, Aj^^j, which can be measured in the 
annihilation of longitudinal polarized e^e~ beams (at the moment avail- 
able only at SLAC's SLD experiment). 
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• The Final Lepton Polarization Asymmetry, (Pf) which can be 
measured only in the r decay ((-Pr))- 

• The Forward Backward Asymmetry in the Lepton Polariza- 
tion, A^^f , which again can be measured only in the r decay . 

The (Pr) and the A^^f asymmetries can be studied in the process e+e" -^t^t^ 
using the energy distribution of the r decay products in the laboratory 
frame @]. 

The mixing angle between the electromagnetic and the weak interaction, 
siia^9\Y, plays a central role in the theory of the Standard Model. Being so 
it is measured by the LEP experiments in several processes and by various 
methods. The Standard Model gives predictions for ApB, {Pr) and A^^f as 
a functions of s = Ecm^ in terms of the Z° parameters (mass and width) 
and its vector (vi) and axial- vector (ai) couplings to the electron and the 
tau. Using the improved Born approximation which accounts for most 
weak radiative corrections, our measurements of (Pr) and A^^f provide a 
test of e-r universality in the neutral current which is independent of lepton 
universality tests performed by studying the line shapes and the forward- 
backward asymmetries of the Bhabha, /i-pair and r-pair cross-sections @. 

Our results supercede OPAL's first r polarization measurements which 
was based solely on the 1990 data. The other LEP collaborations have also 
reported their measurements of the r polarization H]-[|l3l forward-backward 



asymmetry |TTl-|T3[, based on their 1990 and 1991 data. 

The present thesis describes a measurement of the r polarization, (Pr), 

and its forward-backward asymmetry, Ap^f , using the r-pair events collected 

with the OPAL detector at LEP during the period 1990-1992. It is based on 

a sample of 30663 e"*"e~ —^t^t^ events which were detected within the polar 

angle range of[| | cos 6^1 < 0.68. Most of these events (91.5%) were measured 

"'^The coordinate system is defined with z along the e~ beam direction, and (f) being 
the polar and azimuthal angles, respectively. 
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on the peak and the remainder at center-of-mass energies (Ecm) of 1, 2, 
and 3 GeV below and above the peak of Z° resonance. The decay channels 
r^ez/gZ/^, T^^v^Vj. and r— >7r(K)z/T- are used. 



These results combined with the analysis of the r^pv^ decay channel [|I4 



15| were presented as an OPAL publication JT^ |T6|. Therefore for the sake 



of completeness we present here also the combined results. 

The analysis is based on an event-by-event Maximum Likelihood fit to 
the theoretical energy distributions of the r decay product. The energy is 
corrected for radiative effects and detector response, in which all correlations 
between the two tau decays are taken into account. We apply our 'global fit' 
method to the three r decay channels from which the polarization is extracted 
using distributions in simple observables. This method has the advantage of 
taking explicitly into account the correlations between polarization observ- 
ables of the two taus introduced by the selection and identification criteria. 
This is particularly important for the leptonic channels where requirements 
are made on the whole event in order to suppress backgrounds from Bhabha 
and mu-pair events. In addition, the method correctly extracts information 
using the tau-tau spin correlations in those events where both taus have 
identified decays. 

The thesis is organized as follows: A description of the LEP accelerator 
and the OPAL detector are presented in the following two Chapters. A 
brief description of the Standard Model and a discussion of the polarization 
formalism, are given in the Chapter |. It includes the definition of the various 
observables used in this analysis and the relations between them. 

Chapter ^ presents the data and Monte Carlo samples used in the analy- 
sis, the event selection and r-decay identification. Chapter ^ describes the fit 
procedure and all the corrections used in our analysis, starting with radiative 
effects (Sect. |0|) , detector resolution (Sect. |6.3|) , r-pair selection efficiency 
(Sect. |6.4| ), r decay identification efficiency (Sect. |6.5| ), and completing this 
list with background from other r decay channels (Sect. |6.6.1| ) and from non- 
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r events (Sect. |6.6.2| ). Each correction is described along with the associated 
systematic studies. The results, including a summary of the systematic er- 
rors, and a presentation of the consistency checks done on the analysis, are 
included in Chapter |^. Finally the analysis is summarized in Chapter |[ 



Chapter 2 

The LEP Accelerator 



2.1 General Features 

The electrons and positrons which annihilate in the experiment are produced 
in the Large Electron Positron ring LEP, the current worlds largest e^e~ 
collider [^-[^. At the present stage LEP-I can provide colliding beams of 
energies upto ~ 47 GeV. With the machine upgrade planed for the next years 
LEP-II will accelerate the beam to ~ 100 GeV. The circumference of LEP is 
26.66 km, partly in France and partly in Switzerland. The collider tunnel lies 
between about 43 m and 140 m below the surface. As injectors, LEP uses 
existing CERN machines. The LEP ring, has four interaction regions and was 
designed to reach a maximum luminosity of about L ~ 1.7x 10'^^ cm~^s~^ at a 
center of mass energy around Ecm = 90 GeV. The storage ring has eight long 
straight sections, four of which are used for the particle physics experiments: 
L3, ALEPH, OPAL, and DELPHI. Fig. U shows a schematic layout of 
the LEP installations. The construction of LEP and its four detectors were 
completed in July 1989, and the first event was seen in OPAL on the IS*'^ 
of August 1989. 



5 
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2.1.1 Structure 

The LEP storage ring is constructed out of four basic types of components, 
namely: 

• The arc. The arc bends the beam in a circle and focuses it at the 
same time. These regions consists of regular arrays of quadrupoles and 
bending magnets, and most of the beam parameters are determined by 
its optics. 

• Dispersion suppressor. Particles with different energies have different 
orbits in the ring. This can increase the beam size in the interaction 
point which reduces the luminosity. In order to avoid the luminosity 
reduction the dispersion suppressor are introduced in the RF section 
and in the low beta insertion. 

• RF cavities. The function of it is to replace the energy lost by the par- 
ticles through synchrotron radiation and other effects, and to provide 
a phase focussing. 

• Interaction region. In theses regions the two beams collide. In order 
to get a high rate of interactions the beams are strongly focussed at 
these collision points. This is obtained by the low beta quadrupoles 
and a set of additional electro-magnetic lenses which optically match 
the insertion to the rest of the machine. 

The focussing and bending property of the beam optics leads to two effects. 
First, a particle with nominal energy but with an uncertainty in position Ax 
and/or in angle Aip, will get focussed toward the nominal orbit and oscillate 
around it. This oscillation is called the betatron oscillation and it has the 
form 

x{e) = a^JJ{t)cos{(t){t) + 0o). (2.1) 
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Figure 2.1: A schematic layout of the LEP storage ring. Point 1 is the 
injection point, while the even points are the four detectors: L3, ALEPH, 
OPAL, and DELPHI. The circumference of this ring is 26.6 km and it lies 
between 43 m and 140 m under the ground. 
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where £ is the longitudinal position in the ring. The square root of the 
periodic function f3{i) gives the envelope of the betatron oscillation over 
many turns. The P function value is large in a focussing quadrupole and 
small in a defocussing one. 

The additional properties of the focussing and bending is to restrict the 
particles to an energy region AE around the nominal energy . On one hand 
such particles with an excess energy do bend less in the dipole magnets, but 
on the other hand they go off center through the quadrupolcs and therefore 
receive an extra bending. The result is, that the off-energy closed orbit 
particles have a horizontal displacement Axh 



where is called the periodic dispersion, and is a function of the longitu- 
dinal position in the ring, Dj. — Dx{i). The dispersion is typically large at 
the horizontally focussing quadrupoles and typically small at the defocussing 
ones. In LEP it is arranged so that the dispersion vanishes in all the straight 
sections by using proper transition to the arc, namely, the dispersion suppres- 
sors. In the arc section, where there is a finite dispersion, the deviation of the 
circumference is proportional to the energy deviation. The RF acceleration 
works in such a periodic function that a particle with an excess energy, which 
has a longer path length than a nominal energy one, will arrive later to the 
RF voltage area. In its time of arrival the voltage is smaller, and therefore 
it will suffer an energy loss. All this leads to oscillation around a Gaussian 
equilibrium of the LEP particles in their energy, longitudinal position, and 
transverse angle and position. 



Axh = 



AE 



(2.2) 



E 
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2.2 Luminosity 

The reaction rate is proportional to the luminosity, L, and the reaction cross 
section, a. The luminosity is a machine parameter, given by: 

r N'kfo 



(2.3) 



where N is the number of particles per bunch, k the number of bunches, /o 
the revolution frequency and ax and ay the RMS values of the transverse 
radii at the collision points. 

However there is an important limitation imposed by the beam-beam ef- 
fect, due to the electromagnetic forces between the electrons and the positrons 
in the two crossing beams. This force has a linear part which leads to a change 
of the betatron frequency called the beam-beam tune shift. The particles os- 
cillating with small amplitudes compared to a^ and ay at the interaction 
points, suffer from the largest change in their tune. This shift is given by: 

2'Kyax{ax + ai) 

where the index i stands for x or y, rg is the classical electron radius, (5x and 
(5y are the transverse beta function at the crossing points. It turns out that 
for most storage rings experiments, that the upper limit on is between 
0.03 and 0.06. To optimize the luminosity estimation we can choose a tune 
shift as large as possible and this leaves us with 

Therefore, the experimental collision points are arranged in low (3 insertions 
where fi^ and j3y are much smaller than in the rest of LEP, namely that a^ 
and ay become a fraction of a millimeter. 

The cross section cr is a property of the reaction itself, and one often 
compares it to the QED process e"'"e~ ^ A*"*"/^"- The cross section to this 
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process decreases with 1/s. It value is ~ 10^^^ cm^ when Ebeam = Mz/2, 
which means that with a typical LEP luminosity of 10^^ cm~^s~^ one will 
obtain a Z° counting rate of ~ 0.36 hour~^. On the mass the main 
annihilation channel is via the formation of Z boson, e^e~ — > X, with a 
cross section higher by a factor of 5 x 10^ than the propagation via a 7. Thus 
we expect on the peak about 1800 events per hour. 

A typical luminosity for LEP during 1990 was about 0.4 x 10^^ cm^^s"^ at 
the beginning of a fill. This leads to an average of one multihadronic eventQ 
every 8 sec, forward Bhabha event every 5 sec, and decay into charge 
leptons every 50 sec. The TP production rate has further increased in the 
following years as can be seen in Table |2.1| which gives the OPAL integrated 
luminosity in the years 1990-1992. In 1993 the four LEP experiments to- 
gether saw 3 million TP s, making a total of some 8 million since LEP began 
operating. 



Year 


Int. luminosity 


No. of MH 


1990 
1991 
1992 


5.637 

12.236 ph-^ 
21.6319 ph-^ 


147425 
353324 
767156 



Table 2.1: The integrated luminosity and the number of MH events measured 
with the OPAL detector in the years 1990-1992. 



2.2.1 Energy Loss Due to Synchrotron Radiation 

The energy and the number of particles in a bunch are both limited by the 
synchrotron radiation. Electrons and positrons circulating in an orbit with 
radius p lose energy due to synchrotron radiation. This loss per turn is given 

^ TP decays to quark- antiquark pairs produce final states with many hadrons, refers in 
the following as MH 
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by: 

Us = ^i:^^!^ (2.6) 
3 p 

where p is the bending radius of the arc which is 3096 m in LEP. From this 
follows that the energy loss is 140 MeV and 2330 MeV for LEP operated 
at energies of 47 GeV and 95 GeV respectively. This loss has to be replaced 
by the RF system. The RF requirements are dominated by the losses in the 
cavities, which in storage rings increase with the beam energy. Therefore, 
synchrotron radiation causes sharp upper limits on both current and energy. 

2.3 The Injection Chain 



The scheme for injection into LEP is shown in Fig. |2.2| . An electron linac 
(LIN ear Accelerator) which is a high current, 200 MeV machine producing 
an output current of 2.5 A for electron-to-positron conversion in a tungsten 
target. From there a positron current of 12 mA is produced in a form suitable 
for subsequent acceleration in another linac which is used to accelerate these 
positrons up to an energy of 600 MeV . The first linac is also used to provide 
electrons for LEP by detuning its electron gun to produce 110 mA at the 
linac output for acceleration to 600 MeV in the second linac. The next 
stage in this chain is transferring the electrons or positrons, as appropriate, 
into an electron-positron accumulator ring (EPA) where they are stored and 
accumulated before injection into the Proton Synchrotron (PS). Acceleration 
to 3.5 GeV takes place in the PS, followed by transfer to the Super Proton 
Synchrotron (SPS) for acceleration to 20 GeV and finally injection into LEP. 
The 20 GeV electrons and positrons are stacked in LEP in eight (four until 
the middle of 1992) equally spaced bunches for each sign of particles. This 
whole process is repeated successively until the required currents are reached. 
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Figure 2.2: Injection into LEP. The LIL output are 600 MeV electrons and 
positrons. The EPA transfer the e~ and as appropriate into the 3.5 GeV 
PS ring. From the PS the e~ and e+ are injected to the SPS ring where they 
are accelerated to 20 GeV, and finally injected into LEP in eight equally 
spaced bunches of electrons and positrons. 
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2.3.1 Magnet System 

For bending, magnets at LEP have low fields between 0.05 T and 0.1 T. The 
accelerated e+ and e~ trajectories are bent by 3304 dipole cores installed 
over the whole 20 km of the arcs. At each end of an arc, there is a special 
weak-field dipole where the magnetic field has about ^ of its normal value, 
to shield the experiments from synchrotron photons emitted in the main 
bending dipoles. At regular intervals 488 quadrupoles magnets are installed, 
to focus the beam around the arc. These quadrupoles, 1.6 m magnetic length, 
have a maximum gradient of 9.7 Tm~^. In the long straight sections, there are 
288 stronger quadrupoles 2 m long with a maximum gradient of 10.9 Tm~^. 
On each side of the of the four experiments a stronger 2 m superconducting 
quadrupole is installed providing a field gradient of 36 Tm~^ to squeeze 
the transverse beam dimension as much as possible in order to increase the 
luminosity. 

LEP, as other colliders, suffers from chromaticity, which is a phenomenon 
of particles with off momentum values which are focused differently in the 
quadrupoles. In order to correct the machine's chromaticity, 504 sextupoles 
magnets of two different lengths 0.4 m and 0.76 m, with a maximum strength 
B" — 180 Tmr'^, were installed near the quadrupoles in the arc. This produce 
a quadratic field variation across the apparatus. The closest orbit distortions 
resulting from the remaining field, and alignment errors of the dipoles and 
quadrupoles, are corrected by means of 616 x 0.5 m long orbit correctors, 
which have a maximum field between 0.04 T and 0.08 T and are also placed 
near the quadrupoles. Other magnets exist in smaller numbers, mainly near 
the interaction and the injection points. Most of the magnet coils are water- 
cooled to remove the electrical power dissipated. 
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2.3.2 RF System 

There are eight e" bunches turning counter-clockwise, and eight e+ bunches 
turning clockwise symmetrically spaced in the LEP ring (four x four at until 
1992). The beam intensity is built up by multiple injections from the SPS 
into the same eight bunches. Once the required intensity is reached, the 
two beams are simultaneously accelerated to their final energy by the RF 
system and kept in their respective trajectories throughout the run. The 
acceleration structure consists of 128 coupled cavity units, each containing 
a five-cell accelerating cavity of 2.12 m active length, which can give a peak 
gain of 3.1 MeV to an electron or positron, i.e. the whole structure gives a 
maximum accelerating gradient of 1.47 MVm~^. The nominal value of peak 
RF voltage per revolution is 400 MV , and to drive the cavities at their full 
rating, a total of 16 MW RF power is required. 

2.3.3 Vacuum System 

In order to limit the particle loss rate due to beam-gas interactions, LEP 
requires a vacuum of 10~^ Torr when the beam is circulating, which is equiv- 
alent to the requirement of 10~^^ Torr without beam. The main difficulty of 
the vacuum system, is the synchrotron radiation which causes an outgassing 
of gas molecules, and can in this way raise the pressure by orders of magni- 
tude during the beginning of a run. The outgassing is highly reduced with 
the increasing of the beam. This effect, called the beam cleaning, allows 
the installation of a smaller pumping system. The huge 26.66 km vacuum 
chamber of LEP is divided into sectors of up to 474 m length. Each sec- 
tor consists of two pumping systems, one for obtaining the starting vacuum 
of ~ 10~^ Torr, and the second for reaching and maintaining the needed 
ultra-high vacuum. The synchrotron radiation also causes the heating of the 
vacuum chambers. Therefore the chambers are made of aluminum which 
has a good thermal conductivity and these chambers are also surrounded 



CHAPTER 2. THE LEP ACCELERATOR 



15 



by water-cooling channels. The heat which escapes the aluminum can cause 
severe radiation problems and therefore the aluminum is covered with lead 
of up to 8 mm in thickness. 

2.3.4 Electrostatic Separators 

During the injection and acceleration, it is necessary to separate the e"*" and 
the e~ beams. To achieve this, each of the eight collision points is equipped 
with four electrostatic separators. These separators can force a total distance 
of 0.49 mm between the and the e~ bunches at 55 GeV in the low beta 
insertions, and a 3.25 mm at the high beta insertions. The system is designed 
so that at higher LEP energies, sufficient separation can be obtained by 
adding two more separators at each experimental collision region. After the 
needed energy is reached, the beams are brought into collision in the four 
experimental points, while being kept separated elsewhere in LEP. 

2.3.5 LEP Parameters 



Table p.3.5| gives some general details on the LEP accelerator, and its mag- 
netic and RF system. In one column are the LEP-I parameters as they are 
at present, and for a comparison, the LEP-II designed parameters are also 
given. 
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units 


phase 1 


phase 2 


Circumference of the machine 


km 


26.658 


26.658 


Inner diameter of the tunnel 


m 


3.80 


3.80 


Radius of curvature in one dipole 


m 


3096 


3096 


Injection energy 


GeV 


20 


20 


Number of e"*" and e~ bunches 




8(4) 


4 


Bunch length 


m 


0.013 to 0.04 


0.012 to 0.04 


Maximum luminosity 


—1 —1 
cm s 


1.7 X 10^1 


7 X 1031 


Circulating current per beam 


mA 


3.0 


3.0 


Magnetic System 








Dipoles magnetic field at injections 


T 


0.0215 


0.0215 


Maximum dipoles magnetic field 


T 


0.0590 


0.1075 


Number of iron dipoles 




3304 


3304 


Number of dipoles with weak fields 




64 


64 


Number of dipoles in the injection region 




24 


24 


Number of quadrupoles 




816 


800 to 816 


Number of focussing sextupoles 




248 


248 


Number of defocussing sextupoles 




256 


256 


Radio-Frequency System 








Number of equipped straight sections 




2 


4 


Number of klystrons 




16 


20 


RF frequency 


MHz 


352.21 


352.21 


Necessary circumferential voltage 


MV/turn 


360 


3224 


Synchrotron energy loss per particle 


MeV j turn 


260 


2830 


Total synchrotron radiation power 


MW 


1.6 


17.1 



Table 2.2: Design LEP parameters. The third column shows LEP-I param- 
eters, the LEP-II designed parameters are given in the fourth column. 



Chapter 3 

The OPAL Detector 



OPAL (Omni-Purpose Apparatus for LEP) is a multipurpose apparatus 
designed to reconstruct efficiently and identify all types of e"'"e~ events. Fur- 
thermore, OPAL was designed so that it will has a good acceptance for Z° 
decays over a solid angle of nearly 47r. A general layout of the detector is 
shown in Fig. indicating the location and relative size of the various 
components. 

The main features of the detector are: 

• Tracking of charged particles, performed by the central detector, and 
providing measurements of the particles direction and momentum, their 
identification by dE/dx and reconstruction of primary and secondary 
vertices at and near the interaction region. 

• Identification of photons and electrons by the Electromagnetic Calorime- 
ter and measurement of their energy. 

• Measurement of hadronic energy by total absorption using the magnet 
yoke instrumented Hadron Calorimeter. 

• Identification of muons by measurement of their position and direction 
within and behind the hadron absorber. 

17 
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Figure 3.1: A schematic view of the entire OPAL detector. The incident 
beams travel along the central beam pipe and collide at the center of the de- 
tector. In each direction particles traverse the central detector, followed by 
the electromagnetic and hadronic calorimeters and finally the muon cham- 
bers. Particles produced at small angles are detected by the forward detec- 
tors. 
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• Measurement of absolute machine luminosity using Bhabha scattering 
events in the very forward direction with respect to the beam line. 

The OPAL has a 3D cartesian coordinate system whose origin is at the 
nominal interaction point. It has the 2;-axis along the nominal electron beam 
direction (this is anticlockwise around LEP when viewed from above), the 
X-axis horizontal and directed towards the center of LEP, and the y-axis 
normal to the z-x plane. 



Full details of the OPAL detector can be found in Ref. p3|, Only a 
very brief introduction is given in the successive sections describing the parts 
of the detector most important to this analysis. In particular we concen- 
trate on the methods the sub-detectors information are utilized in the event 
selection and the physics analysis. 



3.1 The Central Detector 

The Central Detector (CD) consists of a Silicon Microvertex detector and 
three drift chamber devices, the vertex detector, jet chamber and surrounding 
Z-chambers situated inside a pressure vessel of 4 bar. The central detector 
is inside a solenoid supplying a uniform axial magnetic field of 0.435 T. 
Originally and until 1991 there was no Silicon detector and the inner wall of 
the pressure vessel at 7.8 cm radius formed the beam pipe. This beam pipe 
consisted of 0.13 cm thick carbon fibre with a 100 /im aluminium inner lining. 
In 1991 a second beam pipe was added at a radius of 5.35 cm, consisting of 
0.11 cm thick Beryllium, and the Silicon detector inserted between the pipes. 



3.1.1 Jet Chamber 



The jet chamber (CJ) is a cylindrical drift chamber of length 400 cm with 
an outer radius of 185 cm and inner of 25 cm. The chamber consists of 24 
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identical sectors each containing a sense wire plane of 159 wires strung par- 
allel to the beam direction. The end planes are conical and can be described 
by \z\ = 147 + 0.268 x R. 

The coordinates of wire hits in the r — plane are determined from a mea- 
surement of drift time. The z coordinate is measured using a charge division 
technique and by summing the charges received at each end of a wire allows 
the energy loss, dE/dx, to be calculated. 

3.1.2 Z- Chambers 

The Z-chambers (CZ) provide a precise measurement of the z coordinate of 
tracks as they leave the jet chamber. They consist of a layer of 24 drift 
chambers 400 cm long, 50 cm wide and 5.9 cm thick covering 94% of the 
azimuthal angle at the polar angle range of | cos^|< 0.72. Each chamber is 
divided in the z axis into 8 cells of 50 x 50 cm^, each of which has 6 sense 
wires spaced at 0.4 cm. 

3.1.3 Vertex Detector 

The vertex detector (CV) is a high precision cyhndrical jet drift chamber. 
It is 100 cm long with a radius of 23.5 cm and consists of two layers of 36 
sectors each. The inner layer contains the axial sectors, each containing a 
plane of 12 sense wires strung parallel to the beam direction. The wires range 
radially from 10.3 to 16.2 cm with a spacing of 0.583 cm. The outer layer 
contains the stereo sectors each containing a plane of 6 sense wires inclined 
at a stereo angle of ~4°. The stereo wires lie between the radii 18.8 and 
21.3 cm with a spacing of 0.5 cm. 

A precise measurement of the drift time on to the axial sector sense wires 
allows the r — position to be calculated. Measuring the time difference 
between signals at either end of the sense wires allows a fast but relatively 
coarse z coordinate measurement that is used by the OPAL track trigger 



CHAPTER 3. THE OPAL DETECTOR 



21 



and in pattern recognition. A more precise z measurement is then obtained 
offline by combining axial and stereo drift time information offline. Multiple 
hits on a wire can be recorded. 

3.1.4 Silicon Microvertex Detector 

The Sihcon Microvertex Detector (SI) consists of two barrels of single sided 
Silicon Microstrip Detectors at radii of 6 and 7.5 cm. The inner layer consists 
of 11 ladders and the outer of 14. Each ladder is 18 cm long and consists of 
3 sihcon wafers daisy chained together. There are 629 strips per detector at 
25 /im pitch and every other strip is read out at 50 /im pitch. The detector 
was originally installed in OPAL in 1991 and had r — cf) readout only. In 1993 
an upgraded detector was installed that has r — (f) and r — z wafers glued 
back to back. 

3.2 The Electromagnetic Calorimeter 

The function of the electromagnetic calorimeter (ECAL) is to detect and 
identify electrons and photons. It consists of a lead glass total absorption 
calorimeter split into a barrel and two end cap arrays. This arrangement to- 
gether with two forward lead scintillator calorimeters of the forward detector 
makes the OPAL acceptance for electron and photon detection almost 99% 
of the solid angle. 

The presence of ~2 radiation lengths of material in front of the calorimeter 
(mostly due to the solenoid and pressure vessel), results most of the elec- 
tromagnetic showers to start before reaching the lead glass. Presampling 
devices are therefore installed in front of the lead glass in the barrel and end- 
cap regions to measure the position and energy of showers to improve overall 
spatial and energy resolution and give additional •y/n^ and electron/hadron 
discrimination. In front of the Barrel Presampler is the Time Of Flight de- 
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tector. 

3.2.1 Barrel Lead Glass Calorimeter 

The barrel lead glass calorimeter (EB) consists of a cylindrical array of 9440 
lead glass blocks at a radius of 246 cm covering the polar angle range | cos 9\< 
0.81. Each block is 24.6 radiation lengths, 37 cm in depth and ~ 10 x 10 cm^. 
In order to maximize detection efficiency the longitudinal axis of each block is 
angled to point at the interaction region. The focus of this pointing geometry 
is slightly offset from the e''"e~ coUision point in order to reduce particle losses 
in the gaps between blocks. 

Cerenkov light from the passage of relativistic charged particles through the 
lead glass is detected by 3 inch diameter phototubes at the base of each 
block. 

3.2.2 Endcap Electromagnetic Calorimeter 

The endcap electromagnetic calorimeter (EE) consists of two dome-shaped 
arrays of 1132 lead glass blocks located in the region between the pressure 
bell and the pole tip hadron calorimeter. It has an acceptance coverage of 
the full azimuthal angle and 0.81 <| cos^^|< 0.98. 

As opposed to the barrel calorimeter, the endcap lead glass blocks follow a 
non-pointing geometry being mounted coaxial with the beam line. The lead 
glass blocks provide typically 22 radiation lengths of material and come in 
three lengths (38, 42 and 52 cm) to form the domed structure following the 
external contours of the pressure bell. 

The blocks are read out by special Vacuum Photo Triodes (VPTs) operating 
in the full OPAL magnetic field. 
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3.2.3 Barrel Electromagnetic Presampler 

The Barrel Electromagnetic Presampler (PB) consists of 16 chambers forming 
a cylinder of radius 239 cm and length 662 cm covering the polar angle range 
I cos^|< 0.81. Each chamber consists of two layers of drift tubes operated in 
the limited streamer mode with the anode wires running parallel to the beam 
direction. Each layer of tubes contains 1 cm wide cathode strips on both sides 
at ±45° to the wire direction. Spatial positions can then be determined 
by reading out the strips in conjunction with a measurement of the charge 
collected at each end of the wires to give a z coordinate by charge division. 
The hit multiplicity is approximately proportional to the energy deposited 
in the material in front of the presampler allowing the calorimeter shower 
energy to be corrected with a corresponding improvement in resolution. 

3.2.4 Endcap Electromagnetic Presampler 

The endcap electromagnetic presampler (PE) is an umbrella shaped arrange- 
ment of 32 chambers in 16 wedges (sectors). It is located between the pressure 
bell of the central tracking system and the endcap lead glass calorimeter, cov- 
ering the full azimuthal angle in the polar angle range 0.83 <| cos6'|< 0.95. 

3.2.5 Time- Of-F light Counters 

The time-of-fiight (TOF) system provides charged particle identification in 
the range 0.6 to 2.5 GeV, fast triggering information and an effective rejection 
of cosmic rays. 

The TOF system consists of 160 scintillation counters forming a barrel layer 
684 cm long at a mean radius of 236 cm surrounding the OPAL coil covering 
the polar angle range | cos6'|< 0.82. 
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3.3 The Hadron Calorimeter 

The hadron calorimeter (HCAL) is built in three sections - the barrel, the 
endcaps and the pole-tips. By positioning detectors between the layers of 
the magnet return yoke a sampling calorimeter is formed covering a solid 
angle of 0.97 x An and offering at least 4 interaction lengths of iron absorber 
to particles emerging from the electromagnetic calorimeter. Essentially all 
hadrons are absorbed at this stage leaving only muons to pass on into the 
surrounding muon chambers. 

To correctly measure the hadronic energy, the hadron calorimeter information 
must be used in combination with that from the preceding electromagnetic 
calorimeter. This is necessary due to the hkelihood of hadronic interactions 
occurring in the 2.2 interaction lengths of material that exists in front of the 
iron yolk. 

3.3.1 Hadron Endcap and Barrel Calorimeter 

The barrel region (HB) contains 9 layers of chambers sandwiched between 8 
layers of 10 cm thick iron. The barrel ends are then closed off by toroidal end- 
cap regions (HE) which consist of 8 layers of chambers sandwiched between 
7 slabs of iron. 

The chambers themselves are limited streamer tube devices strung with an- 
ode wires 1 cm apart in a gas mixture of isobutane (75%) and argon (25%) 
that is continually flushed through the system. The signals from the wires 
themselves are used only for monitoring purposes. The chamber signals re- 
sult from induced charge collected on pads and strips located on the outer 
and inner surfaces of the chambers respectively. 

The layers of pads are grouped together to form towers that divide up the 
detector volume into 48 bins in and 21 bins in 6. The analogue signals from 
the 8 or so pads in each chamber are then summed to produce an estimate 
of the energy in hadronic showers. 
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The strips consist of 0.4 cm wide aluminium that run the full length of the 
chamber, centered above the anode wire positions. They hence run parallel to 
the beam line in the barrel region and in a plane perpendicular to this in the 
endcaps. Strip hits thus provide muon tracking information with positional 
accuracy limited by the 1 cm wire spacing. Typically, the hadronic shower 
initialized by a normally incident 10 GeV pion produces 25 strip hits and 
generates a charge of 600 pc. 

3.3.2 Hadron Pole-Tip Calorimeter 

The pole-tip hadron calorimeter (HP) complements the barrel and endcap 
ones by extending the sohd angle coverage 0.91 < | cos^|< 0.99. Here the 
gap between the iron plates, available for detectors, was reduced to 10 mm to 
avoid perturbing the magnetic field. In order to improve the energy resolu- 
tion in the forward direction, where the momentum resolution of the central 
detector is falling off, the distance between samplings was reduced to 8 cm 
and the number of samplings increased to 10. 

The detectors themselves are 0.7 cm thick containing a gas mixture of CO2 
(55%) and n-pentane (45%), strung with anode wires at a spacing of 0.2 cm. 
Again, the chambers have pads on one side (of typical area 500 cm^) and 
strips on the other. Corresponding pads from the 10 layers then form towers 
analogous to the treatment in the rest of the calorimeter. 

3.4 The Muon Detector 

The muon detector aims to identify muons in an imambiguous way from a 
potential hadron background. To make the background manageable, particles 
incident on the detector have traversed the equivalent of 1.3 m of iron so 
reducing the probability of a pion not interacting to be less than 0.001. 
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3.4.1 Barrel Muon Detector 

The barrel region (MB) consists of 110 drift chambers that cover the accep- 
tance I cos^l < 0.68 for four layers and | cos^| < 0.72 for one or more layers. 
The chambers range in length between 10.4 m and 6 m in order to fit be- 
tween the magnet support legs and all have the same cross sectional area of 
120 X 9 cm^. 

Each chamber is split into two adjoining cells each containing an anode signal 
wire running the full length of the cell, parallel to the beam line. The inner 
surfaces of the cells have 0.75 cm cathode strips etched in them to define the 
drift field and in the regions directly opposite the anode wires are diamond 
shaped cathode pads. In all, six signals are read out from each cell namely, 
one from each end of the anode wire and four from the cathode pads and 
these are digitized via an 8-bit FADC. 

Spatial position in the plane is derived using the drift time onto the anode 
and can be reconstructed to an accuracy of better than 0.15 cm. A rough 
estimate of the z coordinate is also achieved by using the difference in time 
and pulse height of the signals arriving at both ends of the anode wire. A 
much better measure of the z coordinate is given by using induced signals 
on two sets of cathode pads whose diamond shape repeats every 17.1 cm 
and 171 cm respectively. This results in a z coordinate accurate of 0.2 cm, 
modulo 17.1 cm or accurate to 3 cm modulo 171 cm. 

3.4.2 Endcap Muon Detector 

Each endcap muon detector (ME) consists of two layers of four quadrant 
chambers (6x6 m^) and two layers of two patch chambers (3 x2.5 m^), for an 
angular coverage of 0.670 <| cos6'|< 0.985. Each chamber is an arrangement 
of two layers of limited streamer tubes in the plane perpendicular to the 
beam line, where one layer has its wires horizontal and the other vertical. 
The basic streamer tube used has a cross section of 0.9 x 0.9 cm^ with the 
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inner walls coated with a carbon-suspension cathode. Each plane of tubes is 
open on one side and closed on the other to rows of aluminium strips 0.8 cm 
wide. The strips on the open side, run perpendicular to the tube anode wires 
and typically have charge induced over five or so strips. By finding a weighted 
average using the recorded pulse heights, the streamer is located to better 
than 0.1 cm. The strips on the closed side run parallel to the tube wires and 
so can only give that coordinate to the nearest wire or 0.9/vT2 cm. 
Within each chamber therefore, with two layers of tubes each with two layers 
of strips, the x and y coordinates of a track can be measured once accurately 
and once relatively coarsely. As with the barrel region, the actual position 
of the strips is known to about 0.1 mm via survey information. 

3.5 The Forward Detector 

The forward detector (FD) consists of an array of devices, listed below, whose 
primary objective is to detect low angle Bhabha scattering events as a way of 
determining the LEP luminosity for the normalization of measured reaction 
rates from decays. 

To achieve this, the forward detector enjoys a relatively clean acceptance for 
particles between 47 and 120 mrad from the interaction point, with the only 
obstructions being the beam pipe and 2 mm of aluminium from the central 
detector pressure vessel. 

• Calorimeter. The forward calorimeter consists of 35 sampling layers 
of lead-scintillator sandwich divided into a presampler of 4 radiation 
lengths and the main calorimeter of 20 radiation lengths. 

• Tube Chambers. There arc three layers of proportional tube chambers 
positioned between the presampler and main sections of the calorime- 
ter. The positioning is known to ±0.05 cm and they can give the 
position of a shower centroid to ±0.3 cm. 
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• Gamma Catcher. The gamma catcher is a ring of lead scintillator 
sandwich sections of 7 radiation lengths thickness. They plug the hole 
in acceptance between the inner edge of EE and the start of the forward 
calorimeter. 

• Far Forward Monitor. The far forward monitor counters are small lead- 
scintillator calorimeter modules, 20 radiation lengths thick, mounted 
either side of the beam pipe 7.85 m from the intersection region. They 
detect electrons scattered in the range 5 to 10 mrad that are deflected 
outwards by the action of LEP quadrupoles. 

3.6 The Silicon Tungsten Detector 

The sihcon tungsten detector (SW) is a samphng calorimeter designed to 
detect low angle Bhabha scattering events in order to measure the luminosity. 
There are 2 calorimeters at ±238.94 cm in z from the interaction point with 
an angular acceptance of 25 to 59 mrad. Each calorimeter consists of 19 
layers of sihcon detectors and 18 layers of tungsten. At the front of each 
calorimeter is a bare layer of silicon to detect preshowering. The next 14 
silicon layers are each behind 1 radiation length (3.8 mm) of tungsten and 
the final 4 layers are behind 2 radiation lengths (7.6 mm) of tungsten. 

Each silicon layer consists of 16 wedge shaped silicon detectors. The 
wedges cover 22.5° in with an inner radius at 6.2 cm and an outer one at 
14.2 cm. The wedges are subdivided into 64 pads (32 in r and 2 in 0) giving 
a total of 38912 channels which are read out individually. Adjacent wedges 
in a layer are offset by 800 /im in z and positioned in such a way that there is 
no gap in the active area of the silicon. Consecutive layers in the detector are 
offset in (f) by half a wedge (11.25°) so that any cracks between the tungsten 
half-rings do not line up. 
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3.7 The Trigger 

Events are only recorded by the data acquisition system if they satisfy cer- 
tain trigger conditions (and since 1992 pretrigger conditions). Prom a bunch 
crossing rate of 45 kHz, in the original 4+4 bunch mode, the trigger system 
reduces the event rate to 2-3 Hz at atypical luminosity of 0.5x10^^ cm~'^s~^. 
The data are subsequently processed by a software "filter" which uses a par- 
tial event reconstruction and some preliminary event classification to reduce 
the event rate by a further ~30%. 

Subdetector trigger signals divide into two categories, "stand-alone" signals 
such as multiplicity counts or energy sums, and lower threshold signals from 
a 6 X 24 bins in 9 and respectively. The trigger processor makes its decision 
by forming correlations in space between subdetectors in. 6 /(f) together with 
the stand-alone signals. 

The original trigger system was designed as a single stage trigger, because the 
time between a bunch crossing and the decision being made plus the time 
required by detector components to clear their electronics ("reset time") 
amounts to ~20 /xs, which is less than the interbunch time of 22 jis in 4+4 
bunch mode. 

Prom the start of data-taking in 1992, the detector has run with a two- 
stage trigger system suitable for operation with more than 4+4 bunches. 
With 8+8 bunches in LEP (and potentially even more bunches), the pretrig- 
ger performs a deadtime-free first level decision without compromising trigger 
efficiency or acceptance, reducing the 90 kHz bunch crossing rate down to a 
rate of positive pretriggers of 1-2 kHz. Similarly to the trigger system, subde- 
tector pretrigger signals divide into "stand-alone" signals from energy sums, 
and lower threshold signals from 12 bins in O/cf). The pretrigger processor 
makes its decision by multiplicity counting and possibly forming correlations 
in 9 /(f) between subdetectors, together with the stand-alone signals. 
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3.8 Online Dataflow 

When a beam crossing is selected by the trigger as containing a potentially 
interesting event, the subdetectors are read out. Each one of the sixteen 
subdetectors is read out separately by its own special front-end readout elec- 
tronics into its local system crate(s) ("LSC") |2^. The subevent structures 



from the different LSCs (eighteen of them, including the trigger and track 
trigger) are assembled by the event builder ("EVB"). 

When the complete events have been assembled by the EVB, they are passed 
in sequence to the filter program. In the filter, the events are checked, ana- 
lyzed, monitored and compressed, where some obvious background, typically 
15-35 % of all triggers, are rejected. At this stage, a MH event typically 
occupies 210 kbytes of data storage and other events 60 kbytes. This is 
reduced by the filter by an average factor of five, and than events are written 
to a buffer disk. The buffer with a capacity of several hours of data taking is 
used to decouple the data acquisition from subsequent event reconstruction 
and data recording. As a backup the events are also copied from the filter 
disk buffer to IBM cartridges tapes. 

The event reconstruction program, "ROPE" based on a network of HP 
Apollo workstations, fetches the event files from the filter buffer disk, and 
records them on optical disks for long term storage. The full reconstruc- 
tion is performed as soon as up-to-date calibration data from the LSCs are 
available, usually within an hour of the events being taken. 

The event data are in ZEBRA format data structure. Each ZEBRA 
structure of a complete event includes a "header" with 64 words of basic 
event information, such as trigger pattern and filter event classification, it 
can be updated during the event processing and is used for fast selection of 
events of a particular type. 

Monitoring is performed at various levels in the data acquisition chain. 
The LSCs perform detailed monitoring of subdetectors. At filter level, com- 
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plete event are given, so correlation between subdetectors can be made and 
event classification can be used. The results are immediately available. The 
ROPE monitors the combined detector performance after full processing has 
been performed. The events are available for analysis when moved to a per- 
manent storage area. Data are often ready for analysis within a few hours 
from the end of a run. 



Chapter 4 

Theoretical Background 



4.1 A Brief Review of the Standard Model 

A common aim of all science is to explain as many facts as possible with a 
few simple principles. This leads to the efforts to relate the known phenom- 
ena and the attempts to reach unification in our theories (see for example 
reference P^]). 

Our experimental evidences nowadays suggest that all matter is composed 
of structureless, point like, quarks and leptons, both obey Fermi statistic 
rules, and the interactions between them are mediated by gauge massive and 
massless bosons. The theories describing these interactions are all gauge 
theories which can be described as a change from global symmetry to local 
symmetry. The nature of transmission of forces is the interaction between a 
gauge field, and a conserved matter current. 

The fermions are grouped in three families (sometimes refers as 

generations), each consist of two quarks, a massive charged lepton, and its 
light (or even massless) uncharged partner, the neutrino. In rising weight 
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order the families are: 
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Quarks 
Leptons 



The forces of nature are traditionally reduced to 4 fundamental forces: 

• Electromagnetic, carried by the uncharged massless particle, the pho- 
ton, 7. 

• Weak nuclear force, mediated by the massive bosons W~^, and W~ . 

• Strong nuclear force, carried by the massless boson, the gluon, g. 

• Gravitation, carried by the massless particle with spin 2 the graviton, gr- 

QED, quantum electrodynamics, is the exact theory of electromagnetism. 
The agreement between the theory and experiment is excellent, in fact QED 
is the most precise theory in the field of physics. 

The weak nuclear force, is the force which is responsible, among other 
things, for the P decay: n ^ peu^ and the muon decay /i — > ez/^z/g- The 
weak force is known to be of a short range (~ 10~^^ m) with a maximal 
violation of parity. 

Quantum-Chromo- Dynamics, QCD, the strong field theory, is the gauge 
field theory which describes the interaction of colored quarks and gluons. The 
principle of "asymptotic freedom" determines that the renormalized QCD 
coupling to be small only at high energies, and therefore only at that region 
high precision tests can be performed using perturbation theory. 

QCD, the weak force, and QED are the three components of the Standard 



Model (SM) SU(3) x SU(2) x U(l) theory. 
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The Glashow-Weinberg-Salam Theory 

The V — A weak interaction is only a phenomenological theory, because it 
is not renormahzable. Attempts to construct renormalizable weak interaction 
theory have failed, until Glashow (1961) ||l|, Weinberg and Salam (1967) H 
constructed a spontaneously broken gauge theory, which unified the weak and 
the electromagnetic interactions. After a large number of different versions 
of this theory were tested during the 1970's, it turned out that it is possible 
to describe all the experimental data using Glashow-Weinberg-Salam theory. 
One of the shortcomings of this model is its many parameters which their 
values are not predicted by the theory. 

The electroweak interactions are mediated by 4 massless vector bosons 
Wi, W2, Ws and B. The W triplet couples to the weak isospin of the fermions, 
and the singlet B couples to the "hypercharge" [Y ) a. combination of weak 
isospin ( / ), and the electric charge ((?/): Y = 2qf — I^. 

The theory formulated in this way describes interactions between massless 
fermions, a postulate which clearly disagrees with the experimental observa- 
tions. It is also involves with massless gauge bosons and therefore long range 
forces, which again is not realistic because it is experimentally known that 
the weak interaction are short-ranged. For these reasons a scalar complex 
"Higgs" field, which spontaneously breaks the SU{2) x U{1) symmetry, was 
introduced to the theory. The spontaneous symmetry breaking leads to four 
physical particles: 

• Two massive charged particles and W~ , which are responsible for 
the charged weak currents. 

• Massive Z'^ which is responsible for weak neutral interactions. 

• The massless particle 7, which is responsible for electromagnetic inter- 
actions. 
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The photon,7, and the Z° are correlated through: 

sinOw cos 9w \ ( \ (4^1) 
cos 6w — sin 6w J \ B J ^ ' ^ 

where the mixing angle 9w (Weinberg angle) is a free parameter. The masses 
of the four real bosons are related through the relations: 

37.3 [GeV] 37.3 [GeV] ^ 

mw = ^ , ruzo = ^— = 0. (4.2) 

sm dw sm dw cos dw 

The insertion of the Higgs field generates masses to the fermions. How- 
ever, the fermions mass values are not predicted by the theory, and are left 
as free parameters. Another value which is almost unlimited is the mass of 
the Higgs boson which generates this spontaneously symmetry breaking. 

There is a slight difference between the quarks listed at the beginning of 
this section which are essentially the mass eigenstates and the eigenstates 
of the elctroweak model. The d'-type quarks (0?', s', h') which arc the eigen- 
states of EW are mixed with the mass eigenstates via the Cabbibo Kobayashi 
Maskawa (CKM) matrix: 

/ d' 

s' \ = \ Vcd Vcs Kb I I 5 I (4.3) 




The CKM matrix elements are not predicted by the SM hence they must be 
input from the experimental measurements. 

In addition to the fermion masses and the CKM elements there are still 
three more degrees of freedom left to be experimentally obtained, and there 
are several parameter combinations one can choose to fix the theory with. A 
common set of three parameters is for example: 



1. Q-w - the Weinberg mixing angle defined with tan^^K = 9' 1 9 the ratio 
between the U(l) and the SU(2) coupling constants. 
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2. e - the electric charge of the electron which before radiative correction 



= g'^ sin^ 9\, 



3. Mw or Mz - the mass of the or Z° bosons where without radiative 
correction hold the relation sin^ 9]y = 1 — M^/M|. 

The weak charged currents couple to fermions hj a.V — A interaction. The 
electromagnetic interaction has only a vector coupling. The neutral currents 
are a mixture of V and A, given by the couplings 

sm Uw COS Uw 

where /^'^ is the weak isospin of the left or right handed fermions. 1^ = 
for all right handed fermions, and for the left handed doublets: 






^' Jl 

the upper components have Is = +| and the lower components have 
h = -l 

The vector, v, and the axial-vector, a, couplings are combinations of the 
left and right couplings: 

v = l{9^ + 9% a = l{g^-g^) (4.5) 

The difference between the left and right coupling causes parity violation and 
helicity effects. 

Table |4.1| presents the SM weak coupling constants for the lepton and 
antilepton doublets. 

The electroweak model is currently consistent with all experimental find- 
ings and is generally accepted as giving the correct description of these in- 
teractions, even though some ingredients are still not well tested. The LEP 
experiments in general, and the r decays tests in particular, represent a po- 
tential field for making rigorous tests of the SM predictions, measurements 
of its free parameters, and searching for possible deviations from the model. 



CHAPTER 4. THEORETICAL BACKCROUND 37 





gH 






a 


e+,r+,/i+ 


sin^ 6w 
— sin^ 9w 




— 1 + sin^ 9w 

i — sin^ 6*1^ 
1 

2 
1 

2 


-\ + sin^ 6^14/ 

i - sin^ Ow 

1 

4 
1 

4 


1 

4 

1 

4 
1 

4 

1 

4 



Table 4.1: The weak coupling constants of the leptons doublets in units of 
1 / sin 9w cos 9w 



4.2 The r- Asymmetry Formalism 

The annihilation of e^e~ into tau pair supplies information on the relative 
strength of the neutral vector and axial- vector couplings to the charged heavy 
lepton. This tallies with the fact that one of the main physics goals of 
the LEP experiments, is a precision measurement of the electroweak mixing 
angle, siii^ 9w, in many different ways. 

The r pairs are produced with a forward backward asymmetry, a longi- 
tudinal polarization and a forward backward asymmetry of its polarization, 
due to the interference of electromagnetic, and neutral currents 



The presence of the r polarzation is a purely parity violating effect (arising 
from the neutral current), while the forward backward asymmetry can arise 
also from higher order QED effects (interference of two photons). 

The r+r^ pair, as a system of two fermions, can be produced in four 
polarization (helicity) combinations, 

(P^-,P^+) = (+,+), (+,-), (-,+), (-,-), (4.6) 

with the corresponding cross sections, 

a(++), a(+-), a(-+), a(— ). (4.7) 

This leads to the following definitions of the average and r+ polariza- 
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tions, 



(Pr-) 



(4.8) 



(^tot 




[<7(++) + a(-+)]-[a(+-) + -a(— )] 



CTtot 



where, 



(Ttot = (7(++) + (7(+-) + (7(-+) + (7( )■ 



(4.9) 



However, in e'^e" annihilations, the t'^t" pair are known to be produced 
through an intermediate state of a spin-one boson which can be a photon or 
a Z° . Looking first at QED scattering of an electron, it is known that the 
helicity of the electron is conserved at high energies, whereas at low energies 
the direction of the spin with respect to a fixed coordinate system in space 
is preserved. The analogous behavior is present in e'^e~ — > t+t~, namely: 

• In the ^ — > region: the helicities of the t+ and the t~ tend to be 
opposite to each other (fj,). 

• While in the ^ 1 region : the helicities of the r+ and the t~ 
tend to be parallel (H). 

In this case, the helicity conservation rule at high energies restricts the 
number of helicity combinations to two, namely (+, — ) and (— , +) or 
(P-7— = — P-7-+), whereas the cross sections for the other two combinations, 

cr{++) and a{ ) are neghgibly small. The r' and r"'" helicity relation 

stands also for their average polarization values. For these reasons we will 
denote in the following : 



(4.10) 



(Pr) ^{Pr-) 



--(Pr^)- 
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The general expression for a differential cross section of the reaction 
e~^e~ T^T~ intermediates by a vector boson is given by, 

da 



-(s, cos ^,Pr 

dcosa 



(1 + cos2^)Fo(s) + 2cos0Fi(s) 
Pr[(l + cos^ 0)F2{s) + 2cos^F3(s)]. 



(4.11) 



Where the Fi{s) are form factors. 



Fig. [4.1| presents the two SM e^e^ annihilation diagrams for the produc- 
tion of a r pair: one via 7 and the other through a Z° exchange (while the 
Higgs exchange channel can be neglected p3| ) 





) r\-/r\-/r\-/ >■ 


T + 

{ 


+ 


e+ 

\ 




T + 

( 


e~ 




T 








T 



Figure 4.1: Feynman diagrams for the r pair production in e^e annihilation 
through the single photon and exchange. 

The tree level equation for the process e~^e~ — ^ r+r" is a linear combi- 



nation of the three terms ^ 

• Pure QED contribution : 
cr^ = 1 + cos^ 9. 

• Electro-weak interference contribution: 

= —2qf[veVr{l + cos^ 9) + 2aear cos9]Rex- 

• Pure weak interaction contribution: 

azo = [(f f + ) (f ^ + ) ( 1 + cos^ 9) + SveVrttear COS 9]\x'^\. 
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where the Z° propagator, 

^■<^> = .-M| + ,.r,/M. - ^'■'''> 

describes the Z° resonant shape without radiative corrections and is sym- 
metric arounnd ^/s = Mz- With radiative corrections a has a smaller value 
at the peak and an unsymmetric shape. Mz and Tz are the Z° mass and 
width. 

The qf,Vf and aj {f = e or r) are the electric charge, the vector and 
axial-vector coupling constants of the fermions to the Z° respectively. The 
strength of the vertices following the SM structure are: 

• 6 7 6 : V/ty g = -iqel^ 

• T T : y-TjT = -iqrY 

• eZ" e : V^^^^ ^ = ' «e7^) 

• ■ Kz^ r = -'1^^^'^ - «-^') 
•rWr: r = -^7^(1 - 7^) 

The corresponding four possible differential cross-sections are: 

|geg. + T-+T^-+x(5)|2 

\qeqr + T-+T^-xis)\' 

\qeqr + T+-T-^x{s)\' (4.13) 
|geg. + T/-T+-x(s)|2 



da 


(eZ4 - 




27ra2 


(l + COS^)2 




dcos6 


s 


4 


da 


[eZ4 - 


- r-RTt) 




(1 -COS^)2 


dcosO 


s 


4 


da 


[e^e^ - 


- ^l4) 




(1 -COS^)2 


dcosO 


s 


4 


da 


(eZ4 - 


- rRrt) 




(l + COS^)2 




dcos9 


s 


4 



Here the geg^ are QED contribution while the coupling of the weak neutral 
current are 
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Tf^ = Vf±l3af. (4.14) 

where (3 is the velocity of / in speed of hght units. 

Therefore according to the SM the form-factors Fi{s) of Eq. [4.11| are given 

by, 



2 

^o(s) = "^[qIqI + '2Rex{s)qeqrVeVr + \x{s)\\vl + al){vl + a^)] 
Fi{s) = -^[2Rex{s)qeqraear + \x{s)\^'^Veae2vrar] 
F2{s) = -j^[2Rex{s)qeqrVear + \x{s)\^{v'^ + al)2Vrar] 

2 

F3{s) = ^[2Rex{s)qeqraeVr + \x{s)\^2veae{vl + al)]. (4.15) 



From Eq. [4.11| , four independent cross sections can be constructed: 



a^{s, P^ = +l) = ^[Fo{s) - F^is)] + [Fi(s) - F,{s)], 

a^is, Pr = -l) = \[Fo{s) + F2{s)] + [F,{s) + F,{s)l (4.16) 

a^(., P. = +l) = ^[Fo(.) - F2(.)] - [F,{s) - F,{s)l 

a^{s,Pr = -l) = \[F,{s) + F2{s)] - [F^{s) + F^is)], 

where the indices F (B) refer to events with scattering into the forward 
(backward) hemisphere, namely, < cos 6 < 1 (— 1 < cos 6 < 0). 

These four independent cross sections can be combined into four other 
cross sections which are simpler to measure, but are no longer independent, 

a^(s) = a^{s,Pr = +l) + a^{s,Pr = -l) = ^Fo{s) + 2F,{s), 
a^(s) = a^(s,P. = +l) +a^(s,P. = -l) = ^Fo{s) - 2F,{s), (4.17) 
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a{s,Pr = +l)=a^{s,P, = +l)+a''{s,Pr = + l) = -[Fo{s) - F^{s)], 
a{s, Pr = -1) = a^{s, Pr = -1) + P. = -1) = ^ [^o(s) + F^i^s)], 
and the total cross section is given by, 

16 

a{s) = a^{s) + a'^is) = a{s,P^ = + l) + a{s,P^ = -l) = —Fo{s). (4.18) 



Summing over both hehcity states of Eq. |4.11| one obtains the following 
cos^ distribution, 

— -^ = 1(1 + cos2 e + fA^BCOs^^), (4.19) 

(Ttot uCOS o 

utilizing the forward-backward asymmetry, ApB-, defined as, 

A,,,.) . ^!(£i_^ . (4.20) 

a{s) 4Fo(s) 

In the same way the polarization asymmetry, (Pr), (= {Pt){s)) can be 
written as, 

a[s) Fq[s) 
consequently, the cross-section can be divided into 

a+ = \(l + {P,))at,t 



i(l - i,Pr))atot. (4.22) 



2 



The polarization asymmetry can be also defined separately for each hemi- 
sphere, 

(P Y ^ o^M^+l) - a^(s,P. = -l) 



(4.23) 

4P2(g)+3P3(g) 

4Po(s)+3Pi(s)' 
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whereas for the backward hemisphere 



(4.24) 

4Fo(s)-3F.(s) 



4Fo(s)-3Fi(s) 



The forward-backward polarization asymmetry is then defined as 

A' 



_ [a^{s,Pr = + l) - a^is,P, = -l)] - [cr^{s,P, = +l) - = -1)] 

^P"' ~ a(s) 



(4.25) 



4Fo(s)' 



Using these definitions one obtains, 

:[(1 + (P.))(l + cos2 6) + UAfb + A£;f )cos^^] 



1 da^ 



^ ^[(l-(P.))(l + cos2^) + |(A^^-A^,f)cos0].(4.26) 



atot dcos 9 



From the expressions above, the average polarization for a given polar angle 
6 is given by, 

(P.)(l + cos^g) + |Ag,fcosg 

(1+C0S2^) + |A^BC0S^ ■ ^ ^ 



which is illustrated in Fig. 

Utilizing the three asymmetries ((Pr), A^r^ and A^^f ) as were determind 
here, the expression in Eq. |4.11 of the differential cross section for an explicit 



r helicity state (Pr = ±1) can also be writen as, 
da 



s,cose,Pr) = Fo{s){{l + cos^ 6) + -Ap b{s) cose (4.28) 



d cos 9 3 

pol 



+Pr[{Pr)is)il + cos' 9) + -A^,f (s) COS^]}. 
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TAU POLARIZATION (COS0) 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



cos9 

Figure 4.2: The average tau polarization as a function of the polar angle. 
The B and F stand for the forward and backward regions. 

The distributions described by Eqs. [4.26| cannot be directly measured, 
because it is not possible to determine the r helicity on an event-by-event 
basis. Instead, since the and the r~ decay via weak interaction, where 
parity conservation is maximally violated, the angular distribution of the de- 
cay products depends strongly upon the spin orientation of the r. As the r~ 
is expected to be produced with a negative polarization and the r"*" with a 
positive one, we expect to be able to measure in the lab. system a deviation 
from an unpolarized distribution, while measuring the momentum spectrum 
of the T decay products. Thus for the determination of the average polar- 
ization we are using the kinematical distributions of the r decay products, 
depending on the helicity and the decay mode. 

The drawings in Fig. illustrate the r~ and r+ decay configurations 
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including their spin orientation. In these we are looking at an arbitrary decay 
of T to lepton, T^iuiUT {i = e, /i). Letting mi 0, leads to negative £, ui, z/^ 
helicities, while the get positive helicities. 






(b) 



Figure 4.3: and decay configurations including their spin orientation, 
a) Is an arbitrary configuration of r~ decays to i. b) Is charge conjugation 
of (a) which is not a physically state because the ui, zT^, have forbidden 
helicities. c) Is an allowed state, the mirror image of (b). 

If CP invariance is valid in these processes then the rate of (a) is equal 
to that of (c). Therefore, the distribution of the decay r+ can be 

obtained by changing the sign of the polarization vector in the distribution 
for i. Following this argumentation one can also relate — > tt^ and 

the — > 7T~ in the same way. 

For T^iiJii'T-, the relevant kinematical variable for measring the r polar- 
ization is the lepton energy, scaled by the beam energy, x = Ei/Ebeam- The 
general form of energy and angular distribution of the lepton decaying from 
an arbitrarily polarized tau, can be written in the r rest frame as ||31|| : 

5 1 

T{T^iueu,) = 1^ I dne dx x^[i -2x- {Pr ■ p}){2x - 1)] (4.29) 

where Pr is the polarization vector of the r, p£ is the unit vector along the 
direction of the lepton. The polarization dependent term Pr ■ pe is due to 
the parity violation. 

Near x=l , the lepton is at its highest energy value, and the relative mag- 
nitude of the parity violation term is maximal. Here the i tends to be emitted 
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opposite to the direction of the r~ spin, whereas the i'^ tends to be emitted 
parallel to the r"*" spin direction. 

Near x=0 . exactly the opposite holds. 

To illustrate this behavior let us look at the following diagrams. 

At X = 1 , kinematics require that the z/, and the V are both emitted in the 



opposite direction of the lepton (see Fig. [4.4|) . Since the component of the 
orbital angular momentum is zero along the i direction, and the two neutrino 
spins add up to zero, the spin of the lepton must be parallel to that of the 
T. Since the i has a negative helicity and the has a positive helicity, the 
i tends to be emitted opposite to the r~ spin. 



X=l 



1+ 



X=l 



X=l 



X=l 



favored 



forbidden 



favored 



forbidden 



Figure 4.4: T—^iv^Vj. decay in the r rest frame for x = 1. 



When a: ^ (Fig. [4.5|) , the kinematics require that the p and V come 
out in opposite directions to each other, hence their net spin is equal to 
unity and points toward the direction of the V. In order to conserve angular 
momentum, the ^ has to move in the direction of the V and the spin of the 
T~ points in the direction of V. Hence, near x = 0, the ^ tends to come out 
along the direction of spin of the r~, which is exactly the opposite to the 

X —— X CcLSG. 

Boosting the angular distribution of the T^ivj/Vr from the r into the lab. 
system one obtains the x distribution as given by [pl[| . 
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T 



T 



x ~ 



x ~ 



x ~ 



a; ~ 



favored 



forbidden 



favored 



forbidden 



Figure 4.5: T^iviv^ decay in the r rest frame for very low x. 



' i(5-9x2 + 4x^) +Pri(l -Qx^ + Sx^^). (0 < a; < 1) (4.30) 



Ti dx 

This distribution, as well as the following expressions for the other decay 
channels, holds for as well as for r+ provided that Pj- is always taken as the 
r~ helicity. Terms of order mi/mr have been neglected. This approximation 
is fully justified for i = e, whereas for £ = /i there is a threshold effect around 
X = 0.005 which has been accounted for in our analysis. 

The angular distribution of nlK)^ from a polarized is easier to un- 
derstand. As it is a two-body decay the energy of each particle in the r rest 

^2 , 2 ^2 2 

frame is nxed: = — 2m ^ ^"^^ ~ — 2m ~ ^5-^)- ^mce the 

helicity of the vt is negative it prefers to be emitted opposite to the r spin 
direction. Hence the ti{K)~ prefers to be emitted in the r~ spin direction 
(Fig. 

The same kinematical variable used in the leptonic case can be used also 
for T^Tx{K)vj. decay. Here it is related to cos 6'*, where 9* is the decay angle 
of the hadron in the r rest frame, 

cos<^-^ ^f,-^-;f"' . (4.31) 
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TT 



TT 



TT 
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forbidden 



favored 
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Figure 4.6: r^7r(K)z/T- decay and the r helicity direction. 



Henceforward, the r-velocity term, /?, in the denominator will be approxi- 
mated to 1 (1 — /9 = 7.6 ■ 10~^ for Ecm = itt-z)- The cos6'* distribution is 
given by pi. 



1 dr 



h 



i(l + Prcosr) 



(-1 < cosr < 1), 



Th dcos 9* 

corresponding to the following x = Ph/Ebeam distributions. 



(4.32) 



1 dr, 



1 



Th dx 1 — m\/rn^ 



2x — 1 — m|/m^' 
I -ml/ml ^ 



1 + Pr "'". ^ o'T^V"" 1 {ml/ml < x < 1) 

(4.33) 

Fig. |4.7| presents the x distribution of r-^iugUj. and r^7r(K)z/^ events for 
positive and negative r helicity states, as parametrised in Eq. |4.3CI| for the 
leptonic decays and in Eq. |4.33| for the tt{K) t decays. In this analysis we 
measure the hadron momentum, ph, and do not distinguish between tt and 
K. Therefore, cos 6* cannot be calculated event-by-event using Eq. [4.31j and 
hence we choose the variable x taking into account the tt/K mixture by using 
the distribution, 



]_dl\ 
Th dx 



1 dT^ 1 dTx 

1 - bK)^; — 1 ^ Ok; 



dec 



K 



dx 



(4.34) 
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Figure 4.7: The x distribution of a) r^iuiUr events, b) r^7r(K)i/i- 
events. The sohd (dashed) hne is for positive (negative) r hehcity. 

where, bx = BR{t-^Ku) /[BR{t^^u) + BR{t^Kv)] = 0.064 ± 0.007 [||]. 
The uncertainty in bk is included in our systematic study (see Section |6^ ). 
The difference between the momentum and energy of the hadron is also 
accounted for in our analysis. 

In all r decay channels considered above, the kinematical distribution is 
linear in P^-, and can be written in the general form, 

^^ = f{x) + Prg{x). (4.35) 

Here, and in the remainder of this section, x is a generic name of the relevant 
kinematical variable. As shown above, f{x) and g{x) may differ from one 
decay channel to another, but they always satisfy the following normalization 
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conditions, 

J f{x)dx = 1 J g{x)dx = 0. (4.36) 

Since in Z° t^t^ the hehcities of the two tau leptons are expected 
(assuming CP invariance) to be fully anti-correlated |^T|, their measurement 
carry less information compared with the case of r's from two different events. 
Until now previous analyses have neglected this effect In order to take 

these anti-correlations into account, one has to analyze the double identified 
T events using the triple differential cross section with respect to cos^, Xi 
and Xj of the two r-decay products. 

The joint distributions of the r-pair production and decay is obtained 



from Eqs. [4.26| by combining them with the corresponding decay distribu- 
tions of the r~ and the r"*" and summing up the two helicity configurations, 
resulting in, 

d^a- 

d cos edx^dx, = ^^^^^ ^^^'=±±^ ± + ^) + ± K<^yos9] 

xif,±g,)if,±g,)}. (4.37) 

Here, aij is the cross section to produce r's decaying into channels (In 
the following we shall drop out the arguments Xj and Xj of the functions 
(gi) and fj {gj). This expression accounts for the correlation between the 
decay distributions of the two r's, as it should be utilized when analyzing 
events where both r decay channels are identified. 

When one of the r-decay channels (e.g. the r"*") is not identified, we have 
to integrate over its kinematical variable [xj) and we are left with, 

d'^aj _ 
dcos 9 dxj ^ 



hel=± 

±[(P,)(1 + cos' 9) + lA^^fcose]g,}. (4.38) 

When both r-decay channels cannot be identified, we are left with Eq. [4.19| 
and those events contribute information only to the forward backward asym- 
metry, ApB- 
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Defining Ae and At as, 

_ 2veae . _ 2vT-ar , , 

+ Cte "^T + 

then on the Z^-peak, the SM in the improved Born approximation, neglecting 
the contribution of the intermediate photon yields the following relations for 

.FB 
^pol ' 



(P.), A^sand Af^ 



AFBis = ml) = -^-^= AeAr 

Foimz) 4 



{PrKs = ml) = -^^ = -At (4.40) 
\FB( _ 2\ _ 3F3(m|) _ 3 



Assuming the SM with lepton universality, one can write within the frame- 
work of the improved Born approximation, 

2(1 -4sin2^VF) 
1 + (1 — 4sm Ow> 

or the ratio between the vector and the axial- vector coupling, vi/ai, can than 
be given as 

v^/ai = 1-4 sin^ 9^ (4.42) 

Fig. shows the relations between the three asymmetries, (P^), A^^^ 
and Apgf , and the mixing angle, siia^ 9w- Fig. |4.9| presents the asymmetries 
dependence on the center of mass energy. 

4.3 The Charged Weak Decay Structure 



In the framework of the SM the decay of the r proceeds via a charged 
— A^ type current. A current of this kind has been choosen to reproduce 
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Tau Asymmetries in SM 
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sme„ 



Figure 4.8: The r asymmetries as a function of the Weinberg mixing angle. 
The (-Pr), ^FB and A^^f are presented by a sohd hne, dotted and dashed 
hnes respectively. 
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Tau Asymmetries in tlie SM 
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Figure 4.9: The r asymmetries as a function of -Ecm. The {Pt)i ^fb and 
Apgf are presented by a solid line, dotted and dashed lines respectively. 
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the observed maximal parity violation in the charged weak interactions such 
as the (3 decays. The V — A characteristic of the current was well tested in 
decays. As the r assymetries are extracted from the kinematical distribution 
of the T decay products, in the discussion of the polarization measurements 
it was implicity assumed that the r decay has also the same pure V — A 
decay structure. However, any deviation from this behavior can directly 
effect our measurement. Thus, in the following we shall present again the 
Born term of the r decay differential cross section in a more general form, 
without restricting ourselves to the V — A assumption. The matrix element 
of a r leptonic decay can be written as a four fermion interaction in the 
following way 

AC 

M = ^ y: aiiim^ni^ri^.K), (4.43) 

where the index 7 labels the different types of interactions: scalar, vector or 
tensor, e and fi stand for the chiral projection of the i and r lepton spinors , n 
and m represent the handness of the U£ and the u-r- There are 19 independent 
combinations which can be measured by experiment. The SM corresponds 
to qYl = 1 (vector interaction of left handed i and r leptons) and all the 
other g]^ being zero. 

In this parametrization the lepton momentum spectrum of the r decay 
particles is modified to 

1 dT, 



Tf dx 



f + i-Pr-9 



where 



/ = /l + P-/2 

g = gi + 6-g2 (4.44) 



and 



/i = 2-6x^ + 40;^ 
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/2 

91 

92 



9(-l + 9x^ 

2 ^ 
-- + 4. 



„ 2 O ■? 

bx H — X 
3 



9 



;i - I2x + 27x^ - 16x^)] 



The p, r], 5 and ^ are the famihar Michel parameters , and are function of 
the couphng constants as given in Ref. Terms proportional to //(m^/mr) 
have been neglected as only a vector interaction will be considered here {rj = 
0). 

Table [4.2| gives the predicted values of the Michel parameters for the 
simple cases where only vector and axial-vector components are allowed at 
each of the tut and the iue vertices. 
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V -A 
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3 
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V -A 


A 


3 
8 


2 


3 
16 





V -A 


V + A 





3 









Table 4.2: The Michel parameters for some different V and A component 
combinations at the r and it lepton decay vertices. 



One can see that there are combinations other than (V — A){V — A) 
which lead to a value of 3/4 for the Michel parameter p. For this reason 
one needs also to measure the sign of the Michel polarization parameter ^ in 
order to verify the {V — A)(y — A) structure. Note that when extracting the 
assymetries looking separately at each r , one really measures the product 
^ ■ (Pr) rather than (Pr)- 
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The Matrix element of an hadronic r decay is given by 

M cx ur^.ig^ + g'^^')T- ■ (4.45) 

where stands for the hadronic current and g^^, g^ are the vector and axial- 
vector coupling constants. The parameter ^had is the chirality parameter 
defined as 2 ^ w 

= ~ WF+tfF' ^^-^^^ 

and —^had can be interpeted as twice the Ur helicity. 

Using this definition of the charged weak interaction, the x distributions 
of the r— >7r(K)z/T- decay (Eq. [4.33|) is modified to 

' ^+a„,.P,H^fiz^^ (4.47) 



h ax 1 — niyjm^ \ 1 — 

For the most general case one needs to have the following modifications 
in Eq. [O? 

• Incorporate and where both can stand for either ^ or ^^ad depend- 
ing on the decay channel. 

• Use the modified / and g terms in a leptonic decay channel including 
explicity the p and 6 parameters as in Eq. [4.44| . 



Rearranging expression [4.37] , one gets the following form, 

jj^^ - ^.,^Ej[(l + cos^^) + |A..cos.][l + (^eO-W^ 

-[(P.)(l + cos^^) + |A^^icos^][(±e.) + {±m 
x{fi{x,p)±gi{x,6)) ■ {fj{x,p)±gj{x,6)) (4.48) 

giving the possibility to extract the three asymmetries A^^^, {Pt) and A^^f , 
as well as the ^, C,had, P and 6 parameters by fitting it to the data. 
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The p parameter was aheady measured precise enough to exlude inter- 
actions different from (V" — A) at least in the iue vertex. Therefore, one can 
take a. {V — A) at the iug vertex and allow a combination of {V — A) and 
{V + A) at the r one. This assumption reduces the number of parameters 



for the charged weak interaction and implies 



3, , X 
P = g(l - hur) 

16(1 + ^) 
^ = 2 + K^ 

We see that the only way to extract the three r asymmetries and the 
charged weak parameters all together, is to measure the triple differential 
distribution of the variables X'ij Xjj COS 9 using events where both r decays are 
identified. However, one usually takes ^ = ^had = l,P = <^ = 3/4, then the r 
asymmetries can be obtained also from the double differential distributions 
of the variables Xj (or xj) and cos 6', using events where only one of the r 
decays was identified. 

It is evident from Eq. [4.48| that C, can be obtained also at low energies 



where the r asymmetries vanish. This was done by the ARGUS collaboration 
for events where one r decays to electron and the other r decays into muon, 
with the preliminary result ^ = i/^e^/i=0.90±0.13(s?/s.) it 0.13{stat.). 



In the following unless it is explicity stated we shell assume (V — A)(y — A) 
couplings, and hence the ^ and C,had as well as the other Michel parameters 
(p, 6 and rj) we will substitute by their SM values. 

4.4 Correlations Between A^^^, (P^) and A^^ 

At first sight, the measurement of the r-polarization carried via the momen- 
tum distribution of the r-decay products is independent of the measurement 
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of the r-forward-backward asymmetry based on the r- angular distribution. 
This statement is however wrong, since the angular and momentum distri- 
butions are not independent. It is a direct consequence of the r-polarization 
dependence on the r-scattering angle, which can be seen for example by com- 
paring the expressions for the polarization asymmetries in the forward and 
backward hemispheres (Eqs. [4.23| and |4.24| ). 

In order to obtain the correlations between Aps, {Pt) and A^^f , they must 
be expressed in terms of other variables which are mutually independent. We 
choose (Pt)^, (Pt)^ and the fraction ( of the forward scattering events. (Pr)^ 
and (Pt)^ are clearly mutually independent since they relate to different 
event samples. They also do not depend on the relative size of each sample 
which is determined by (. The expressions of ApB, (Pt) and Ap^f in terms 
of (Pt)^, (Pt)^ and ( are straightforward 

Afb = 2C-1, 

(Pr) = c{py+ii-c){py, (4.50) 

Using the standard error propagation procedure, the non-diagonal ele- 
ments of the covariance matrix can be calculated in terms of the variables (, 
(Pt), Apof and their statistical errors, 

(AApsAiPr)) = 2((P.)^-(P.)^)AC^ 

(AAp^AAp^^f) = 2mf+{P^f)A(', (4.51) 

(A(p.)AAp^,f) = [mfr-mmAC 

+C'(A(P.)^)2-(1-C)^(A(P.)^)2. 

Dividing those covariances by the corresponding statistical error products 
yields the correlation coefficients. 



PAMPr) = (AApbA(P.))/(AA^sA(P.)), 
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-^A.^A^^ = (AApsAA;„f)/(AA^BAA;„f), (4.52) 

and in order to obtain their numerical values, we can take the SM values for 
the asymmetries, using the sin^ 6w value of 0.2337 yielding, 

{py = 0.225 (Pr)^ = 0.032 C = U^FB + 1) = 0.506. 

These numbers correspond to the full solid angle and they must be modified 
if a cut on | cos6'| is performed on the data. 

It is straightforward to calculate the statistical uncertainty in ( using the 
binomial distribution expression, 

AC^ = (4.53) 
where N is the total number of events. Hence, 



AApB = 2AC = 2^ (4.54) 

On the other hand, the statistical errors of the polarization asymmetries are 
not so easy to calculate and they depend on the way that these asymme- 
tries are measured. The calculations will be performed in the next section, 
however a rough estimate can be done already with former values obtained 
in the experimental analyses. For example we made before our analysis (see 
Ref. 1^) an estimate of the correlation using the OPAL 1990+1991 values [0 
of A(P^)^, A{Pr)^, A{P^) and AA^^f obtained within the angular region of 



I cos^^l < .68 The correlation results are given in Table show that all the 
correlation coefficients turn out to be at the level of few percents. 

4.5 The Fit Method 



In order to choose and decide about the method of extraction the r asym- 
metris and to show the strength of our selected fit tecnique we present in the 
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r-decay mode 


number of r's 


PkFB{Pr) 
(%) 


(%) 


P{Pr)Kol 
(%) 




1809 


1.77 


2.76 


1.82 




1729 


1.64 


2.39 


1.82 


r^7r(K)z/^ 


553 


4.00 


6.10 


1.91 



Table 4.3: The SM correlation coefficients between r- asymmetries, calcu- 
lated using the OPAL 1990-1991 results and their uncertainties. 



following results of a comparison between several approaches (for details see 
Ref. BTI). 



4.5.1 Events with only one Identified r-decays 



In order to obtain {Pr) and A^^f from the experimental data, the total 
(normalized) momentum, x, of the r-decay products (excluding the i/'s) has 
to be measured. However, as shown in Eq. 4.38| the r-asymmetries A^^^, 



{Pr), A^poi ^"^^ their full covariance matrix can be obtained by fitting that 



expression to the data. One can also reduce Eq. |4.38| to a one parameter 



expression, namely sin Ow, which is related by Eq. 4.40 and Eq. 4.41 to 



{Pr) andA^i. 

For an estimate of the theoretical expected errors prior to our measure- 



ments, we have normalized both sides of Eq. [4.38| to unity, and wrote it in 
the shorthanded form, 

1 

■(s, cos6', x) = Aijfi{x)hj {cos 9), 



(fa 



a d cos 9dx 



(4.55) 



i,j=0 



where, 



and. 



A, 



00 



A 



01 



A 



FB 



/io(cos6') 



^(1 

8^ 



COS 



Aio = {Pr) An 
■ 6) hi(cos9) = cos 9 



\FB 

-^pol ! 
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In order to estimate the covariance matrix, we divide the interval of x 
into n bins and define fj' as the integral of fi{x) over the bin k {k = 1 . . .n). 
In the same way, the interval of cos 6 is divided into m bins and is defined 
as the integral of hj{cos6) over the bin / (/ = l...m). The quantity to 
minimize is then, 

n m 1 

W = Y: Etl/" - E A.ffh'^ViAy'r, (4.56) 

fc=l 1=1 i,j=0 

where y^'' is the fraction of events in the two-dimensional bin (k,l) and its 
statistical error, Ay'^', can be approximated for the case of large number of 
bins (namely, y'^' ^ 1) by. 



Ay'' = ^y'^'/N, (4.57) 

with N being the total number of r's. The inverse of the covariance matrix 
is a 3 X 3 matrix given by, 

1 dHv ^^ mm^ (4 58) 

where i = j = and i' = j' = are excluded, since Aqo = 1. Apart from the 
constant A^, the matrix V~'^ can be calculated even without data, by using 
the SM estimates for y'^K 



The results listed in Table [4.4| are obtained from the elements of the 
covariance matrix V, using 100 x 100 bins. The Table shows that the errors 
of ApB are smaller than those of (Pr) and A^^f . By comparing between the 
two decay modes, one can see that the errors on {P^-) and A^^f are smaller 
for the r— >-7r(K)z/T- decay mode. The factor between the two errors is 2.7. As 
expected, the errors on Aps do not depend on the r-decay mode. 

All the errors listed in Table |4.4| are for the ideal case. They were calcu- 
lated with the implicit assumption of no background and a uniform efficiency 
over all x- and cos^- values (0 < a; < 1, | cos^| < 1). In practice however. 
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decay mode 










0.935 
0.934 


4.58 
1.72 


4.28 
1.61 



Table 4.4: The statistical errors of the r-asymmetries obtained from a two- 
dimensional fit using the SM estimates. 



this is not the case and therefore, our calculated errors should be considered 
as lower limits. 



Table ^3] lists the correlation coefficients obtained from the non-diagonal 
elements of the covariance matrix. They are slightly higher than the values 
obtained in the last section using a simpler method and the experimentally 
measured statistical errors. This is mainly because those correlation coeffi- 
cients are inversely proportional to the statistical errors, and the values used 
here are for the ideal case. 



decay mode 


PAFBiPr) 
(%) 


Pa \fb 

^FB\oI 

(%) 


(%) 




2.24 


2.77 


4.38 


r— ^7r(K)z/^ 


5.96 


7.44 


0.18 



Table 4.5: The correlations in the r-asymmetries obtained from a two- 
dimensional fit, using the SM estimates. 



An alternative method to the two-dimensional fit is the use of two inde- 
pendent one-dimensional fits for the forward and backward scattering events. 
This method is based on the integration of Eq. [4.3^ over cos 6, separately in 
the forward (0 < cos 6 < 1) and in the backward (— 1 < cos 6 < 0) regions 
to obtain (Pr)^ and (Pr)^ ■ From (Pr)^ and (Pt)^ one calculates then (Pj.) 
and Apof • This approach has been adopted in the previous published OPAL 
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analysis [0. 

A comparison between the two methods done in Ref. |^] has used 1241 
T—>-fii?^i'T decays from the OPAL 1990 data |^. Performing a two-dimensional 



Maximum Likelihood (ML) fit, the result reached was, 

sin^ = 0.2306 ± 0.0085. (4.59) 

On the other hand, separate one- dimensional fits on forward and backward 
scattering events gave ||4^ , 



sin2^|^= 0.2342+°:r- (4.60) 

The statistical error of the last result is higher than the uncertainty of the 
value from the two-dimensional fit (Eq. [4.59|) . As one could expect, some 



information in the data was lost by changing from a two-dimensional fit to 
two separate one-dimensional fits. 

For a study of the dependence of the error on the grid size, we repeated 
the two-dimensional fit by varying the number of bins. The results are listed 
in Table [4.6|. The sin^ error values were obtained by combining the errors 



from ApB, {Pt), and A^^f , taking into account their correlations. 

As expected, all errors decrease with increasing number of bins. Using 
only two bins in cos 6 leads to errors which are 6-7% higher in ApB and (Pt), 
corresponding to a 3-4% higher error in the combined value of sin^ 9^. The 
error on (P^) is unaffected by varying the number of cos6'-bins, as expected. 
10 bins in cos 6' and x seem to be sufficient, since the improvement obtained 
by using 100 bins in each variable is marginal. The number of possible 
bins is of course limited by statistics. However, replacing our least square 
method with the (event by event) ML approach, one should get effectively an 
accuracy corresponding to infinite number of bins. This was also tested by 
performing this ML fit on 50000 Monte Carlo r-pair events. The results are 
listed in Table H^. They turn out to be very close to those of the least square 



CHAPTER 4. THEORETICAL BACKGROUND 



64 







T— 








X 


cos^ 




V 






l)iiis 


l)iiis 




A(P.) 




Asm f/,{ 


10 
10 
100 


2 

10 
100 


1.000 
0.937 
0.935 


4.61 
4.61 
4.58 


4.61 
4.32 
4.28 


0.390 
0.376 
0.373 


ML 


0.936 


4.39 


4.09 


0.359 


T— s>7r(K)z/T- 


X 


cos^ 




V 




Asin^ ef 


bins 


bins 


AAfb 






10 

10 
100 


2 

10 
100 


1.000 

0.937 
0.934 


1.73 

1.73 
1.72 


1.73 

1.62 
1.61 


0.173 

0.168 
0.167 


ML 


0.932 


1.73 


1.62 


0.169 



Table 4.6: The statistical errors in the r-asymmetries as obtained from a 
two-dimensional fit, using the least square method for different numbers of 
bins and the ML method. 

method with 100 x 100 bins, except for the errors in {Pt) and A^^f from the 
T— >£P^i/,--decay channel which are lower. Hence, using a ML fit for this decay 
channel, may reduce the statistical uncertainty in sin^ 9^ even further, to a 
value which is about 10% below the value obtained by the one-dimensional 
fits. This improvement was also obtained at the beginning of this section, 
using fits on real data. 

4.5.2 Events with Two Identified r-Decays 

With a probabihty of 30% to identify a r-decay (to e, ji, 7r{K)), about 9% 
of the events have two identified r's. Since both hehcities are fully anti- 
correlated, their measurement carries less information compared with the 
case of two r's of different events. This effect was neglected in the previous 
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LEP analyses. In order to take the anticorrelation into account, one has to 
analyze these events using the triple differential cross section of Eq. [4.48 



With this expression, the covariance matrix of the three r-asymmetries can 
be estimated from a three-dimensional fit, in an analogous way to the two- 



dimensional fit of the last subsection. The results are listed in Tables ^ 
and |4.8| . Comparing these with the results of events with one identified r 
(Tables |4.4| and [4.5|) , one can see that the errors in (P^) and Ap^f are smaller 
by a factor of from that expected for the case of two uncorrelated r's. 
The errors in Ai?^ are the same in both cases, but the correlations between 
ApB and (P^) or A^^f become now more significant. 



decay modes 




Vn X 








AAfb 






n 














0.937 


3.31 


3.10 




7r(K)z/^ 


0.937 


1.66 


1.56 


7r(K)z/, 


7r(K)z/^ 


0.937 


1.37 


1.28 



Table 4.7: The statistical errors in the r-asymmetries for events where both 
r-decay are identified. 



decay 

Tl 


modes 


PApsiPr) 
(%) 


Pa \fb 

Kfb%oI 

(%) 


P{Pr)A^^S 
(%) 






3.09 


3.84 


3.71 






6.15 


7.69 


0.31 


7r(K)z/^ 




7.48 


9.36 


0.05 



Table 4.8: The statistical correlations in the r-asymmetries for events where 
both r-decay are identified. 



All our results lead to the conclusion that the data should be analyzed 
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using a simultaneous fit to all the relevant kinematical variables which can 
be measured. The only way to realize it, especially for data with limited 
statistics, is the event by event ML method which has the following advan- 
tages: 

• All the information included in the data is used and nothing gets lost 
through binning. 

• All the correlations between the three asymmetries are obtained by 
this method. If one chooses to assume the SM and lepton universality, 
taking sin^ 0^ as the fit variable, all the correlations between the three 
asymmetries are implicitly taken into account. 

• For events which are not exactly on the Z^-peak, the asymmetry values 
vary as a function of s (see Fig. |4^ ). In this case it is advantageous 
to express the asymmetries as a function of sin^ 9^ which can then be 
fitted over the whole line shape. 

• It is possible to use in the same fit events with different information 
content. For example, events where no r-decay mode was identified 
have only the cos^- information. Events where only one r-mode was 
identified will have the cos 6 and x variables. Events with both r-decay 
modes identified will have the cos 6', Xi and xj variables. 



Chapter 5 
Data Selection 



5.1 Data 

This analysis uses r-pair events selected by the Tau Platform (TP) program 
(TP103 version)!] using data taken by OPAL in 1990 (PASS 4) and 1991-1992 
(PASS 5)0. The detector and trigger status requirements are summarized in 
Table |5.1| . They are identical to those required by TP 103 except for the 
status cut on FD which is not imposed here, since the FD is not used in any 
stage of the event selection and the luminosity measured by this detector, is 
not needed in this analysis. 

In Figures |5.1| and |5.2| we compare the data samples of 1992 with the com- 
bined 1990-1-1991 by looking at the distributions of those measured variables 
which are used to select r-pairs and identify their decays. No statistically sig- 
nificant difference was found between the two samples. We therefore combine 

^ The Tau Platform (TP) is an OPAL set of programs that create and manipulate "mini 
DST" (Data Summary Tapes) files, designed for the t decay analysis. The first version 
contributors are: K. Riles, M. Sasaki and E. Etzion. The version used in our analysis is 
TP103. A primer which gives an overview of the programs, its details, the variables and 
the event selection is available for version 1.02 ||5^. 

^PASS4 and PASS5 are stages of the events reconstruction done with ROPE, a collec- 
tion of modules designed to reconstruct events from digits produced by the OPAL detector, 
or by the OPAL simulation program, GOPAL. Its aim is to produce the DST files, and it 
also provides the event viewing facility through its graphic processor GROPE. 



67 



CHAPTER 5. DATA SELECTION 



68 





CV 


CJ 


TB 


PB 


EB 


EE 


HS 


MB 


Detector 


3 


3* 


3t 


2 


3 


3 


3 


3 


Trigger 




2 






2 


3 







* CJ>2 for runs 2001-2293. 
t only for 1991-1992 data. 

Table 5.1: The detector and trigger minimal status requirement. The trigger 
and detector status of each subdetector is coded in the following way: 

- Trigger / detector status unknown. 

1 - Trigger / detector status not on. 

2 - Detector partially on (low voltage or readout problems). 

3 - Detector in a good operation mode. 

them into one sample. 

5.2 Monte Carlo and Control Samples 

In order to simulate the various processes which potentially contribute to 
the selected tau-pair data sample, several Monte Carlo (MC) data sets were 
used. The response of the OPAL detector to the generated particles in each 
case was modelled fully using GOPAL, a simulation program based on 
the CERN GEANT3 package. We have further smeared the particles 
momenta and their electromagnetic energies of these MC events in order to 
account for effects not well modeled in GOPAL. In all cases, the MC and 
real data were reconstructed and analyzed in an identical manner. 

Initial estimates of the efficiencies and purities of the selections for the 
individual tau decay channels were obtained using the KORALZ 3.8 MC 
generator |Q. This takes into account initial state bremsstrahlung (with 
exclusive exponentiation) up to (9(a^), final state bremsstrahlung and elec- 
troweak corrections up to 0{a), It accounts also for single bremsstrahlung 
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Figure 5.1: (a) Comparison of 1990+1991 (histograms) with 1992 data sam- 
ples (error bars). The variables compared here are: ptcl - transverse mo- 
mentum using ECAL clusters, ptch - transverse momentum using charged 
tracks, cost - cos 9, acop - acoplanarity, ncls - no. of ECAL clusters, ntrk - 
no. of charged tracks, nhhc - no. of HCAL strip hits assigned to a cone, xttr 
- total track momentum/E^beam, xtcl - total ECAL cluster momentum/E'^eam- 
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comparison data 90-1 with data 92 
30000 




muon hits 



HCAL hits 



5 10 

HCAL hits (last 3 layers) 



Figure 5.2: (b) Comparison of 1990+1991 (histograms) with 1992 data sam- 
ples (error bars). The variables compared here are: conm - cone's invariant 
mass, xels - ECAL energy excluding the clusters associated to track/E'^eam, 
dphm - 5(f)max, maximum difference between the track and the prcsampler 
signals, eovp - E/P, nb90 - no. of blocks containing 90% of ECAL en- 
ergy, pclc - ECAL energy, MUON hits, HCAL hits, HCAL hits (last 
3 layers). 
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(in the leading logarithmic approximation) in tau decay for the decay modes 
used in this analysis. Tau polarization and its effect on the decay spectra, 
as well as the correlation between the two r's are included (see Sect. §3). 
372K MC tau pair events were generated for these studies, which is a sample 
approximately six times larger than the data sample. 72K events were gen- 
erated with the 1991 detector setup (GOPAL127/ROPE401) and the rest, 
300K, were generated with the 1992 setup (GOPAL129/ROPE403). The 



MC samples illustrated in Figs, and |5^ show no significant difference 
between these two samples. Hence, also the MC samples are combined. 
In order to estimate the "non-tau background" coming from the multi- 



hadronic decays of the Z°, the JETSET MC was used with the parameters 



tuned to fit the global event shape distributions of OPAL MH data . Back- 
grounds from radiative dimuon and Bhabha events were estimated by using 
events generated with the KORALZ 3.8 MC program |^6[ and the BABAMC 



MC program , respectively. Non-resonant t-channel two-photon processes 



were simulated with the generator described in reference ||60|| . 

Other event samples were also used for further corrections in the analysis. 
Some of them were measured data samples of other processes such as //-pairs, 
Bhabha events, single electrons and photon-photon interactions. For this we 
used events recorded in 1990-1992 corresponding to our r-pair sample. The 
events were used as control samples to investigate some corrections needed 
in our analysis, as will be later described. MC simulations of these processes 
were used to test the quality of the detector simulation and to estimate the 
background to our signal. More details about these event samples will be 
given in the next chapters. 

5.3 Tau Pair Selection 

In this analysis we initially select r-pair candidates from which samples of 
T^eUeT^T, 'T^fJ'i^fj.^T and T—^niK)^^- candidates are subsequently identified. 



CHAPTER 5. DATA SELECTION 72 



10000 




Figure 5.3: (a) Comparison of MC samples with 1991 (histograms) and 1992 
(error bars) detector setups. The variables shown are identical to those given 
in Fig. 5.1 
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Figure 5.4: Comparison of MC samples with 1991 (histograms) and 1992 
(error bars) detector setups. The variables shown are identical to those given 
in Fig. 5.2 
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Detailed description of most of the r selection and identification criteria 
can be found in the OPAL Tau Platform manual |Q. At -E'cm=91 GeV, 
e~^e~ —>-T^T~ events have sizeable missing energy and unbalanced trans- 
verse momentum appearing with a topology characterized by a pair of back- 
to-back, narrow jets with low particle multiplicity. In the terminology of r 
recognition we do not restrict ourselves to the conventional "jet" definition, 
in place of that it is used as a general name which stands for the low particle 
multiplicity jets of the r hadronic decays, as well as for the single tracks 
characterizing the r leptonic decay channels. These characteristics are ex- 



ploited by dividing the event into jets as defined in Ref. ||6T| where charged 
tracks and ECAL clusters are assigned to a cone of 35° half-angle . Signals 
occurring in the presampler, HCAL or MUON subdetectors within the solid 
angle defined by a jet cone are also assigned to the jet. Only presampler 
signals associated to an ECAL shower are included. 

Charged tracks and electromagnetic clusters are defined as being "good" 
in this analysis if certain criteria on the quality of the reconstruction are met. 
A good track must have 

• at least 20 hits registered in the jet chamber, the first one at a radius 
smaller than 75 cm, 

• a distance of closest approach to the nominal beam axis in the r-0 plane 
(do) of less than 2 cm, 

• the z at this point {zq) within 20 cm of the nominal interaction point, 
and 

• a momentum transverse to the beam axis of greater than 100 MeV. 

A good barrel ECAL cluster, which is a contiguous group of lead-glass 
blocks registering the deposition of some energy, must have a minimum en- 
ergy of 100 MeV. A good endcap ECAL cluster has more than 200 MeV 
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deposited energy contained within at least two lead-glass blocks, no one of 
which may contribute more than 99% to the cluster's energy. The energy of 
clusters used in the r^ez/g^r analysis is corrected for the expected energy 
loss of electrons in the material in front of the calorimeter. The visible en- 
ergy of the cone is defined as the maximum between the (scalar) sum of the 
charged track momenta and the sum of the ECAL energies associated to the 
cone. 

A I cos^l is calculated for each tau jet and the average of the two jets, 
|cos^|, must be less than 0.68. Background from two-photon processes is 
then suppressed by requiring that the event has a minimum total energy and 
significant missing transverse momentum when the total energy in the event 
is low. After removing cosmic ray backgrounds, the events which remain 
are almost entirely Icpton-pairs. The tau-pair events arc further isolated by 
identifying and removing the mu-pair and Bhabha events using their high 
Kisibie/Ecm characteristics. 

An event is selected as a e^e' — >r+r~ candidate if each of the following 
r (Tj) requirements is met. 

Tl.] In order to reject MH the event is required to contain between one and 
six good charged tracks and no more than 10 ECAL clusters. 

T2. The event must have exactly two tau jets and the average |cos^^| of 
the jets has to be below 0.68. This geometrical requirement, although 
introducing a 42.0% loss in efficiency, removes regions of non-uniform 
calorimeter response where the separation of different tau decay modes 
is degraded. 

T3. Cosmic rays events are removed by requiring that 

(i) at least one track in the event has |cio| and |^o| below 0.5 cm and 
20 cm respectively, 
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(ii) the average value of \zq\ of all the good charged tracks at their 
point of closest approach to the beam is below 20 cm, 

(in) there is at least one TOF signal within 10 ns of the expected 
value assuming particles in the event to originate at the e"'"e~ 
interaction point, 

(iv) the event is discarded if all pairs of TOF signals which are sep- 
arated by more than 165° in the azimuthal plane have a time 
difference of more than 10 ns. 

T4. Events from two-photon processes and events with extreme radiation 
are suppressed by requiring that the acollincarity between the two tau 
jets is below 15°. The two-photon background is further reduced by 
rejecting events if: 

(i) the sum of the visible energies of the jets (defined as the maximum 

of the scalar sum of the charged track momenta and the sum of 
the ECAL energies in the jet) is below 3% of Ecm, 

or, 

(ii) the sum of the visible energies of the tau jets is less than 20% 
of Ecm and the missing transverse momentum of the event is be- 
low 2 GeV. This requirement on transverse momentum must be 
satisfied by both the transverse momentum calculated from the 
charged tracks and that calculated from the ECAL clusters. 

The events passing these T1-T4 criteria are almost entirely leptonic de- 
cays of the Z'^ . The last stage in the tau-pair selection is the identification 
and removal of e+e" ^e+e" and e+e" -^ii'^/i" events through the following 
cuts. 

T5. An event is identified as e+e" — ^e+e" if 
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(i) the total ECAL energy in the event exceeds 80% of Ecm, 
and, 

(a) the total of ECAL energy plus 30% of the scalar sum of the charged 
track momenta in the event is greater than Ecm- 

T6. An event is classified as e"'"e~ if each of the two jets is a muon 

candidate and if the scalar sum of the track momenta in the event plus 
the sum of the energies of the most energetic ECAL cluster in each 
jet is greater than 60% of Ecm- A muon candidate jet must have at 
least one charged track and satisfy one or more of the following three 
requirements: 

(z) there are at least two signals in the MUON chambers which are 
associated to the highest momentum track in the jet, 

(a) there are more than three signals in different HCAL layers in the 
jet, at least one of which is in the last three layers, and the total 
number of HCAL signals divided by the number of HCAL layers 
hit in the jet, '^hiZ/iayer^ is less than two, 

(iii) the ECAL energy of the jet is less than 2 GeV. 

We select 30663 tau-pair events by these criteria. The e"'"e~ — +r+r^ se- 
lection is estimated from MC studies ^ ^ to have 93% efficiency within 
the geometrical acceptance used in our analysis (which is about 54.1% of all 
the Z° T'^T' events) and a background level of 1.7%. These selected 
events are spread over 7 center-of mass energy bins on the Z" peak and 1,2 
or 3 GeV above or below the resonance. The Erm distribution of the 



selected events is plotted in Fig. |5.5| which shows that most of the events are 
measured on the peak bin itself. 

Details of the background in each r-decay channel are discussed in the 
following chapters. Note that the criteria designed to remove two-photon 
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Figure 5.5: A logarithmic plot of the Ecm distribution of the r-pairs selected 
in 1990-92. The numbers on top of each bin present the exact number of 
events. Most of the events (91.5%) are at the central energy bin (the 
peak). 
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(T4), Bhabha (T5) and mu-pair (T6) events introduce kinematic biases which 
must be studied for the measurement of the polarization. The treatment of 
the above selection criteria and the systematic uncertainties associated with 
them are considered in Chapter 

5.4 Tau Decay Identification 

After selecting a r-pair event, an attempt is made to classify the r decay 
channels of its two jets. Our r^euei^r identification algorithm selects jets 
containing a track which deposits almost all of its energy within a highly 
localized region of the ECAL. The r^^v^Vr selection requires an isolated 
charged track with energy deposition in the ECAL and HCAL consistent 
with the passage of a minimum-ionizing particle and signals in the MUON 
chambers associated with the track. 

Unfortunately, the signature for a r— s>7r(K)z/,- decay is not as distinct 
as that for a r^euei^r or a t-^^v^Vt decay because hadronic interactions 
may begin in the magnet coil, ECAL or HCAL. While electrons and muons 
can be efficiently removed from the t^'k{K)vt- sample, the semi-leptonic tau 
decays containing neutral pions are problematic because at Ecm=91 GeV 
the boost of the tau causes significant overlap in the ECAL of the charged 
pion hadronic shower and the electromagnetic showers from the 7r° decay to 
photons. We exploit the fact that electromagnetic interactions begin in the 
coil with a higher probability than hadronic interactions. The presampler 
thereby provides an effective veto against decays containing neutral pions. 

5.4.1 r^eUei^T identification 

Typical candidate for the r^euei^r decays is characterized by a single track 
which deposits its energy in a narrow cone before the end of the ECAL. A tau 
jet is identified as a r^eve^r candidate if each of the following requirements 
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are met. 

El. The jet has one or two good charged tracks. If the jet has two tracks, 
the higher momentum track is taken as the electron candidate. 

E2. The ECAL energy associated to the candidate track divided by the track 
momentum, E/p, is between 0.7 and 2. 

E3. The electrons shower is expected not to be too wide. Hence, it is required 
that at least 90% of the energy of the ECAL cluster associated with 
the candidate track is deposited within three lead-glass blocks. 

E4. The ECAL energy not associated to the candidate track, Eg^cess, is 
less than 4% of the beam energy. This removes tau decays containing 
neutral pions and highly radiative tau pair events. 

E5. The difference between the track and the prcsampler signal furthest 
away in (f> but still assigned to the jet, Scpmax, is less than 5°. This 
effectively removes tau decays containing neutral pions. 

E6. To reduce background coming from muons or from r hadronic decays, 
the number of HCAL layers containing a signal which are assigned to 
the jet is less than two. 

E7. If the opposite tau jet contains a single charged track with a momentum 
greater than 75% of the beam energy then the acoplanarity between 
the two jets must be greater than 0.1°. The acoplanarity is defined 
in the v-cf) plane. This suppresses residual background from Bhabha 
events where one of the two electron energies is mismeasured. 

E8. The variable used for the polarization study of r— >ePeJ^T decays is the 
total electromagnetic energy associated with the jet divided by the 
beam energy, Xe — Eds/ Eheam- E^s is corrected for the energy lost 
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in the material in front of ECAL where the correction uses presampler 
information. The x-parameter used in the fit is expected to be between 
and 1. It can exceed 1 due to finite detector resolution. Values far 
above 1 might happen due to bad track reconstruction which is not 
fully taken into account in our parametrization of the resolution (see 
Sect. |67^ ). Tau candidates with very high x are few, but still may affect 
the result of the ML fit because of their extremely low probability. In 
order to restrict ourselves to a kinematic region with high and well 
understood electron identification efficiency, we require measured Xe to 
be above 0.05 and below 1.1. 

E9. Fiducial requirements are imposed to remove candidates entering ECAL 
regions where the four subassemblies of the calorimeter come in contact. 
Good EB performance is important not only for the measurement of 
the electron energy, but also for electron pion identification using E/p 
cuts, and for Bhabha events rejection. The small gaps in the calorime- 
ter coverage in these contact regions allow some Bhabha events to be 
accepted as r^ez/eZ^r candidates. Hence, a candidate is also removed 
if the track in the opposite jet points to one of these regions. 

The PB subdetector is extremely important for the pion identification, 
but it is used also to correct the electron energy (see below in Sect. |673|) . 
If the track enters a region of the presampler which is not fully opera- 
tional, the tau jet is also removed. 

ElO. The jet is not identified as a r^7r(K)z/T- candidate according to the 
criteria listed below. 

5.4.2 t^^v^Vt identification 

The signature of this decay is an isolated charged track characterized by a 
small energy deposition in the ECAL and HCAL, and with associated signals 
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in the MUON chambers, consistent with Ere of a minimum-ioniziner 

particle. The t^iiu^j^Vt candidate jet must satisfy each of the following re- 
quirements. 

Ml. There is exactly one good charged track in the tau jet. 

M2. At least two of the following three criteria which characterize the pas- 
sage of minimum ionizing particle must be met: 

[i) two or more hits in the MUON chambers are associated to the 
track, 

[ii) the number of HCAL layers containing associated hits, N^^^^, is 
greater than three and ^ hits /layer ^^^^ than two, 

[iii) the ECAL energy in the jet is less than 2 GeV. 

M3. In order to further ensure the minimum ionizing particle characteristic, 
the jet is discarded if N^^^^^ is greater than two and ^hus/iayer greater 
than or equal to three. 

M4. The combined invariant mass of the track (assuming a charged pion) 
and all the ECAL clusters in the jet (assuming 7) is less than 0.3 GeV. 
This is calculated after subtracting 0.5 GeV from the cluster nearest 
to the track in order to account for energy deposition from a minimum 
ionizing particle. This suppresses contamination from semi-leptonic 
tau decays accompanied by neutral pions. 

M5. Residual background from e'*'e~ -^ii'^/i" events is suppressed by re- 
jecting jets in which the opposite jet is consistent with a muon, and 
the sum of track momenta and ECAL energies, excluding the expected 
energy deposition of 0.5 GeV, is greater than 80% of the beam energy. 
The opposite jet is considered to be a muon if: 

(i) one of the first two conditions of cut M2 is satisfied. 
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or if, 

(ii) the track extrapolates to a region not fully covered by any of the 
HCAL or MUON chambers. 

M6. A cut on the variable used for the polarization study of r-^^v^Vr 
events. It is defined as = {ptrk + E^s - 0.5)/Ef,eam where ptrk + E^is 
is the sum of track momentum and ECAL energies associated to the 
jet, and 0.5 GeV is the expected energy deposition of the muon. As 
in the electron case (see cut E8) the muon candidate is required to 
have above 0.05 and below 1.1 in order to ensure a reliable and well 
understood identification efficiency. 

M7. An additional fiducial requirement is imposed to suppress background 
from /i-pair events having tracks close to an anode or cathode wire plane 
of the jet chamber. This source of background is not well simulated. 
A cut rejecting muon tracks near the anode wire planes has been in- 
cluded in the 1990 analysis [0]. Later this cut was removed from the TP 
standard cuts because it was believed that the momentum resolution 
near the anode wires has been improved and the remaining deficiency is 
properly described in GOPAL. However, a recent study ||62| has demon- 
strated that there is still a problem, near the anode wires and even near 
the cathode wires. Due to the deteriorated resolution in those regions, 
/i-pair events which should have been rejected by the cuts T6 and M5, 
still infiltrate to our sample more than expected from the MC. Con- 
sequently, a candidate track that is within 0.5° of an anode plane or 
within 0.3° of a cathode plane is removed if both tau jets are identified 
as muons or if one jet is identified as a muon and the opposite tau jet 
is not identified at all. 

M8. The jet is not accepted as a r^^u^Ur candidate if it also satisfies the 
r-^7r(K)z/,- criteria listed below. 
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5.4.3 r^7r(K)z/T- identification 

The difficulty in the identification of this channel arises from the fact that its 
signature is not highly distinct from those of other channels, and in particular 
hadronic channels with neutral pions. 

A T^7r(K)i/^ candidate must satisfy the following requirements. 

PI. There is exactly one good charged track in the jet. 

P2. The total cluster energy associated to the track is less than 80% of 
the reconstructed track momentum. This removes electrons and other 
contaminants having large deposited energy in ECAL. 

P3. Eexcess IS Icss than 2% of the beam energy, and d(f)max is below 0.5°. 
This rejects semi-leptonic tau decays containing a neutral pion using 
the same variables described in the r^eugi^r selection criteria E4 and 
E5. 

P4. In order to suppress the contamination from t^^VjjVj. decays, only 
tracks extrapolating to a geometrical region in which there is active 
MUON chamber coverage are selected. Furthermore, if a track extrap- 
olates to a region where HCAL is active, the jet is rejected if 

(i) the number of layers hit out of the four MUON chamber layers 
and the last three HCAL layers, '^u^irs^ ^ more than two, 

and or, 

(^«) KayH^s^ = 2 and ^t^s /layer bclow 3. If thc track extrapolates 
to a region which is not fully covered by HCAL then the jet is 
rejected if the number of hit MUON layers is greater than 2. 

P5. The variable used in the polarization study of t^'k{K)i't is Xh = Ptrk/Ebeam- 
For the same arguments given for the electron and muon (cuts E8 and 
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M5) the pion candidate is required to have 0.05 < Xh < l-l- This 
ensures a rehable and well understood identification efficiency. 

P6. Fiducial requirements are imposed to remove candidates entering ECAL 
regions where the eight subassemblies of the calorimeter come in con- 
tact (See E9). Tau jets are also removed if the track enters a region of 
the presampler which is not fully operational. 



5.4.4 T^pUr identification 

The T^pVr channel is not included in the present analysis. We utilize how- 
ever, the results obtained by OPAL in other analyses (see Ref. []T5[ and |T^) 
to combine with our findings. Hence, here we describe only briefly its se- 
lection procedure. Due to the fact the p decays into a charged pion and 
a neutral pion, the signature of the t—^pUt decay is a single charged track 
accompanied by ECAL energy deposition consistent with the interaction of 
the two photons from the neutral pion decay. 

Two analysis methods to identify rho decays were used in OPAL. The 
two t—^pUt- identification algorithms (for Sample I and II) are similar to 
those used to measure the T—*h 7r° branching ratio [Q. The selection of 



sample I is based on an attempt to reconstruct from the ECAL clusters 
the neutral pion which is then combined with the charged track to form a 
p-meson candidate. 

The second r— >pz/^ selection procedure uses a different 7r° identifica- 
tion algorithm. It finds clusters of activity in the ECAL which are restricted 
to a maximum size of 2x2 blocks in 6 and (p plane. Such clusters contain 
99% of the energy of an electromagnetic shower and 95% of the energy of a 
hadronic shower. These clusters which are limited in size differ from those 
used in the r+r~selection or the r— s^ez/eZ^r, T^pv^Vr and r— s>7r(K)z/^ identifi- 
cations. Unlike the sample I selection, this algorithm accepts as 7r° candidates 
single clusters which are not associated to the track. 
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5.4.5 Summary of r Decay Identification 

The typical global efficiencies obtained by the selections described above are 
65.8%, 72.1% and 30.4% for the r— ^ez/gZ/^, r^fiu^u^ and r^7r(K)z/T- , respec- 
tively. These numbers include the approximately 93% preselection efficiency 
within the fiducial region of the analysis and ignore decay or detection cor- 
relations. The corresponding background levels are 5.0%, 2.5% and 7.2%, 
respectively. Plots of the global efficiencies as a function of x are shown in 
Fig. p. 

The results from the r^pur decay are based on two samples of r^pur can- 
didate events. Each sample is independently selected and values for (P^) and 
A^gf are determined separately. Both selections have an efficiency of approx- 
imately 50% and about one third of the combined number of candidates is 
common to both samples. Sample II candidates are selected only from the 
1991 and 1992 data which contain 27353 tau-pair events. The efficiency for 
selecting sample I is 47.5% where the background is 21.2%. For sample II 
the efficiency is 57.5% and the background 22.5%. 

In Table |5.2| we list the number of r-pair events in the data, according 
to their r decay mode identifications. The percentage efficiencies and back- 
ground are also presented for each event classification within | cos^^| < 0.68. 
The efficiencies include the contributions from the tau-pair selection. The 
label '0' refers to the case where the r decay was not identified at all. The 
events in which none of the two jets was identified (23% of the whole sample) 
are not used in the analysis. As is evident from this table, the identification 
efficiencies of tau jets are not independent. The requirements that remove 
two-photon, Bhabha and mu-pair events introduce these correlations. Such 
correlations are, however, taken into account in the analysis. 

Fig. |5.7| shows a "typical" t~^t~ selected event in the detector^. In the 

^ The plot is produced with the event viewing facihty, GROPE, which is the graphic 
processor of the OPAL reconstruction program, ROPE. 
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Figure 5.6: Global efficiencies as a function of x for a) r^eUeUj-h) 
T^liUfj^UT-smd c) r— >7r(K)^'T- decays. The efficiencies include the 93% tau- 
pair selection efficiency and are after the fiducial cut of | cos^| < 0.68. The 
complex behavior in the low and high x regions of a) and b) are caused by 
requirements designed to remove two-photon, Bhabha and mu-pair events. 
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figure we see an e+e~ t^t~ event selected by our criteria. One of tlie tau 
decays into a one cliarged particle while on its back side we find a coUimated 
3 charged particles jet. 
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Table 5.2: Classification of the r-pair events with | cos^| < 0.68 where the 
first figure in each box represent their number. The efficiency (e) and back- 
ground (P) given in percentages are presented underneath. The efficiencies 
include the contributions from the tau-pair selection ~ 93%. The efficiencies 
in the first column (i.e. when only one tau is identified) are subdivided into 
two classes: the first is for the case where the opposite tau does not decay 
into one of the four channels; the second (in brackets) is where the opposite 
tau decays into one of the four channels but is not identified. The label '0' 
refers to the case where the r decay was not identified. 
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Figure 5.7: A display of typical r'^r~ event in the OPAL detector. Tracks in 

the central tracking system are shown in red. Small magenta boxes show hits 
in the TOF system. Clusters of energy in the lead glass ECAL are shown as 
blue boxes, of size proportional to their energy. Similarly, clusters of energy 
in the HCAL are drawn in green. Penetrating charged particle tracks, which 
are candidates for muons, are shown as light-blue arrows. In this event ,one 
of the taus (on the right side) has decayed to a penetrating muon, and the 
other to an electron. 



Chapter 6 
Physics Analysis 



6.1 The Fit Procedure 

For the analysis of the r— *-ez/e^r; ^/U^'^'^r and r-^7r(K)z/T- decays, {Pr) and 
Apgf are determined using an event-by-event ML fit to the data of the the- 
oretical distribution (Eq. ^.371 ). In this method, the following expression is 
minimized, 

where C is the likelihood function. The sum in Eq. |6.1j runs over all N 
selected r-pair events, where at least one of the two r's has been identified 
as a r— ^ez/gZ^T, t^^v^Vt or r— >7r(K)z/T- candidate. The term in the logarithm 
is the differential cross-section for a production r-pairs where one r decays 
into channel i and the second one decays into channel j ( Eq. |4.37|) . It 
is normalized to one, and corrected for radiation, x resolution, efficiency 
and background effects, hence it will be denoted as primed a. This term 
corresponds to the case where both r decays are identified. When only one 
decay is identified, it must be replaced by [1/ a' i)d^a' i/ dcosO dxi. 

When both r decays are identified, the corrected cross-section can be 
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written in the following form, 

X S^{xi,Xj,cos6)[hf{xi)ei{xi,cos6) + (3f{xi,cos6)] 

X [hf{xj)ej{xj,cos9) + P^{xj,cos9)]} 

+ A7"-"(xi,x„cos^). (6.2) 

Here £^{xi, Xj, cos 9) is the efficiency for selecting a r-pair event in which one 
tau decays via channel i and the other via j, ei{xi, cos 9) is the efficiency for 
identifying channel i which for some of the decay combinations depends on 
the decay mode in the opposite jet and on the helicity. The (3f{xi,cos9) is 
the distribution of background events from other r decay channels when ap- 
plying the channel i selection procedure and P^°^~^{xi, Xj, cos 9) corresponds 
to background contributions from events which are not tau-pairs. The func- 
tion hf{xi) is the x distribution for channel i, after applying corrections for 
detector resolution and radiative effects to fi{xi) ± gi{xi). It is given by the 
following convolution function 

hfix,) = /'[/,(x,°) ± (7i(x.°)]rf (x,)7^,(a;,^x,)c^x,°. (6.3) 



where the normalization factor Cj was introduced to satisfy, 







oo 

hf{xi)dxi = 1. (6.4) 



The is the Xj before resolution effects are taken into account. The func- 
tions fi and gi are the Born distributions listed in Chapter ^ Eqs. [4. 30] , [4.33| , [4.35| , |4.36 . 
The correction term rf(xi) is due to radiative effects and TZi{xi^,Xi) is the 
momentum resolution function. 

If only one tau decay is identified, the corrected differential cross-section 
takes a simpler form, 

TT-W - Ej|(l±(P.))(l + cos^.) + |(A,«±A-)cos.] 
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X £^ {xi, cos 9)[h^ {xi)ei{xi, cos 9) + {xi, cos 6)]} 

+ A7"-"(xi,x„cos^), (6.5) 

where 008 6*) is given by the following expression, 

8-'^{xi,cos6) = S^{xi,cos6) (6 
+ Bj I dxj£^{xi, Xj, cos 9)hf{xj)[l — ej{xj, cos 6 

Here, the first term is the efficiency for selecting a r-pair event in which 
one tau decays via channel i and the other via a channel not included in 
our analysis. The second term is the contribution of events where the other 
channel is included in our analysis but is not identified due to inefficiency. 
The coefficients Bj are the r branching fractions. They were taken from the 
OPAL measurements [0 . The summation in Eq. ^]6| is over the three decay 
channels considered in this analysis. 

The normalization factor in Eq. |6.1| , a' ij (or a' i in the case of only one 
identified tau) are the total corrected cross sections evaluated in every step 
of the fitting procedure. 

Our fitted parameters are {Pr) and A^^f . In a former work of OPAL 
the values of Apb{s) have been determined using the distribution of cos^ 
alone. In the present analysis we have used these reported values as an input. 
Since ApB has a strong s dependence (see Fig. ^^ ), we divide the available 
cm. energy range as shown in Fig. |5.5| into 7 intervals. For each interval we 
calculate the average of the OPAL ApB measurements taken during 1989- 
1992. In our global fit we use for each event the corresponding Aps value 
according to its cm. energy. The energy intervals, the average energies, the 
sample sizes and the corresponding A^^^ values are listed in Table |S.1| . The 
uncertainties in the Aps values were taken into account in the calculation of 
the overall systematic error. 

In performing the fit, the functions kf^^.^ e^^j), (3^^^^, S^^--^ and as 
well as the normalization factors are calculated for each event. The hf,-^ func- 
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rn 


( -^crn ) 


#of 




{Gc 


V) 


{GcV) 


cn'Oiits 


(<X) 


87.7- 


88.7 


88.398 


149 


-29.0 ±5.1 


88.7- 


89.7 


89.372 


337 


-8.4 ±4.0 


89.7- 


90.7 


90.228 


555 


-7.8 ±3.2 


90.7- 


91.7 


91.276 


28034 


1.2 ±0.5 


91.7- 


92.7 


92.046 


790 


5.1 ±2.8 


92.7- 


93.7 


93.072 


478 


10.8 ±3.4 


93.7- 


94.7 


93.914 


320 


14.6 ±4.0 



Table 6.1: Forward-backward asymmetries used in the fit and numbers of 
events for 7 cm. energy intervals. 



tions are obtained from the theoretical decay spectra corrected for radiative 
effects using large MC r^r~ event samples generated without detector simu- 
lation and corrected for resolution effects using resolution functions measured 
in data samples of c+c^ — ^c+c^ events, low energy "single electron" events 
from highly radiative Bhabha scattering, and e^e~ ^fi^fi^ events. The 
1^%) ^^"^ functions arc parametrized in terms of simple 

functions, in most cases low order polynomials, using few parameters. The 
parameters arc determined from fits to the pertinent MC samples and sub- 
sequently corrected using appropriate control samples from the data. The 
uncertainties in the parameters, their correlations and the corrections are 
taken into account in the treatment of systematic errors which arc discussed 
in detail below. For this we are using the standard error propagation rule, 



AA 



where A is the final result ((Pr) or A^^f ), pi, pj are the parameters used to 
describe the correction function and pij is the correlation matrix between 
them. The partial derivatives ^ are calculated numerically, by modifying 
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the parameter by its uncertainty, by repeating our fit and calculating the 
change in the final result A. 

As suggested in Sect. |4.^ in order to check the parameters characterizing 
the structure of the charged weak decay of the r one has also to treat them 
as free parameters. In that case we use Eqs. |4.44| and |4.48| to modify Eq. ^]2| 
in the following way: 



a 



dcos 9 dxi dxj 



'-^ — =lix{ /.+/.;[c(i+a.)-ne. + e,)] 



+h-hj[Cil+^,Q+V{^, + Q] 
+hthj[Cil-^,^,)-Pi^,-Q] 
+hthj[Cil-^,.^,)+Vi^.-Q] (6.8) 
+2(p - 0.75) [/^(0i(C -V)+ ct>'-{C + V) 

+m<l>X{c-v) + <pi{c + v)] 

+2(6 - 0.75) [e.^7^(0i(C - P) - 0i(C + V) 
Hj9imC-V)-<j>l{C + V)]} 



where C, V, hf^j^^ and 0^^.^ are defined as following: 

C = 1 + COs2 6' + |Ai;'ijCOs6' 

V = (P,)(l + cos2^) + |A^„fcos^ 

hf = fl{x,) + ^n{x,)±[gl{x,) + ^g^,{x.)] (6.9) 
0i = flixi) + pfi{xi) ± ii[g\{xi) + 5g\{xi)] 



In the first two terms of Eq. |6]^ (first two lines) there is an exact repetition 
on Eq. [4.37| where aV — A structure is assumed for the r decay. These terms 



are corrected for radiation and detector (resolution, efficiency and purity) 
effects as in Eq. For the other six terms, which vanish in the case when 
V — A is exact, we cannot use the standard correction procedure as our 
calculated correction are helicity dependent. Therefore in order to correct 



CHAPTER 6. PHYSICS ANALYSIS 



95 



these terms we used corrections averaged over the two hehcity states. For 
example, rj(xj) = 0.5 ■ [rf{xi) + r^^xi)]. 

Eqs. and as it is written is apphcable to the r— >£z/^z/^ events 
through which and 6. In the case of r— i>7r(K)i/^ decay only ^had can be 
measured. In order to generalize Eq. |6.8| and |6.9| we equate for the hadronic 
decays fig)^""^ = f {g)i"''^ + f [9)2""^ where their values, not to change Eq. |4.33 
relations, are: 




^2-^ (6.10) 

1 — mji/m^ 

2a; — 1 — m\/m^ 
(1 - ml/mlY 

6.2 Radiative and Threshold Effects 

For a given helicity state, initial and final state radiation as well as radiation 
in the decay of the tau affect the spectra of the tau decay products. The 
radiative distortions to the spectrum for the decay mode i are taken into 
account using a radiative correction function, rf{xi). 

Radiative effects are usually separated into photonic (QED) and weak 
interaction effects. The QED effects are due to all the diagrams involving 
real or virtual photons. They do not depend on unknown parameters but 
they usually depend on the experimental cuts. The other effects are due to 
the weak corrections and depend on unknown theoretical parameters but not 
on the experimental cuts. 

Since the r-polarization results are part of the TP OPAL lineshape pa- 
rameters and are used with them in all sorts of overall electroweak fits, they 
are handled in the same way. Therefore, following the procedure used in all 
LEP experiments for the line shape analysis, we do not apply corrections due 



fhad 
J2 

had 
92 

/had 
1 



had 

9l — 
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to QED and weak radiative effects on the measured (P^) or A^^f . However, 
even for a given lielicity state, initial and final state radiation modify the 
spectrum of outgoing leptons and pions. 

The functions rf(xi) are determined from the ratio of the spectrum con- 
taining radiative effects to the spectrum where only the Born level cross- 
section is considered. A high statistics run of the KORALZ 4.0 MC genera- 
toif], is used to create the radiative corrected spectrum. 

The QED effects, and in particular the initial state radiation, depend on 
the cm. energy. We use in this analysis data taken on the Z° peak as well 
as at 1, 2 and 3 GeV below and above the peak (see Table |6.1| ). Therefore 
the KORALZ 4.0 is used to generate events distributed at the same seven 
cm. energies as in the data. We actually generated three samples of 420000 
events each, corresponding to the three decay channels investigated here. In 
each sample, all the r's were forced to decay into a single channel. 

For the r-polarization measurement on the Z° , the most important QED 
effect is the final state radiation from the r or its decay products . This 
would modify the decayed particle momentum and distort its distribution 
from the Born term expectation. In order to minimize this effect in our anal- 
ysis, the electromagnetic calorimeter is used to detect the radiated photons 
in a 35° cone around the track. The energy of these photons is added to 
the track energy and in this way, the effect of the final state radiation is 
minimized. This procedure is done only for lepton final states, since in the 
case of pions, events with extra photons are rejected anyhow, because they 
resemble r^pVj. decays. The same treatment is done also for the generated 
events in order to investigate the residual effect of radiative corrections. 



^The KORALZ 4.0 is the latest version of this MC program in which the radiative 
corrections were improved compared to previous versions. It was used to calculate our 
correction to the data. However, as all our other MC tests were done with KORALZ 3.8, 
we have also used radiative corrections calculated with the 3.8 version when we have cross 
checked the analysis with MC events. More details about KORALZ are given in Sect. 5.2 
and Ref. M] 
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Figure 6.1: x distribution of the ratio between generated MC events with 
radiative corrections and the Born term expectation, plotted with our 
parametrization for the radiative effects (solid lines). Figs, a) and b) are for 
the T— >ePei^T channel for positive and negative helicity states respectively. 
Figs, c) and d) correspond to the r— >//P^i/,- channel and Figs, e) and f) are 
for the T^7r(K)i/,- channel. 
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We calculated separately the radiative corrections, rf (xj), separately for 
positive and negative helicity events. In this way the corrections are inde- 
pendent on the polarization measurement. The MC distributions are divided 
by the corresponding lowest order diagram expectations. The obtained ra- 
tios are the radiative corrections rf(xi) which are plotted in Fig. These 
are slightly dependent on the helicity state. The reason of this dependence 
originates partly from the different slopes of the Born spectra. For instance, 
in the channel r— >-7r(K)z/^ the spectrum is increasing (decreasing) for positive 
(negative) helicity. The initial and final state radiation tend to decrease the 
energy of the outgoing pion. Therefore, one should expect the low x range 
be enhanced and distorts the distribution. 

The solid lines in Fig. |6.1| represent our parametrization. We needed 3 
to 4 parameters to describe each distributions. The values between the 
distributions and their parametrization are indicated in the figures. The 



first bin in Figs. |6.1| c and |6.1| d is low due to the finite /i-mass effect, which 
will be discussed in the following, and was not incorporated in the Born 
expression. Therefore, this bin was not included in the calculation. For 
the r— i>'7r(K)z/T- channel the first bin starts at the kinematical threshold and 
therefore no depletion is presented. 

The major uncertainty associated with these corrections arises from the 
fact that the 0{a) in the QED corrections to tau decay are included only in 
the leading logarithmic approximation. The effect of this correction on the 
polarization asymmetries is quantified by analyzing the MC event sample 
after full simulation as if it was measured data. The event sample is analyzed 
twice, once for the full event sample and also for a sample which excludes 
events with radiated photons in the tau decay process. The difference in the 
results is used as a measure of the magnitude of the correction. Following 
Ref. [^], the correction is assumed to have a systematic error of the order of 
l/ln(m^/mj) of the correction. An additional but small contribution to 
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the error on rf(xi) arises from the parameters used to describe the radiative 
corrections due to the MC statistics which is large but nevertheless limited. 
These parameters are obtained, together with their correlations, from the 
fits of the parametrization functions. Their effect on the final (Pr) and A^^f 
results was calculated with the method described in the last section. 



Table summarizes the systematic errors associated with the radiative 
corrections for the various decay channels and for the whole sample. 



source 


all 




HPr) ^Kol 


r^7r(K)i^T- 


MC 
LL 


0.15 0.01 
0.25 0.03 


0.44 0.06 
0.16 0.03 


0.56 0.04 
0.15 0.04 


0.14 0.01 
0.33 0.02 


Total 


0.29 0.03 


0.47 0.07 


0.58 0.06 


0.36 0.02 



Table 6.2: Systematic errors (in %) associated with radiative effects. MC 
is the systematic error introduced by the MC finite statistics and LL is the 
main contribution coming from the leading logarithm approximation. Their 
quadratic sum is given in the last line. 



As can be seen the treatment of initial and final state radiation in KO- 
RALZ 4.0 introduces a smaller contribution to the error than the leading 
logarithm approximation. The total error assigned to the radiative correc- 
tion uncertainties of {Pr) and A^^f are 0.29% and 0.03%. 

We also investigated the possible influence of the radiative corrections on 
the cos 6' distributions. Since we found this to be negligible we have further 
taken rf{xi) to be independent of cos6'. 

To account for the finite muon mass in the r^^v^v-r channel we have 
multiplied the expression under the integral of Eq. |6.3| by the term K'^i^x) 
defined as ^ 

K^{x) = [^-{'^Y forx>x± ^g_^^) 
[ for X < 
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The Xm and q parameters are given in Table |6.3| . They were derived from 
a large MC event sample. Fig. p]^ shows the x distribution of these MC 
events at the low x region together with our parametrization. 



helicity 






(+) 
(-) 


0.0024 
0.0047 


1.27 
1.56 



Table 6.3: Parametrization of the finite /i-mass effect. 
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Figure 6.2: Low x distributions of positive and negative helicity T^fiu^Ur 
events obtained from high statistics MC sample (histograms), compared with 
our parametrization (solid line) of the finite /i mass effect. 



The effect of the threshold for tau decays to kaons (Eq. [4.34|) introduces an 
error on (Pr) and A^^f associated with the uncertainties in the tau branching 
ratios to pions and kaons and leads to errors below 0.1%. There is also a very 
small error (less than 0.1%) introduced from the measurement uncertainty 
of ApB which enters the analysis (see Eq. [4.27|) . 
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6.3 Detector Resolution 

As in every experimental observation the measured quantities have a natural 
spread which depends on the detector resolution. The main influence of 
this in our analysis comes from the resolution in x which is essential for 
the polarization measurement. We are less effected by the inaccuracy in the 
measurement of cos 9 since our fit to A^^f is sensitive mainly to the numbers 
of forward and backward scattered particles and not to the specific shape 
of the cos ^-distribution. Our results are also not effected by the charge 
misassignment effect, which occurs only in a very small number of events. 

A non-neghgible effect on (P^) and A^^f can rise from the x-smearing, 
since the fit is sensitive to very small deviations from the true shape of the 
x-distribution. In particular, due to this smearing, some particles may have 
an unphysical x value above one, which requires a special care. 

To describe the resolution effect we introduce the x resolution function, 
TZi{xi^,Xi), which is for a given x° the measured x. This function depends 
on the identification of the particle. 

For pions, x is the particle momentum normalized to the beam energy. 
For muons, we add to the particle momentum also the energy deposited in the 
EM calorimeter associated with the r-jet. From this EM energy we subtract 
0.5 GeV, which is the expected energy deposited by the muon itself. This is 
done in order to include the energy of a photon which might be radiated from 
the tau or the muon. This correction is very small in most events, therefore, 
also for muons the resolution in x is essentially the momentum resolution. 

For electrons however, we are using only the EM calorimeter to define x, 
summing up all the clusters in the r-jet. Therefore, the resolution function 
for electrons is completely different than that of the // or tt. Both resolution 
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functions, for e, jj, or vr, are investigated in two steps. In the first step we 
determine a characteristic variable y, which takes into account the resolution 
dependence on x and possibly also cos 6*. In the second step we are using a 
control sample to investigate the exact shape of the y distribution. In the 
following we shall discuss each resolution function separately. 

6.3.1 The Resolution Function for Muons and Pions 

First we examine only the shape of Ax/x where as before x = Ptrack/Ebeam 
and Ax is the difference between the true and measured normalized momen- 
tum, X — x^, for each decay particle respectively. This was investigated by 
M. Sasaki using r^^u^Ur and r^7r(K)z/T- MC events which were ini- 
tially generated by KORALZ 3.8, and then processed by GOPAL129 and 
ROPE403. For this study a Gaussian distribution of the kind 

P{x,x) = ^^ e"2^ (6.12) 

v27rcr(x) 

was taken allowing also a possible shift in the average x parametrized as 

X = ai + a2X° + a3(x°)^ (6.13) 

The distribution of x{x^) was investigated separately in 40 intervals between 
to 1. Neglecting the non Gaussian tails the x shift and the width (t{x) 
were determined for each x interval. The relation between and the track 

X 

momentum was taken to be 

^ = V^i + bblack (6.14) 



and from this fitted values for bi and 62 could be obtained (see Table |6.3.1|) . 



Fig. |6.3| illustrates the peak (x — x) shift and the Gaussian width as a func- 
tion of together with the fitted curves estimated from the MC. For the 
T^^v^Vj. decays, x^ has a more complicated definition, namely 

Vt^„^h — E^i^ — 0.5 GeV 
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Figure 6.3: MC T^7r{K)i'r resolution parametrization with the peak {x — x) 
shift (top) and the resolution ax/x (bottom). The solid lines present the fit 
curves. 
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ai 


0-2 


as 


h 


b2 


r— i>7r(K)i/^ 


0.0029 
0.0000 


1.0000 
1.0000 


0.0000 
0.0000 


1.0610 
0.0131 


1.1270 
0.0012 



Table 6.4: The parameters for the detector response function determined 
from T^fiu^UT- and r-^7r(K)z/T- MC investigation. 

where depends on both the E^is and the track momentum. Therefore the 
resolution function for the r—^fiu^Uj. is a combination of the ECAL energy and 
the charged track momentum resolution. This combination introduces tuning 
parameters to weight the components assuming that this decay deposits the 
minimum ionizing energy of 0.5 GeV in the ECAL. Hence, the width term 
has the form 

Cr{x^.) = 7^ \E,t,=0.5GeV (6.16) 

-L-'beam 

Fig. |6.4| illustrates {x — x) shift and the Gaussian width as a function of , 
together with the fitted curves estimated from the MC. Table |6.3.1| presents 



the results of the fitted parameters. As can be seen, a shift in the peak 
appears only in the r^/iz/^i/^ decays. This means that the energy of 0.5 GeV 
assumed to be deposit by the muon in the ECAL was underestimated. 

To obtain the final resolution function for muons and pions we use real 
yU-pair data events to investigate the exact shape of the y^^^ distribution 
[yfi,Tr = (x — x°)/cr(x)]. These events were selected according to the standard 



criteria of the TP103 package [53| with two additional requirements: 



The event acollinearity must be below 0.2°. This cut was imposed in 
order to remove events with initial state radiation where the outgoing 
muons have less than the beam energy. 

The track must be away from an anode (cathode) wire plane by at least 
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Figure 6.4: MC t^ixv^Vt resolution parametrization with the peak {x — x) 
shift (top) and the resolution (Tx/x (bottom). The solid lines present the fit 
curves. 
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0.3° (0.5°). 



Fig. and ^]B| show the resulting and distributions respectively. 
These distributions contain non-Gaussian tails. Therefore, they are parametrized 
as sums of three Gaussians with different width, mean and normalization pa- 
rameters 



7^(x°, X, cos 0) = '^pi exp 



i=l 



y - y-; ^ ^ 



(6.17) 



The values of the fit parameters and their errors are listed in Tables |0 
and ^]6| . The / D .0 .F. values of the fit are 61/53 and 54/64 for the 
t^'k{K)vj. and t-^^v^Vj. respectively. 



i 


1 


2 


3 


norm. 


Pi 


1. 


0.292±0.031 


0.016 ± 0.012 


mean. 




-0.023±0.020 


0.528±0.049 


0.51±0.11 


width. 


<yi 


0.818±0.010 


1.22±0.072 


2.57±0.44 



Table 6.5: Parameters of the three Gaussians describing the exact shape of 
the pion x resolution function. 



i 


1 


2 


3 


norm., pi 


1. 


0.2732±0.0048 


0.0098±0.0037 


mean, |/j 


-0.057±0.010 


0.427±0.030 


0.36±0.14 


width, (Tj 


0.7266±0.0071 


1.127±0.029 


2.60±0. 



Table 6.6: Parameters of the three Gaussians describing the exact shape of 
the muon x resolution function. 



We considered the following contributions to the systematic error associ- 
ated with the muon and pion resolution function. 
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Figure 6.5: The number of events as a function of given in a logarith- 
mic scale. The solid line is our parametrization which represent the energy 
resolution 71. 




Figure 6.6: The number of events as a function of given in a logarith- 
mic scale. The solid line is our parametrization which represent the energy 
resolution IZ. 
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Mu-Pair (data) Momeiitum vs. cost 




s vs. cost. 



Figure 6.7: x vs. 008 6* for data /x-pairs. 
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The uncertainties in the parameters of the resolution functions coming 
from the hmited statistics of the /i-pair sample. We used the errors 
listed in Tables |6.5| and |6.(j| including their correlations in order to 



propagate them to the final results. 

We investigated a possible dependence of the momentum measurement 
on cos6'. Such dependence could result from a possible twist of the 
two end-plates of the jet chamber, that will result in a tracks distortion 
so that the measured momentum would be larger or smaller than the 
true value, depending on the charge and cos^^ values. In this case, 
there would be a net bias in the A^^f measurement. Fig. shows a 
scatter plot of x vs. cos 6 for /i-pairs data which looks uniform, dividing 
the cos 6' range into 14 intervals, average over x for each interval, and 
fitting X to a(l + 6cos6') we obtain b = (-1.99 ± 0.94) ■ 10"^ with 
X^/ D.O. F.=11.8/12. We introduced this small b value in our ML fit 
the effect on the final results was found to be very small (see Table | 



We investigated whether there is any dependence of the y^^T^ distribu- 
tion on I cos^l, by considering separate distributions from MC r^fiu^UT- 
decays for | cos 6*1 below and above 0.35. The obtained RMS values 
were 0.0716 and 0.0759. Assuming a linear dependence of the RMS 
on I cos 6*1, we multiplied y^^^r by 0.996 + 0.011| cos 6*1 and the obtained 
deviations in (Pt) and A^^f are listed in Table 



Table summarizes all the systematic errors associated with the muon 
and pion momentum resolutions. This amounts a total of 0.38% and 0.21% 
for {Pr) and A^^f , respectively. 

Before concluding this sub-section, one should mention two other sources 
of systematic errors associated with the track measurement, which are related 
to the Apof determination. 

• Uncertainty in cos 6 measurement may affect the A^^f determination, 
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source 


whole sample 






/x-pair statistics 
X — cos 6 dep. 
resolution cos 6* dep. 


0.10 0.01 
0.21 0.21 
0.30 0.02 


0.47 0.03 
1.03 1.17 
0.98 0.11 


0.10 0.01 
0.32 0.72 
0.11 0.00 


Total 


0.38 0.21 


1.50 1.18 


0.35 0.72 



Table 6.7: Systematic errors (in %) associated with the muon and pion res- 
olution function. 



since events near cos ^ = can migrate from positive to negative cos 9 
value and vice versa. It also introduces an inaccuracy in the cut on 
I cos 6'| at 0.68. The | cos 9\ value is obtained in our analysis by averaging 
the corresponding values of the two r-jets, taking into account tracks 
and EB clusters. This is a measure of the | cos^| of the tau and not 
just the track. Therefore, part of the inaccuracy in cos 9 results from 
the invisible i/'s in the event. We investigated the cos 9 resolution using 
MC events and obtained an uncertainty of 0.0225. We introduced in 
our ML fit a resolution function in cos 9 using an error in cos 9 twice 
as large, in order to account for possible data-MC differences. The 
resulting variation in A^^f is smaller than 0.05%. 

• Wrong charge assignment of the r-jcts, resulting in a wrong sign of 
cos^. The charge assignment is done by the following algorithm: The 
charge of each r-jet is defined as the sum of its track charges. If both 
r-jets have opposite charges, there is no ambiguity. If one jet charge 
is zero, only the second jet charge is used. If both jets have the same 
charge and this charge is not zero, then the correct charge is taken from 
the jet with a smaller number of tracks. In the remaining cases where 
the two jets have the same number of tracks, we consider only the most 
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energetic track from each r-jet and take the charge of the track with 
the smaller momentum. As a result of this procedure, only 0.35% of 
the events have a wrong charge assignment, as obtained from the MC 
sample. This would reduce the value of A^^f by a negligible amount of 
0.007 X . If we assume that in the data the situation is twice as 
bad, we end up with a systematic error in A^^f of 0.16%. 

In total, the systematic error in A^^f associated with these two sources is 
0.17%. This uncertainty is common to all decay channels. 

6.3.2 The Resolution Function for Electrons 

The electron energy is taken from the EM calorimeter. We are using the 
corrected EB energy and we include a further correction depending on the 
total presampler multiplicity in the cone. In order to investigate the possibil- 
ity of using the PB data, we used Bhabha events and tested the EB energy 
normalized by the beam energy [xeb) as a function of the PB multiplicity 
(number of hits). A decrease in the EB energy for increasing PB multiplicity 
is found as expected due to electrons which started showering in the material 
in front of The result of our study is following correction to Xe, 

X, = + 0.02 X Max{NpB - 10, O)/^^^™, (6.18) 

where Npg is the PB multiplicity. 

In order to investigate the energy resolution dependence on its energy, 
we are using real Bhabha electrons data, single electron events and r— ^ei^e^T 
decays. For r— i>ei>eZ^r electrons we are using our analysis sample and the 



other events were selected using the standard cuts of the TP103 package [§3 
For Bhabha events we required also an acollinearity below 0.2°, to reduce 
radiative effects similarly to our procedure with /i-pairs. 

The Bhabha electrons events are assumed to have the beam energy. We 
do not know the real energy of the electrons from the other sources but we 
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can utilize their E/p and Ap/p distributions at different energies in order to 
obtain the energy resolution, 



|A,£)/,£)|._(^). (6.19) 

X E \ p p p 

The width A{E/p) is obtained by a Gaussian fit to the distribution of E/p 
at the different energy points. For the evaluation of Ap/p we use Ap/p = 



C^x^ + (0.3x)2 + 0.02x/x where the factor C = 0.073 ± 0.006 is obtained 
from Bhabha and /i-pair events. 

Fig. |6.8| presents the AE/E dependence on x. The resolution in the last 
point is estimated from the Bhabha events, whereas the first three points, 
were calculated with the single electrons. As a consistency check we also 
plot the results of the electrons from r^eVg.Vj. decays, using their E/p dis- 
tribution. The error bars of all these results derived from the statistical 
uncertainties of the E/p widths and the error factor C. 

The ^ values were fitted by the function a + h/ y/x, 

Ax AE 

= -— = (0.0192 ± 0.0008) + (0.0139 ± 0.0008)/v^. (6.20) 

X E 

This is shown in Fig. having a y/^ / D.O.F. of 4.7 for 6 D.O.F. 



As for the muons and pions, we define Ue for electrons as, 

= ^-^ =. (6.21) 

0.0192x0 + 0.0139^x0 

As before also here the second step is the study of the exact resolution 
shape, including the tails of the distribution. For the r— ^ez/gZ^r it is investi- 
gated by using Bhabha events describing them by a sum of three Gaussians 
as in Eq. |6.17 . 



Fig. |6.9| presents the ye distribution and our parametrization. The values 



of the fit parameters and their errors are listed in Table 6^. The 'x^ / D.O.F. 
of the fit is 74/58. 
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Figure 6.8: AE/E dependence on x. The last point is calculated with 
Bhabha events and in the first three points it is done with single electrons. 
For consistency we plot also the E/p of the r— ^ePgi^r events from our data 
sample. The line represents our parametrization to this dependence. 
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Figure 6.9: The number of events as a function of i/e given in a logarith- 
mic scale. The solid line is our parametrization which represent the energy 
resolution 71. 
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i 


1 


2 


3 


normalization, pi 


1 


0.192 ± 0.032 


0.0067 ± 0.0018 


mean, yi 


0.120 ± 0.010 


0.00 ± 0.039 


-1.72 ± 0.52 


width, (Tj 


0.571 ± 0.016 


1.317 ± 0.067 


4.65 ± 0.67 



Table 6.8: Parameters of the three Gaussians describing the electron energy 
resolution functions. 



Two sources of uncertainties contribute to the systematic error associated 
with the electron energy resolution. 

• The errors on the energy dependence of the resolution (see Eq. |6.20|) . 

• The errors on the parameters describing the exact shape of the resolu- 



tion function, shown in Table 6.8 



These errors were propagated properly to the final results, taking also their 
correlations into account. 

A summary of the systematic errors associated with the electron energy 
resolution is presented in Table |6.9| . The uncertainty in the calorimeter re- 
sponse introduces total errors of 0.48% for {Pj.) and 0.12% for A^^f . 



source 


whole sample 
A(P.) AA^,f 


r^ez/gZ/T- 


energy dependence 
resolution function 
energy scale 


0.01 0.00 
0.32 0.10 
0.37 0.07 


0.04 0.01 
1.46 0.45 
2.23 0.17 


Total 


0.48 0.12 


2.67 0.48 



Table 6.9: Systematic errors (in %) associated with electron resolution func- 
tion. 
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6.4 Tau-Pair Selection Efficiency 

The efficiency of r-pair selection, £f(xi,cos9), is investigated using the r- 
pair MC events and comparing the kinematic variable distributions before 
and after the selection requirements. This is done separately for events with 
positive and negative r~ helicities and for the various combinations of the 
decays of the two taus. 

When both r's decay via one of the modes considered here (e, fi or 7r(K)), 
the efficiency is parametrized as a function of the kinematic variables Xj, Xj 
and cos 6'. When the selection requirements involve two or even all three 
kinematic variables in a correlated manner then these are explicitly taken 
into account in the parametrization of the efficiency. 



The efficiency used in our fit (Eqs. |6.2| and |6.5| ) is expressed as a function of 
the measfirec? variables, Xj and Xj, on which all the cuts are applied. However, 
we cannot obtain the efficiency as a function of the measured variables since 
for events which were not identified as r-pairs, very often those variables are 
not defined. Therefore, we have to use the corresponding generated variables 
Xi^ and Xj^. As long as the variation of efficiency with the variables is weak, 
it should not make any difference. However, for those r-pair selection cuts 
which depend on xt and xj variable one has to be more careful and include 
this fact in the efficiency function. The three cuts which depend on Xj and 
Xj are T4, T5 and T6 and these are treated in the following way: 

• The efficiency X'i^ X j tj COS 6) to select r-pair events where both taus 
decay into the same channel is set to be if Xj + Xj < 0.4 and (xj — 
Xj) sm9 < 0.04. This corresponds to the cut against two photon events 
(T4). This cut rejects events with low visible energy namely, if both 
missing transverse momentum calculated from the tracks and the from 
the EM clusters are below 2 GeV . It therefore affect mainly r-pair 
events where both taus decay to the same channel. 
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• The efficiency £ij{xi,Xj, cos 9) to select events where both taus decay 
into electron is required to vanish when Xi + Xj > 1.6. This is due 
to the rejection of Bhabha events with cut T5 which influence r-pair 
events of this type. 

• The rejection of mu-pair events by cut T6 affects mainly events where 
both taus decay into muons. For such events, we require that 
S^{xi,Xj,cos9) — when Xi + xj > 1.2. A possible residual effect of 
this cut on events with one T^/iP/j^i/r and one T^7r(K)i/^ decays was 
taken into account by allowing a step function change in £^{xi, Xj, cos 9) 

B;t ~\~ X j ^— X . 2 . 

To determine the r-pair selection efficiency we investigate r-pair MC 
events which satisfy the three cuts mentioned above. We consider the dis- 
tributions of the kinematic parameters before and after the r-pair selection. 
The ratio between the distributions after and before the selection gives the 
efficiency of the selection excluding the cuts T4, T5 and T6. This procedure 
is done separately for positive and negative r~ helicity events. The |cos^| 
dependence of the efficiency turns out to be constant over the whole region 
from to 0.68 and that for all decay channel combinations. 

In the case where both taus decay into electron, muon or pion, we have to 
consider in addition to | cos^| the two kinematic variables, Xi and xj. There- 
fore, we consider two-dimensional efficiency distributions in a 6x6 regions in 
Xi and Xj and fit them using expression of the following form 

£ij{xi, Xj) = ai(l -|- a2Xi + 04X4^) (1 -|- a^Xj + a^Xj'^) 

+aQXiXj + a-jXi'Xj + a^XiXj^ + a^Xi'Xj' . (6.22) 

For Xi < .05 or Xj < .05 the efficiency is reduced by a factor, while if both Xi 
and Xj are < .05 then the efficiency is set to 0. The parameter a\ is expected 
to be close to 1 and all other parameters are expected to be close to 0. If after 
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doing the fit any of tlie parameters turns out to be consistent within error 
with its expected value it is set to that value as a constant in the fit, and 
the fit is repeated. In this way we are using just the number of parameters 
necessary to obtain a reasonable fit. 

When only one r decays to e, fi or 7r(K) and the other r decay is not 
identified, the selection efficiency has two components. The first accounts for 
the case where the other r does not decay via e, fi or 7r(K) while the second 
accounts for the case where it does decay via these channels but could not 
be identified. The first component is parametrized as a function of Xj and 
cos^, checking that these two kinematic variables are not correlated. The 
second component, which arises from identification inefficiencies, can have 
correlations which were taken into account. 

In the cases where only one r decay is identified, the selection efficiency 
Si{xi) is parametrized as a quadratic function of just one variable, Xj. For 
Xi < .05 the efficiency, Si{xi), is further reduced by a factor. Also here, when 
any parameter is obtained from the fit as consistent with its expected value 
(1 for the normalization and for the others) it is set to the expected value 
and the fit is redone. 

In Fig. |6.10| we plot the r-pair selection efficiencies as function of x 
for these cases where only one tau decays to e, fi or 7t{K) together with 
our parametrization. In Table |6.1CI| we list the values obtained in our 
parametrization fits. 

The error coming from the r-pair selection efficiency is composed of the 
MC statistics and the error related to the resolution and scale uncertainties 
in the ECAL and tracking detector. 

Systematic uncertainties in the efficiency resulting from the MC statis- 
tics were investigated by propagating the statistical errors in the parameters 
obtained from the parametrization fits to the final results as explained in 
Sect. |6TT| . 

As for possible effects in the data not fully modeled by the MC, we were 
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Figure 6.10: Tau-pair selection efficiencies as functions of x for events where 
only one tau decays into electron (Figs. a,b), muon (Figs. c,d) or pion (Figs. 
e,f). Figs. a,c,e (b,d,f) correspond to positive (negative) helicity events. 
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Decay channel combination 


positive helicity 


negative helicity 


a) T^eUe^T, none 


13.6/16 


15.2/16 


b) T—^fiUfj^Ur, none 


12.6/15 


12.6/14 


c) r— >-7r(K)z/T-, none 


11.4/14 


9.0/14 


d) T—^eUel^T, T— ^ez/gi^r 


15.9/29 


24.9/32 


e) T^eUe^r, T-^filJf,Ur 


35.9/24 


41.6/24 


f) T—>-eUel^T, T— >7r(K)z/T- 


18.9/27 


41.7/24 


g) r^/iZ/^Z/^, T^^V^Vr 


19.4/29 


29.4/32 


h) T^iiu^Ur, r^7r(K)z/^ 


32.4/25 


44.9/24 


i) r— i>7r(K)i/T-, r— s>7r(K)z/T- 


12.1/33 


20.0/28 



Table 6.10: 'x^ jD.O.F. coming from the parametrization fits for r-pair selec- 
tion efficiency. 

concerned with the momentum and energy resolutions being broader in the 
real data as compared with the MC. In order to investigate this effect on the 
final results we used our MC event sample and smeared the generated track 
momenta and electron/photon energies according to our resolution functions 
obtained from real data control samples (see Sect. After replacing those 
smeared values with the corresponding MC simulated ones, we recalculated 
the r-pair selection efficiencies and parametrized them in the same way as 
before. The different tau selection efficiency parametrization, calculated with 
and without additional smear to the MC, has a negligible effect to the overall 
uncertainty in (P^) and in A^^f of the order of 0.03% and 0.05% respectively. 
Nevertheless, we have used in the analysis the MC sample with an additional 
smear. 

Table |6.11| summarizes the systematic errors associated with the r-pair 
selection efficiency for each decay channel separately. 

Thus, the contributions from the tau selection efficiency to the overall 
uncertainty on (P,-) and A£,f are 0.39% and 0.19%, respectively. 
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A(P.) AAf« 


r— s>7r(K)i/T- 
A(P.) AA^,f 


error 


1.10 0.80 


0.59 0.06 


0.36 0.08 



Table 6.11: Systematic errors (in %) associated with r-pair selection effi- 
ciency. This is essentially based on the MC statistics. 

6.5 Efficiencies of the r Decay Identification 

The efficiencies to identify the various r decays arc obtained in two steps. 
In the ffist step, the MC events are used in order to obtain the efficiency 
distributions for each decay channel, separately for events with positive or 
negative r~ helicity. These distributions depend on the x variables of the 
identified r and on cos 6. However, some of the identification requirements, 
such as those designed to remove Bhabha and mu-pair events, also introduce 
a dependence on the x variable of the r on the opposite side. This dependence 
is taken into account as described below. 

For r^eUei^T decays we refer to the cut against residual Bhabha events 
(cut E7) which rejects electron candidates when the other side has a single 
track with x > .75, and its acoplanarity is below 0.1°. Therefore, in evalua- 
tion of the r— >cPeZ^r identification efficiency, the events where the other tau 
decay into e, /j or 7r(K) were divided to those with the x of the opposite side 
below and above 0.75. We did not find any statistically significant difference 
between the efficiency distributions of those event subsamples. In the case 
where the other r do not decay into these three channels our criteria are 
independent of its x value. 

For the t^^v^Vj. decay, the against yU-pair events (cut M5) introduces a 
dependence on the opposite r properties. This cut rejects muon candidates 
if the opposite jet is consistent with being a muon and has x > 0.8. Hence, 
for events where both taus are identified as muons we require both x-values 
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to be below 0.8. This condition is applied already in the parametrization of 
£^{xi, Xj,cos9) and therefore does not have to be repeated. However there 
might be still some residual effects of this cut on events where the opposite jet 
contains a charged pion, which resemble a muon. Therefore, we investigated 
the T—>-fiUfj_i'r identification efficiency for the following cases. 

• events where the other tau does not decay into e, fi or vr(K); 

• events where the other tau decays into a lepton (e, /i); 

• events where the other side is a r-^7r(K)z/^ decay with x < 0.8; 

• events where the other side is a r— s>7r(K)z/^ decay with x > 0.8; 

It turned out that for each case the identification efficiency was different. 
Moreover, for the third case, there was even a difference between positive and 
negative r~ helicity events. This is due to the fact that the pions in positive 
T-^fiijfj_i>r helicity events are more energetic and thus are more sensitive to 
the cut. We had, therefore, to parametrize each case separately. 

The efficiency plots for T^eUe'^r, '^^fJ'i^tJ.'^T and r— s>7r(K)z/,- are given in 
Figs. |6.11[ - |6.13| along with their parametrization curves. To allow a possible 
correlation between the efficiency dependencies on x and cos 6' we use in 
the analysis a two-dimensional distribution and fit it to an expression which 
depends on both. The list of 'x^ / D.O.F. values of all the efficiency fits are 
given in Table |6.12| . 

Control samples 

A shortcoming of efficiency evaluation based on MC distributions alone 
comes mainly from possible effects not modelled correctly in the detector 
simulation. Therefore, events from control samples of the data are used 
to correct the MC derived r decay identification efficiency. The selection 
procedures of these control samples are largely independent of the r decay 
identification requirements. The efficiencies obtained from those real data 




Figure 6.11: Efficiencies of T^eug^r (a,b) and r— >7r(K)i/,- (c,d) identification 
as functions of x and | cos^|. 




Figure 6.12: Efficiencies of T^jjiu^Ur identification as function of x and 
I cos6'|. In (a,b) the decay mode of the second hemisphere is not identified, 
where in (c,d) the second r decays into a lepton. 
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Figure 6.13: Efficiencies of r— identification as function of a; and | cos 6*1 
for different cases (see also text) . 
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Decay channel combination 


X^/D.O.F. 


r— s^ez/g'^T 
r^7r(K)z/T- 

T^fiu^Ur, none 

r^fiu^Ur, r^7r(K)z/T-, x < 0.8, pos. helicity 
r^fiu^Ur, r^7r(K)z/^, x < 0.8, neg. helicity 
r^yuz/^z/^, r— >7r(K)z/T-, x > 0.8 


17.7/27 
26.5/28 
41.0/292 
14.5/24 
19.9/30 
28.6/28 
25.5/28 



Table 6.12: jD.O.F. of the parametrization fits for the r decay identifica- 
tion. 



events are compared with those from MC simulating the same process and 
the same selection criteria. In this way, any bias due to the selection of the 
control sample itself is minimized. 

For r— ^^ez/ei^r decay, Bhabha and single electron events are used. Fig. |6.14| a 
shows the ratio between the efficiencies of the data and MC control samples, 
where the low energy points are from the single electron events and the 
high energy point is from the Bhabha events. The straight line fitted to 
these points parametrizes the correction which must be multiplied with the 
r— >-ez/eZ/r MC derived efficiency. 

For Bhabha events, namely for high energy electrons, the efficiency cal- 
culated from the data is lower by 4.1% than the efficiency calculated with 
the MC Bhabha events. For single electron events, which are low energy 
electrons, the data efficiency is lower by 5% than the MC derived efficiency, 
leading to an energy dependent correction of (0.949±0.005) + (0.009±0.005)a;. 

In order to investigate the muon detection efficiency we used /i-pair events 
for high energy and 77— ^/i"*"/!" events for low energy muons. The result- 



ing correction for effects not modelled in the MC, shown in Fig. |6.14j b, is 
(0.961 ± 0.010) + (0.026 ± 0.010)x. The uncertainty of these corrections, for 
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Figure 6.14: Corrections to the MC derived efficiencies for the tau identifica- 
tion using data control samples: a) for r^ez/e'^r decay, using single electron 
and Bhabha events; b) for r^iiu^u^. decay, using e"'"e~—>-e"'"e~/.i and /x- 
pairs; c) for T^7r(K)i/,- decay, using r^pur events (see text). The sohd fines 
are our parametrization. The dotted lines in b) are alternative parametriza- 
tion used to estimate the systematic uncertainties. 
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both electrons and muons, is increased beyond the quoted errors of the fit 
parameters in order to allow for possible non-linear interpolations between 
the extreme x regimes of the control samples. 

For single charged pions or kaons there are no clean non-tau control sam- 
ples available over the momentum range of interest. Instead, decays of the 
tau into hadronic states containing neutral pions and a single charged hadron 
are used. For example, the 7r(K) identification requirements that remove 
muons (see P4) can be controlled using the r—^pur sample since these re- 
quirements are not used in the r-^pVj. identification. The resulting correc- 
tion for these requirements is approximately uniformed equal to 0.98, but 
can also be parametrized as a second order polynomial. This ratio between 
data and MC is shown in Fig. |6.14| c. We use the flat correction function, but 
investigate the effect of using the other alternative for the systematic study. 
We investigated this correction also separately for positive and negative cos 9 
values and did not find any statistically significant difference. 

In order to study the requirements designed to remove background con- 
taining neutral pions (see P3) an alternative pion selection based on a low 
jet mass requirement is used. This sample has larger background (12%), but 
is clean enough to estimate the systematic dependence of the results on the 
requirements on the presampler signal and unassociated ECAL clusters. The 
conclusion of this study is that the MC modelling of the x dependence of the 
efficiency is in excellent agreement with the data. This is verified using a sub- 
sample of the T^pi>T- events in which the neutral pion is well separated from 
the track and with a sample of very high momenta r^'7r(K)z/^ events which 
had little background from the other hadronic decays. The statistical error 
obtained from this control sample investigation is taken to be a systematic 
error. The r— >pz/^ events are also used to demonstrate that the requirement 
that removes electrons (see P2), introduces a negligible contribution to the 
overall systematic error. 

To conclude, three different sources of systematic errors associated with 
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the decay identification efficiencies were considered. 

• The MC statistics, i.e. the statistical errors of the parameters obtained 
by the ffis to the MC efficiency distribution. 

• Control sample statistics, i.e. the statistical errors is the parameters 
describing the corrections to the MC efficiencies obtained from the con- 
trol samples. For the t^tt{K)i>t- channel we made an alternative fit to 
a non-flat line, a + bx, and we used the deviation of the b parameter 
from and its statistical error. 

• Uncertainty in the interpolation of the corrections obtained from the 
control samples. This uncertainty refers only to the r— ^ez/gZ^T and 
t^^v^Vt- channels where the control samples cover only very low and 
very high lepton energies. In order to estimate the effect of this un- 
certainty, we applied two alternative interpolations, in addition to the 
straight line between the two regions. These interpolations are de- 
scribed by the dotted lines in Figs. |6.14|a,b. 



The effect of each uncertainty was propagated to the final {Pr) and A^^f 
results using the same method as used in the previous sections, and the results 



are summarized in Table |6.13| . For each decay channel we list separately the 
effect on the {Pr) and on the A^^f results of the whole sample, and on the 
results of a subsample containing the individual decay modes. 

These various efficiency corrections contribute systematic errors on {P^) 
and A^„f of 1.79% and 0.19%, respectively. 

6.6 Background 



As seen from Eqs. ^]2| and |6.5| , the correction for background contamination in 
channel i is performed separately for background from other r decay sources 
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whole sample 


corr. subsample 


source 


T 


decay 


A{Pr) 








MC statistics 


T- 


-^ez/gZ/,- 


0.52 


0.14 


2.36 


0.64 




T- 




0.20 


0.08 


0.85 


0.34 




T — 


^7r(K)z/^ 
all 


0.77 
0.95 


0.07 
0.18 


1.42 


0.13 


Control sample 


T- 


— i>eZ/eZ/T- 


0.15 


0.00 


0.68 


0.01 




T- 




0.24 


0.01 


1.02 


0.04 




T — 


^7r(K)z/^ 
all 


1.33 
1.51 


0.06 
0.06 


2.47 


0.10 


Interpolation 


T- 


-^ez/gZ/,- 


0.04 


0.00 


0.18 


0.01 




T- 


all 


0.11 
0.11 


0.00 
0.01 


0.61 


0.04 


All sources 


T- 


-^ez/gZ/r 


0.52 


0.14 


2.46 


0.64 




T- 




0.38 


0.08 


1.46 


0.34 




T— 




1.54 


0.08 


2.85 


0.16 




Total 


1.79 


0.19 





Table 6.13: Systematic errors (in %) associated with tau decay identification 
efficiencies. 

(Pf) and background from non-r sources {P]jj)~^)- There are two reasons 
for this separation of background sources. 

• The background from other r decay channels depends on the r helic- 
ity of the event and therefore it must be calculated for each helicity 
separately. 

• The background from non-r events depends on the whole event, whereas 
the background from other r decays depends only on the r-jet consid- 
ered. 



For these reasons we had to handle those two background types differently 
(see Eqs. |6^ , |6.5| ). In Section |6.6.1| we describe only the first background 
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and in the next Section we shall present the second one. 

6.6.1 Background from Other Tau Decay Channels 

The background from other tau decays was initially investigated using our 
r-pair MC sample and separating between the two r~ helicity states and be- 
tween five decay channels, r^eUgi^T, T^f^i^fi'^T, t^tt Ur, r— »/i7r°z/ and com- 
bining all other channels together. The separation between the sources was 
done in order to be able to correct for the MC branching fractions of the dif- 
ferent decay channels and to investigate systematics (see below). We consid- 
ered separately the dependence of the background on the kinematic variables, 
Pf{xi,cos9), and its overall normalization, (3f = J dxidcos 9 /3f{xi, cos 9). 

To obtain the overall normalization, we calculate the background fraction 
to channel i from source j by, 

n*^^ r ■ 

= Jlcf^ (6-23) 

where r?-^'" is the number of MC events from source j, identified as coming 
from channel i. The correction factor q is the ratio between the branching 
fractions to channel i in data and MC. For the branching fractions in the 
data of the T^ev^^ri '^^^fJ'^fi'^r T^7r(K)z/^ and T^hTc^u channels we are 
using the published OPAL results 0, No correction was applied to the 
background from the fifth source (all other channels). Table |6.14| lists the 
MC and data branching fractions. 

The overall background fraction for channel i is obtained by summing up 
the contributions from the five sources, 

Pt = E5=i A^- Table |A5| lists the 



various contributions to the background fractions and their sums. 

The dependence on the kinematic variables of the background fractions 
was also taken from the MC events. The x and cos 6' MC background dis- 
tributions for each helicity and decay channel, and from each source were 
fitted to third order polynomials. The normalization of each background 
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Figure 6.15: Background from other tau decays as functions of x to the 
channel T—^eu^i^r (Figs. a,b), T^^v^Vr (Figs. c,d) and r^7r(K)z/^ (Figs. 
e,f). Figs. a,c,c (b,d,f) correspond to positive (negative) hehcity events. In 
Figs. a,b contamination contributions from r— >7r(K)t/^, r— ^pz/^- and other 
T decay channels (including t^ixV^Vt) are presented by open, hatched and 
cross-hatched histograms correspondingly. In Figs. c,d contamination con- 
tributions from r— *7r(K)z/T- and all other r-dccay channels are shown by open 
and hatched histograms correspondingly. In Figs. e,f contamination contri- 
butions from T^p^T, T^Qi'e^r-i 'T^l^^n^r and other r decays are drawn with 
open, hatched, black and cross-hatched histograms correspondingly. The 
curves are our parametrization. The vertical scale is in arbitrary units. 
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Figure 6.16: Background from other tau decays as functions of | cos^^|. to 
the channel r^ez/e^^r (Figs. a,b), r— >yuz/^z/^ (Figs. c,d) and r^7r(K)z/^ (Figs. 
e,f). Figs. a,c,e (b,d,f) correspond to positive (negative) hehcity events. In 
Figs. a,b contamination contributions from r— >7r(K)i/^, r— ^pz/^- and other 
T decay channels (including t^hP^Vt) are presented by open, hatched and 
cross-hatched histograms correspondingly. In Figs. c,d contamination con- 
tributions from r— *7r(K)z/T- and all other r-decay channels are shown by open 
and hatched histograms correspondingly. In Figs. e,f contamination contri- 
butions from T^p^T, T^Qi'e^r-i 'T^l^^n^r and other r decays are drawn with 
open, hatched, black and cross-hatched histograms correspondingly. The 
curves are our parametrization. The vertical scale is in arbitrary units. 
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r— s>7r(K)z/^ 






18.25 


17.76 


12.59 


24.77 


J^data 


17.62±0.47 


16.80±0.41 


12.22±0.59 


26.27±0.61 


^ j^data 1 j^MC 


0.965±0.026 


0.946±0.023 


0.971±0.047 


1.061±0.025 



Table 6.14: Tau branching fractions (in %) used in our analysis. 



r~-liel. 


positive 


negative 


r decay 


r^ez/gZ/,- 




r— >7r z^T- 


r— ^efgi^T 






r— ^ez/eZ^r 




0.00 


0.92 




0.01 


0.92 




0.03 




0.26 


0.02 




0.27 


r^7r(K)z/T- 


2.34 


1.01 




2.68 


1.76 






1.52 


0.03 


5.55 


1.83 


0.03 


3.68 


other 


0.19 


0.02 


1.57 


0.17 


0.01 


1.96 


overall 


4.08 


1.06 


8.30 


4.70 


1.81 


6.83 



Table 6.15: Background fractions (in %) to the three decay channels 
T—>-ei'e^r, 'T■^^J'V^l'^-, T— ^>7r(K)z/^ (columus) froui the various sources (rows). 



source was then corrected by multiplying with the corresponding branching 
fraction correction c. The contributions from all sources were then summed 
up and normalized to j3f. The resulting functions, l3f{xi,cos6), are plotted 
in Figs. |6.15| and |6.16| , together with the corrected MC distributions. The 
values for the agreement between the MC distributions and their polynomial 
parametrization are indicated on the plots. 

For T^eUei^T, the background level from other r decays is ~ 4.7%, differ- 
ing slightly between positive and negative r~ helicity events, and originating 
mainly from r^7r(K)z/^ and r^pUr decays. For T^pv^Vr this background 
is ~ 1.5% coming from r— s>7r(K)z/T- decays, and the background to r— s>7r(K)z/T- 
is ~ 7.8%, mainly from r—^pUr, t^ttK^u and r-^eVe^r decays. The MC 
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expectations for these background levels are checked against the data by 
studying the tau jets that are identified in more than one decay channel and 
from studies of contaminations with Bhabha and muon-pair control samples. 

The following contributions to the systematic error associated with the 
background from tau sources were considered. 

• The uncertainties in the parameters of the background functions, origi- 
nating from the statistics of our MC sample. An error on the t^ttK^u 
background is estimated by varying the MC expectation by ±100%. 

• A small (<0.1%) additional contribution to the errors arising from the 
experimental uncertainties in the tau branching ratios has been in- 
cluded to correct the MC numbers. 

• Possible disagreement between the data and MC simulation. This was 
investigated in two ways. 

1. We used our Bhabha and /x-pair control samples and counted the 
number of particles identified as pions. We compared these num- 
bers for data and MC, and the ratios were 0.72±0.48 for electrons 
and 1.17±0.39 for pions. These ratios represent correction fac- 
tors to be multiplied with the MC estimates of the background to 
T^7r(K)z/T- from T^eUeUr and T^ixv^Vr decays. 

2. We counted tau decays which were identified by two decay algo- 
rithms and compare it to MC. Almost all those identified as a 
pion and electron or as a pion and muon are pions (from MC). 
Therefore we can get a correction of the pions background to the 
electrons and muons. These corrections turn out to be 0.41±0.09 
and 0.58±0.08 respectively. 50% of those identified as a rho and 
pion or electron or muon are p's. Therefore we can get a correction 
to the background from rho to electron, muon and pion which are 
1.30±0.10, 5.0±1.0 and 1.09±0.09. 
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These corrections are applied in the systematic studies. We use either 
their uncertainties or their difference from the factor one (whichever 
greater) and see the corresponding effects on the final results. 



Table |6.16| summarizes the systematic errors associated with the back- 



ground from other tau decays. The contribution of the branching fraction 
uncertainties actually affect not only the background estimates but are used 
also in the expression for cos6') (Eq. and in the estimate of the 

tt/K mass threshold effect. 



source 


all 


r^ez/gZ/^ 
A(P.) AA^^f 


r->/i z/^z/^ 
A(P.) AA,^i 


MPr) AA^„f 


MC 
Data 


0.29 0.04 
0.58 0.09 


0.50 0.09 
1.49 0.06 


0.44 0.06 
0.25 0.16 


0.45 0.05 
0.10 0.01 


Diff 


0.32 0.02 


0.15 0.08 


0.22 0.01 


0.49 0.02 


Total 


0.72 0.10 


1.58 0.13 


0.55 0.17 


0.67 0.05 



Table 6.16: Systematic errors (in %) associated with background from other 
tau decays and tau branching fractions. MC, Data, and BR stand for MC 
statistics, the difference between data and MC and Branching ratio fraction 
respectively. 



These background sources contribute errors of 0.72% and 0.10% on (P^) and 
Apof > respectively. 

6.6.2 Background from non-tau Events 

Concerning the residual non-r background after the selection criteria, the 
following sources are considered. 

• e'^e~—>-fi~^fi~. This source mainly contaminates events where one r is 
identified as r—^^v^Vr decay and the other as t^ti{K)vt- or when one 
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tau is identified as a T^jiu^Ur and the other is not identified, at the 
level of 0.4% and 1.3%, respectively. 

• e"'"e~— >-e"'"e~/i~^yU^. This source contributes only to events with two 
identified muons or one identified muon and an unidentified tau, at the 
level of 1.6% and 0.3% respectively. 

• e+e"— ^-e+e". This source contributes background to events where only 
one r decay is identified as e or 7r(K) and the other tau decay is uniden- 
tified, and events where both decays are identified as r— s>ei>eZ^r- The 
corresponding contamination levels are 0.3%, 0.2% and 0.3%, respec- 
tively. 

• e"'"e~— >-e"'"e~e"'"e~. This source contaminates only events with one elec- 
tron where the opposite tau decay is unidentified or two identified elec- 
trons, at the level of 0.3% and 2.6%, respectively. 

• e"'"e~— The background from this source is negligible for all channels 
considered here. 

The j3'^"^~^ {xi,Xj, cos 6) correction functions for each source were initially 
determined using the corresponding MC event samples. In practice the cor- 
relations between and cos 6 are small and the /^j^j""^ functions factorize 
into products of simple functions. These were adjusted by factors obtained 
from comparisons between the data and MC when the tau selection require- 
ments designed to suppress these backgrounds were loosened. 

Table |6.17| lists the five sources which were considered where for each 
source the size of the MC sample is quoted. These numbers were estimated 
from the integrated luminosity collected by OPAL during 1990-1992. In 
addition, we list the expected background fraction for each combination of 
decay channels. 



CHAPTER 6. PHYSICS ANALYSIS 



139 



source 


ee— 


ee^qq 


ee— i>ee 


ee— s>ee/i/i 


ee^eeee 


total 


MC 


200000 


702400 


46187* 


50.23 pb-i 


42.98 pb^i 




e, none 


0.05±0.02 


0.00 


0.31±0.12 


0.00 


0.25±0.08 


0.61±0.15 


/i, none 


1.29±0.09 


0.00 


0.00 


0.26±0.08 


0.00 


1.55±0.12 


IT, none 


0.00 


0.00 


0.18±0.13 


0.00 


0.00 


0.18±0.13 


e, e 


0.00 


0.00 


0.28±0.28 


0.00 


2.60±0.65 


2.88±0.71 


e, /i 


0.14±0.05 


0.00 


0.00 


0.00 


0.00 


0.14±0.05 


/i, /i 


0.05±0.02 


0.00 


0.00 


1.60±0.40 


0.00 


1.65±0.40 


/i, 71 


0.40±0.16 


0.00 


0.00 


0.00 


0.00 


0.40±0.16 



* numbers correspond to | cos6\ < 0.9 

Table 6.17: Non-tau background sources and their relative contributions (in 
%) to the various combinations of r decay channels. The combinations not 
mentioned in this table have no background. 



The X and cos^ dependence was taken from the MC distributions and 
parametrized by simple functions. Fig. |6.17| shows these distributions for 



those decay channel combinations where the overall non-tau background is 
above 0.5%. 

Two sources of systematic errors associated with the non-tau background 
were considered. 

• Finite statistics of our MC event samples, illustrated by the given errors 
in Table |6.17|. To estimate their effect, we varied both the background 



levels and the parametrization of their x ans cos 6' dependence within 
the statistical errors and examined the effect on the final results. 

Possible incompatibility between MC and data. This was investigated 
by J. Clayton |^ using background enriched subsample where he ob- 



tained correction factors for all the non-tau background sources. As 
we did for the tau background, we apply these correction factors only 
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Figure 6.17: Background from non-tau sources to events with one identified 

T^ePeT^T s-nd the other side not identified (a,b); events with two identi- 
fied T^eVf.Vj. (c,d); events with one identified t^^v^Vt and the other side 
not identified (e,f); events with two identified r-^fiu^Ur (g,h). Figs. a,c,e,g 
(b,d,f,h) show the x (cos 6*) dependence. In Figs, a-d, the open, hatched and 
cross-hatched histograms represent the background from Bhabha, ee— >eeee 
and //-pairs respectively. In Figs, e-h, the open and hatched histograms are 
the background from yU-pairs and ee^eefi^ respectively. The dashed lines 
are our parametrization. The vertical scale is in arbitrary units. 
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for the systematic study. We use either their uncertainties or their 
difference from the value one (whichever greater) and observe the cor- 
responding effects on the final results. 



Table |6.18| summarized the systematic errors associated with the non-tau 



background. 



source 


all 


r^ez/eZ/r 
A(P.) AA,^i 




A(P.) AA^i 


MC 
DiflF 


0.48 0.16 
0.48 0.08 


1.56 0.60 
0.72 0.26 


1.20 0.17 
1.97 0.26 


0.33 0.15 
0.07 0.03 


Total 


0.68 0.18 


1.72 0.65 


2.31 0.31 


0.34 0.15 



Table 6.18: Systematic errors (in %) associated with the non-tau background. 
MC and Data stand for the MC statistics and the difference between data 
and MC. 



The contribution to the errors on {Pr) and A^^f from the non-tau back- 
ground corrections are 0.68% and 0.18%, respectively. 



Chapter 7 

Results and Cross Checks 



7.1 Maximum Likelihood Fit Results 

The ML fit, applied to the T^eUe^^T, T^fJ'^^ii^r and r^7r(K)z/^ decay channels 
yields the following results, 

{Pr) = (-13.5±2.9(stat))% 

APB = {-11.0 ± 3. 5{stat))% (7.1) 
where the correlation between the two numbers is 0.03. 



The r W z^r Coupling 

Since the r decays via weak interaction, where parity is not conserved, the 
angular distribution of its decay products in the r rest frame depends strongly 
both on the r spin orientation and the characteristics of its charged weak 
decay process. The results given in the previous paragraph were calculated 
using the theoretical Eq. [4.37| which assume a V — A structure for the r 



decay. However, following Sect. [4.3| one can use the more general expression 
presented in Eq. [4.48| , and extract not only the r production parameters but 
also its decay parameters. Our results of a six parameters ML fit to the two 
polarization asymmetries and the Michel parameters p, 6, ^ and C,had are: 
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(P^) = (-11.8±3.2(stat))% 
A^^f = (-11.8±3.7(stat))% 

^ = 0.83±0.18(stat) (7.2) 
Chad = 1.13±0.16(stat) 

p = 0.77 ±0.01{stat) 

S = 0.87±0.13(stat) 

(7.3) 

which are consistent with a {V — A){V — A) structure (see Table [4.2|) . 

The existing measurements of the Michel parameters are already precise 
enough to exclude interaction different from V — A in the iui vertex (£=e,/i), 
therefore one can use Eqs. [4. 50] to evaluate the helicity of the Uj. from its 
relations with all last four parameters. This was further checked by fixing 
the Wi vertex structure and allow only one free parameter to determine the 
characteristic of the rur vertex. The fit results are: 



(-12.3±2.8(stat))% 
(-11.4±3.4(stat))% 

-1.06±0.02(stat) (7.4) 

(7.5) 

The measurement of the chirality parameter and the leptonic Michel param- 
eters are all consistent with a h^^ = —1 which is expected from a pure V — A 
coupling at the Ti/f vertex. It is also an indication for the vr being a left 
handed particle. 



{Pr) 

\FB 
^pol 
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7.2 Summary of Systematic Studies 



Table [7.1| summarizes the systematic errors in the ML analysis. Each error 



listed in Table is a combined result of several related contributions. Since 
are no correlations between entries on different rows of the table, they are 
combined in quadrature to give the overall error listed in the last row. Thus, 
the values of the systematic errors of (Pr) and A^^f using the r— ^-ez/gZ^T, 
T^liv^Vr and r^Ti{K)Vr decays are measured to be 

A(P,) = ±2.2%{sys.) 

AAj,f = ±Q.h%{sys.). (7.6) 



Source 


(%) 


A AS 

(%) 


r decay identification 


1.79 


0.19 


backg. from other r decays 


0.72 


0.10 


radiative effects 


0.29 


0.03 


backg. from non-r events 


0.68 


0.18 


r-pair selection efficiency 


0.39 


0.20 


calorimeter response 


0.48 


0.12 


tracking response 


0.38 


0.21 


ApB and K mass threshold 


0.10 


0.07 


cos 9 and charge measurement 


0.00 


0.17 


Total 


2.19 


0.43 



Table 7.1: Summary of systematic uncertainties in the determination of 
{Pr) and Apgf using r— ^ez/eZ/,-, r^fiu^Ur and r^7r(K)z/^ channels . 



The largest contribution to the systematic error of (Pr) is derived from 
the uncertainties in the efficiency of the r decay mode identification. For the 
most part these arise from the limited statistics of the control samples, and 
in particular, the T^pVj. sample used to correct the efficiency of r— i>7r(K)z/^ 



CHAPTER 7. RESULTS AND CROSS CHECKS 



145 



identification. Therefore, one can expect tliat tlie systematic errors will de- 
crease as more data is collected. 

7.3 Monte Carlo checks 

7.3.1 MC test of polarization biases 

In order to check for potential biases in the fitting method, the analysis was 
performed on the r-pair MC sample. The ML fit is performed on the MC 
events, removing the corrections to the r identification efficiencies derived 
from the data, and using resolution parameters extracted exclusively from 
MC events. Non-tau background correction is also switched off, since our 
sample contains only r-pair events. The input values for (P^) and A^^f are 
— 14.0% and —10.5%, respectively and the fit returns values of (— 13.9±1.0)% 
and (-8.9±1.3)%. 

The fit results are listed in Table which is similar to Table 17731 for 



the data. Although these values are consistent with the input values, we 
have changed the input parameters from purely negative to purely positive 
r~ helicity and studied the dependence of the measured asymmetries on the 



input ones. Fig. |7.1| illustrates this investigation for the three individual 
channels as well as for the whole sample. We find no evidence for a bias in 
our measurement method. 

7.3.2 Check of the various parametrization using MC 

For an overall check of our parametrization of the various corrections ob- 
tained from MC, we plot the MC x distributions separately for each decay 
combination and compare it with the corresponding corrected theoretical 
curves (Figs. |7]^, [f73| , [f^ , |7.5|) . The values indicated on the figures show 



that the agreement between the MC distributions and the theoretical curves 
is good. 
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Figure 7.1: Test of Bias in the {Pr) measurement. The error in the measure- 
ment of {Pr) as a function of the MC input polarization. Fig. a) - c) are 
for the T— >ez/e^^r) t^I^^h^t and r— >7r(K)z/^ respectively and d) present the 
results of all three channels together. 
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Figure 7.2: MC x distributions for the r^eUei^T decay channels and for 
various decays on the opposite side, compared with the theoretical predictions 
(solid lines) with all corrections included, taking the (Pr) value obtained 
from the global ML fit. The filled area represent the background part of the 
spectrum. The dotted lines are the uncorrected theoretical curves. 
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Figure 7.3: MC x distributions of the T^jdv^Vr decay channels and for var- 
ious decays on the opposite side, compared with the theoretical predictions 
(solid lines) with all corrections included, taking the {Pr) value obtained 
from the global ML fit. The filled area represent the background part of the 
spectrum. The dotted lines are the uncorrected theoretical curves. 
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Figure 7.4: MC x distributions for the r^7r(K)z/^ decay channels and for 
various decays on the opposite side, compared with the theoretical predictions 
(solid lines) with all corrections included, taking the (Pr) value obtained 
from the global ML fit. The filled area represent the background part of the 
spectrum. The dotted lines are the uncorrected theoretical curves. 
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Figure 7.5: MC x distributions for the r—^euei^r, T^fii^^i^r and T—^TT{K)h'r de- 
cays channels, summing over all decays on the other side and compared with 
the theoretical predictions (solid lines) with all corrections included, taking 
the (Pr) value obtained from the global ML fit. The filled area represent 
the background part of the spectrum. The dotted lines are the uncorrected 
theoretical curves. 
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(Pr) 


\FB 
^pol 


Number of Events 


Decay Channels 


(%) 


(%) 




All decays 


-13.9 ± 1.0 


-8.9 ± 1.3 


93462 


1 r identified 


-13.9 ± 1.3 


-8.5 ± 1.6 


67199 


(a) T^eUgi^T, none 


-11.9 ± 3.2 


-5.5 ± 3.6 


27058 


(b) T^fiu^UT, none 


-15.6 ± 3.2 


-8.1 ± 3.6 


27120 


(c) r— >7r(K)z/T, none 


-13.9 ± 1.7 


-9.6 ± 2.0 


13021 


2 r's identified 


-13.9 ± 1.7 


-9.5 ± 2.1 


26263 


(d) T^ePeJ^T, T^ePe^T 


-8.9 ± 5.7 


-15.0 ± 6.6 


4684 


(e) r—>-ePe^T, T^iiV^Vr 


-13.7 ± 3.6 


-8.3 ± 4.2 


10328 


(f) r^ez/gZ/^, T^Ti{K)vr 


-16.7 ± 3.2 


-8.2 ± 4.0 


3183 


(g) T^flU^U^, T-^flV^Vr 


-19.3 ± 9.4 


0. ± 11. 


4289 


(h) T^nu^Ur, r^7r(K)z/^ 


-13.4 ± 3.3 


-9.1 ± 4.0 


3195 


(i) T— >7r(K)i/^, T— >7r(K)i/,- 


-9.4 ± 6.3 


-19.5 ± 8.2 


584 


Average 


-13.9 ± 1.0 


-8.9 ± 1.3 




probability (%) 


3.2/8 


4.4/8 




92.4 


81.8 





Table 7.2: r polarization results (in %) from our MC sample. The average 
values the / D .0 .F. and the probabilities correspond to the combined 
results of the 9 different channels (a) to (i). The input values are {Pr) — 
-14.0% and AH = -10.5%. 



7.4 Cross Checks 

7.4.1 Statistical error examination 

In order to convince ourselves that the statistical errors obtained from the 
fit are correct, we used a large sample of generated events without detector 
simulation and perform a simple ML fit to each decay, ignoring the other side 
and ommiting the corrections. We only applied the x > 0.05 and | cos^| < 
0.68 requirements. The errors obtained, after scaling them for the different 
numbers of events, were consistent with the statistical error obtained in our 
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global fit within 3%. 

7.4.2 Separate checks with the different decay topolo- 
gies 

The data has also been classified into nine independent subsamples, corre- 
sponding to all possible combinations of both r decays excluding an identified 
T^pi^T- decay. Doing this, nine independent results for {P^.) and A^^f are ob- 
tained. These are listed in Table |7.3| . The weighted mean of these values 
are within 0.3% of the global fit results and the probabilities for the con- 
sistency between the nine results for (Pr) and A^^f are 17.2% and 81.6%, 
respectively. 

7.4.3 Graphical evaluation of the fit quality 

In order to evaluate the fit quality, we plot the x distributions separately 
for each decay combination and compare it with the corresponding corrected 
theoretical curves (Figs. [^6] , [777| , [7.8| ). The bins corresponding to a; > 1 
are not used in the fit, since they have low statistics and are dominated 
by the non-tau background which has large uncertainties. We note that 
the complex behavior of the theoretical distributions in the low x region 
where both taus decay to electrons or both decay to muons is the result of 
the need to remove two-photon events. The effect of the kaon threshold in 
the r^7r(K)z/T- decays is also evident on these plots. The x distributions for 
T^eUeT^T, T^f^i^fi^T and T-^niK)^^-, summing over all decays of the opposite 
tau are presented in Fig. \I.9\ The normalization of the theoretical curves 
is obtained by a fit to the experimental distributions. The low values 
verify the good agreement between the experimental distributions and the 
theoretical curves. The filled regions in the figures represent the background 
and the dotted lines are the uncorrected theoretical distributions. One can 
see that the background contributions are small, and the overall corrections 
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Decay Channel combination 


(%) 


(%) 


Number of Events 


Gl()l)al fit value 


-13. 5± 2.9 


-11.0± 3.5 


12397 


1 r identified 


-17.5± 3.5 


-10.2± 4.3 


9250 


(a) T^eiJeUr, none 


-24.7± 8.9 


-17.±10. 


3585 


(b) T-^fiUfj^Ur, none 


-8.7± 7.9 


-10.8± 9.0 


3897 


(c) r— >7r(K)i/T, none 


-18.4± 4.4 


-8.1± 5.5 


1768 


2 r's identified 


-5.5± 5.0 


-12.6± 5.9 


3123 


(d) T— >ePei^T, T^ePe^T 


9. ±16. 


-28.±18. 


565 




-16.±10. 


3. ±12. 


1288 


(f) T^eUel^r, T^7l(K)Ur 


-9.5± 9.7 


-23. ±12. 


346 


(g) T^ixv^Vr, r^ixv^Vr 

(h) t^iiu^Ut, T^7r(K)i/^ 


7. ±24. 
6.1 ± 8.9 


2. ±27. 
-13.±11. 


484 
385 


(i) T^7r(K)i/T-, T— >7r(K)^'T- 


-29.±20. 


-21.±26. 


55 


A\'erage 


-13. 2± 2.9 


-11. 1± 3.1 




xyD.o.F. 

probability (%) 


11.6/8 


4.4/8 




17.2 


81.6 





Table 7.3: r polarization results for the nine (a to i) independent subsamples, 
corresponding to all possible combinations of both r decays to e, /i or 7r(K). 
The average values, the / D.O.F. and the x'^ probability correspond to the 
combined results of these 9 different r decay channels. 



to the theoretical curves, excluding the two lowest x values, are smooth, and 
do not depend strongly on x. 

7.4.4 A fit consistency check 

Another check of the ML fitting method is a fit of the x distributions 
to a linear combination of the corrected theoretical curves for positive and 
negative helicity. The x"^ fit is meaningful only for the case where there is 
enough entries (at least 10) in each of the bins. Therefore, it was only applied 
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Figure 7.6: x distributions for the r—^eUeJ^T decay channels and for various 
decays on the opposite side, compared with the theoretical predictions (solid 
lines) with all corrections included, taking the (Pr) value obtained from the 
global ML fit. The filled area represent the background part of the spectrum. 
The dotted lines are the uncorrected theoretical curves. 
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OPAL 




Figure 7.7: x distributions for the r^^iv^Vr decay channels and for various 
decays on the opposite side, compared with the theoretical predictions (solid 
lines) with all corrections included, taking the (P^) value obtained from the 
global ML fit. The filled area represent the background part of the spectrum. 
The dotted lines are the uncorrected theoretical curves. 
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Figure 7.8: x distributions for the r— >7r(K)z/T- decay channels and for various 
decays on the opposite side, compared with the theoretical predictions (solid 
lines) with all corrections included, taking the (Pr) value obtained from the 
global ML fit. The filled area represent the background part of the spectrum. 
The dotted lines are the uncorrected theoretical curves. 
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Figure 7.9: x distributions for the r— ^ez/gZ^T, T^^v^Vr and T^'K{K)vr decay 
channels, summing over all decays on the other side and compared with 
the theoretical predictions (solid lines) with all corrections included, taking 
the (Pr) value obtained from the global ML fit. The filled area represent 
the background part of the spectrum. The dotted lines are the uncorrected 
theoretical curves. 
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to the case where just one tau was identified. The results which are hsted in 
Table [7.41 are consistent with the ML fits listed in Table 1773. 



Decay Channel 


(Pr) 


T—>-eUel^T 
r—>-TT{K)l'r 


-21.1 ± 9.1 
-4.2 ± 8.1 
-20.1 ± 4.6 



Table 7.4: Tau polarization results calculated by fits to events where just 
one tau was identified. 



7.4.5 Measurements of polarization as a function of 

cos^ 

The results were also checked by dividing the cos 6 range into five intervals to 
determine the tau polarization in each of them separately using the ML fit. 
The results are plotted in Fig. [7.10| . The solid line in the figure describes re- 
sult of a fit of the tau polarization expected dependence on cos 9 (Eq. |4.27|) . 
The OPAL result for the Z° peak Afb = 1-2% @ is used (see Table 
We obtained ((P^) = -13.9±2.9(stal))% and Aj„f = (-10.8±3.5(stat.))%, 
which is in excellent agreement with our nominal results. The fitted is 0.7 
for 3 degrees of freedom and the correlation between {Pr) and A^^f is 0.027. 

7.4.6 A test of the Aps influence 

We also modified our fitting procedure allowing the seven Aps values to 
be free parameters in the fit. They are consistent with the values which 
have been used in our standard fit, where the largest deviation is 1.8a. Our 
statistical errors for Aps are twice as large since we use only events in the 
barrel region and only those events where at least one tau was identified. 




Figure 7.10: Measured tau polarization results as function of cos 6'. The solid 
line is our fit result (see text) . 

The corresponding results for {Pt) and A^^f were (—13.5 ± 2.%{stat.))% 
and (—10.9 ± ?>.b{stat.))% respectively, in an excellent agreement with our 
standard global fit values. 

7.5 {Pr) and Aj^f from r-^pUr decays 

In order to combine our results with those measured from the r—^pv^ de- 
cays, we quote here the OPAL (P^) and A^^f measured from this decay 
channel. The two T^pvr analyses done by OPAL has used a least squares fit 
to determine (Pr)6» in five cos^ intervals, and than extracted the {Pr) and 
Apgf parameters in a least squares fit of the {Pt)0 dependence on cos^. 

The results of these two T^pv^ analyses were combined, taking care of 
the fact that these two data sets have an overlap of 40%. The T^pvr OPAL's 
result is then: 
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{Pr) = (-15.7 ±2.4 ±1.5)% 

AFB = (-7.1 ±2.8 ±1.5)% (7.7) 

It is worthwhile to know that there is no overlap between the events used 
in our analysis for r decay to e, /i or 7r(K) and the events used by OPAL for 
the study of r-^pv^-. Specifically, in our analysis the events of the type £-p 
(i.e. one tau decays to a lepton and the other tau decays to a p-meson) are 
discarded. On the other hand in the r^pvr analyses events of the type vr-p 
were discarded. For events of the type p-p each identified r^pur decay was 



weighted in the r^pu^ analyses by 0.5. Table |8.1| summarize the results of 
the r — >• e, p or 7r(K) channels (ML fit), the results of the T^pVj. channel 
(x^ fit) and the combined results. 

7.6 Comparison With Previous Measurements 

7.6.1 Previous determination of the Michel parame- 
ters 

Muon decay has been used for determination of the structure of the charged 



weak interaction since the 1950's |148| , pq] . The Particle Data Group |^ 



values for the Michel parameters measured in muon decays are p = 0.7518 ± 
0.0026 |7|, ^ = 1.0027±0.0079±0.0030 l?!, and 5 = 0.755±0.009 0]. The 



ARGUS collaboration has recently published a measurement of two of the 
Michel parameters evaluated in leptonic tau decays |j73|. For pr^fj. and pr-^e 



they obtain from the e and p momentum spectra the values of 0.76 ± 0.07 ± 
0.06 and 0.79 ± 0.08 ± 0.06 respectively. Assuming C,T~>fi = ^r-»e they have 
measured C, = 0.90 ± 0.15 ± 0.10. The ALEPH collaboration have used the 
correlation in r pairs where both r's decay into hadrons (vr or p) to measure 
^had- Their values are 0: = 0.95 ± 0.11 ± 0.05, = 1.03 ± 0.11 ± 0.05 



or assuming = they obtained 0.99 ± 0.07 ± 0.04. 
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T decay 


Exp. 


{Pr) (%) 


Alf (%) 

pot V ' ^/ 


T — >PU„U 


L3 


-12 7±9 7±6 2 






DELPHI 


-12±27±8 






ALEPH 


-22 5±8 5±4 5 






OPAL 1 


20±13±8 


-16±19 




OPAL 2 


-8.5 ± 5.8±4.5 


-10.4 ± 6.6±1.3 


T — V 


L3 


-2 0± 10 1±5 5 






DELPHI 


-5±18±7 






ALEPH 


-15 4±6 5±2 9 






OPAL 1 


-17+16+10 


-88±22 




OPAL 2 


-8.0 ± 5.4±3.3 


-10.1 ± 6.2±1.3 


r— ^>7r(K)z/T- 


L3 


-14.8± 9.7±6.2 






DELPHI 


-35±11±7 






ALEPH 


-13.3±3.1±1.8 






OPAL 1 


-8±10±7 


-34±16 




OPAL 2 


-14.3 ± 3.7±3.0 


-10.9 ± 4.5±0.8 


AU 


L3 


-13.2± 2.6±2.1 






DELPHI 


-24±7 






ALEPH 


-14.3±2.3 


-9±2 




OPAL 1 


-9.7±8.4 


-22±10 




OPAL 2 


-11.3 ± 3.3 


-10.6 ± 3.3 



Table 7.5: The LEP experiments r polarization published results for the 
r— >-ei>eZ/^,r— i>/iz/^z/^ and r— i>7r(K)i/^ channels. When two errors are quoted 
the first is the statistical and second is the systematic error, when one error 
is quoted it is the combined statistical and systematic error. The L3 [|1^] and 
ALEPH [|ri| published the study of 1990-1991 data using five decay channels 
of r (e, /i, 7t{K), p and ai). DELPHI'S publication ||^ is based on the analysis 
of four decay channels (e, fi, tt{K) and p) from the years 1990-1991. The 
OPAL 1 is the OPAL first publication ^ based on the analysis of three decay 
channels (e, /i and vr(i^')) during 1990. OPAL 2, is the present analysis using 
1990-1992 e, /i and tt{K) decay channels. Note that for the comparison 
with the other numbers the individual results quoted assume, incorrectly, 
that there are no correlations between the taus. These correlations are fully 
taken into account in the global ML fit result. 
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7.6.2 Results to compare with other measurements 

In the analyses published by OPAL in the past [0] and by other LEP col- 
laborations 0-[|ri|], one did not look at the whole r-pair event. Instead, one 
did consider each decay separately, ignoring the opposite side. In order to 
be able to compare our results with those of other analyses, we repeated 
our ML fit adopting their policy. For that, we had to recalculate those cor- 
rections depending on the whole event, such as r-pair selection efficiencies 
and non-tau background. Each event, where both r's were identified, was 



handled as two independent events. In Table 1?^^ we list the results for each 
decay channel separately. One should note, however, that the results for the 
T^eVei^T) T^f^'i^fi'^T and r— s>7r(K)z/^ decays in this table cannot be interpreted 
as independent polarization measurements because the same tau-pair event 
can contribute to two channels. 

7.6.3 LEP tau polarization results 



The following Table |7.6| , presents the up to date but preliminary {Pj.) and 
Apgf results of all other 3 LEP detectors presented at the La Thuile and 
Moriond Conferences in March 1994 The OPAL figures are those of 



the present analysis combined with the r^pur figures [0. In averaging the 
results we repeat the procedure discussed by the LEP tau polarization group 
in Ref. Three common effects, which affect all the tau decay mode 



equally, were taken into account: 

• A dependence of the results on Ecm due to a direct 7 and 7 - Z° 
interference contributions to the Z° intermediate state. These were 
corrected so that the results quoted correspond to Ecm = Mz- 

• Electromagnetic initial and final state radiative corrections. 

• Born level mass terms leading to helicity flip configurations. 
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As shown by the theoretical calculations in Refs. the 
combined effect of these three is to reduce the strength of the measured 
(Pj.) (Ap^f) by about 0.3% (0.2%). Therefore the experiments results (in- 
cluding ours) were corrected by adding this amount to the measured (Pr) 
and A^of . The four experiments results were combined without extracting a 
common systematic error. The x^/D.O.F. for the average is 2.2/3 and 1.7/3 
for {Pr) and A^^f respectively. The results of the four experiments and their 
average are summarized in Fig. [7.11. 



Experiment 


{Pr) 


^pol 




(%) 


(%) 


OPAL 


-15.3 ±2.3 


-9.2 ±2.4 


L3 


-17.4 ±2.1 


-9.8 ± 1.9 


DELPHI 


-15.4 ±2.3 


-6.1 ±2.5 


ALEPH 


-13.7 ± 1.4 


-9.5 ± 1.3 


LEP Average 


-15.0 ±0.9 


-9.1 ±0.9 



Table 7.6: LEP results for (P^) and A^^f . The OPAL figures are the updated 
results where the figures of all other three experiments are preliminary results 
presented at MORIOND 94. 
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Figure 7.11: LEP results for (P^) and Aj,f . The ALEPH, DELPHI and L3 
are preliminary results given at the MORIOND 94 conference, whereas the 
OPAL figures are the updated results sent for publication. The shaded areas 
are the averages of the 4 experiments. 



Chapter 8 
Summary 



We have measured the average tau polarization and its forward backward 
asymmetry using the tau decay channels r^ez/eZ^r, t—^^v^Vt and r— >7r(K)z/T-. 
The results of these measurements (taking into account all corrections) are 

(P^) = (-13.5 ±2.9 ±2.2)% 

Apof = (-11.0 ±3.5 ±0.5)% (8.1) 

where the first error is statistical and the second systematic. Our Maximum 
Likelihood fit method takes into account correlations between the decay dis- 
tributions of the T~ and the r+ in the same event due to their opposite he- 
licity. The theoretical expectation had to be corrected for radiative effects, 
detector resolution, efficiencies of r-pair selection and tau decay identifica- 
tion, and background originating from other r channels and non-r events. 
These correction were calculated from Monte Carlo events and whenever 
possible, checked with control measured data. 

A summary of the measurements from all channels is presented in Ta- 
ble pTT| . The results are quoted for a/s = mz- A very small correction (1% 
of the (Pt)), which takes into account the fact that some data was collected 
off the peak of the Z° resonance, is made. Our final result, combined with 
the OPAL T^pUr results are: 
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Decay Channel 


(Pr) 


\FB 




(%) 


(%) 


T^eUei^T, fJ'i^fi^T,'^{K)h'T present analysis 


-13.5±2.9±2.2 


-n.0±3.5±0.5 


OPAL T^pUr 


-15.7±2.4±1.5 


-7.1±2.8±1.5 


OPAL combined results 


-14.9±1.9±1.3 


-8.9±2.2±0.9 



Table 8.1: OPAL combined r polarization results. 



{Pr)= -Ar= (-14.9 ± 1.9 ± 1.3)% 

iFB _ _3 

Following these results we obtain: 



Apof= -jAe= (-8.9 ±2.2 ±0.9)%. (8.2) 



Ar= 0.153 ±0.023 (8.3) 

Ae= = 0.122 ±0.032 



where the errors include both statistical and systematic uncertainties. These 
are already corrected, (using ZFITTER [^), for the contribution of the 
intermediate 7, 7-Z° interference and other photonic radiative corrections. 
Therefore they can be interpreted in terms of the vector and axial-vector 
couplings within the improved Born approximation. In the context of the 
SM Ar and Ae yield a measured values for the coupling ratios of 

Vr/ttr = 0.077 ±0.012 (8.4) 
Ve/ae = 0.062 ±0.016 



where as before the error includes both statistical and systematic contribu- 
tions. 
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These measured values of {Pr) and A^^f correspond to the following values 
of sin^ (see Eq. 



sin"^ 9 yf{{Pr)) = 0.2308 ±0.0030 (8.5) 
sin2^^(A^„f) = 0.2345 ±0.0040. 

These values are in a good agreement with the hypothesis of lepton univer- 
sality. Therefore they can be averaged (with = 0.56) to give , 

v/a = 0.072 ±0.010 (8.6) 
sin^^^ = 0.2331 ±0.0023 



Fig. |8.1| presents the result of various determinations of sin^ ^^measured 



at LEP and SLD. These were calculated from the observed forward backward 
asymmetries of leptons {Afb{^^)) and quark pairs {ApBipb), Afb{cc) and 
(Qfb)- the inclusive forward backward asymmetry of negative quarks), the 
tau polarization and its forward backward asymmetry. 

The figures of the charge asymmetries are those which were presented as 



LEP results in the Moriond 94 conference The results given for the {Pr) 
and Ap^f are the modified OPAL numbers using the results of the present 
work. The obtained LEP average (with x^/D.O.F. = 5.9/5) is 

sin^ e^{LEP) = 0.2322 ± 0.0004 (8.7) 

For a comparison, also the SLD value of sin^ O^j} derived from their measure- 
ment of the left-right asymmetry, A^ji is given. This asymmetry is defined 
in a fully polarized electrons beam as 

A,n = = A (8.8) 

which means that it probes the same couplings as the A^^f . One can see 



from Fig. 3.1 that while there is a good agreement betwen all LEP results, 
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the present SLD value for sin^ 6^/ is shifted from the LEP numbers. 
Therefore as expected, the obtained when combining the measured weak 
mixing angle of LEP and SLD there is no good agreement [x^ / D.O.F. = 
13.3/6). The average value being 

sin^ e^{LEP + SLD) = 0.2317 ± 0.0003 (8.9) 
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LEP & SLD Results 



OPAL <Pt> 
OPAL A""^, 



pol 



LEP <Px> 



FB 



LEP A 

<Qfb> 

A''^(CC) 

Ai'^Cbb) 
a''^(II) 



pol 



-•—I 



0.2308+0.0030 
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Figure 8.1: Comparison of several determinations of the weak mixing angle 
from asymmetries as was presented in the MORIOND 94 conference. The 
values of LEP were modified using the results of the present analysis. Com- 
bined results of the LEP experiments and a combined LEP-SLD results are 
obtained as weighted averages assuming no correlations. 
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